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Abstract 


The  diffraction  of  shock  waves  (2  < Ms  < 8)  in  perfect  and 
imperfect  nitrogen  and  argon  by  sharp  compressive  corners  (2°  < Qy  < 60*) 
were  investigated  analytically  and  experimentally.  It  is  shown  that 
seven  shock-wave  diffraction  domains  exist  in  nitrogen  and  six  in  argon 
in  the  ranges  1 < Ms  < 10  and  0 < 0w  < 90°.  The  domains  consist  of  the 
four  well-known  shock  wave  reflections,  i.e.,  regular  reflection  (RR), 
single-Mach  (SMR),  complex-Mach  (CMS)  and  double-Mach  (DMR)  reflections. 
All  the  transition  boundaries  between  these  regions  were  established 
analytically  and  substantiated  by  the  present  experimental  results  as 
well  as  the  data  from  other  sources.  Over  100  experiments  were  conducted 
in  the  ITT  IAS  10  x 18  cm  Hypervelocity  Shock  Tube  equipped  with  a 23- cm 
dia  field  of  view  Mach-Zehnder  interferometer  equipped  with  a dual- 
frequency  laser  light  source . It  is  shown  that  real-gas  effects  have  a 
significant  influence  on  the  size  of  the  regions  and  their  transition 
boundaries.  Some  comparison  between  steady  and  nonstationary  reflections 
are  made  and  discussed.  Isopycnics  (lines  of  constant  density)  as  well 
as  density  distributions  along  the  wedge  surface  are  presented  for  the 
various  diffraction  processes  and  their  differences  and  similarities  are 
discussed. 
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Notation 


Ar 

E(a) 

He 

v. 

"i 

L 


Li 


^c 

% 


n 


P(2) 


r 


Ti 


t 


X 


y 


argon 

absolute  error  in  measuring  quantity  a 
helium 

specific  enthalpy  of  the  flow  in  region  (i) 

■ VLs 

distance  between  first  an4  second  triple  points 

distance  travelled  by  incident  shock  wave  from  the  wedge  corner 

flew  Mach  number  in  region  (i) 

incident  shock  wave  Mach  number 

change  over  Mach  number 

flow  Mach  number  in  region  (2)  for  a DMR  to  f arm 

nitrogen 

isopycnic  number 

oxygen 

flow  pressure  in  region  (i) 
partition  function 

distance  along  wedge  surface  or  shock  tube  wall  measured  from 

the  wedge  corner 

flow  temperature  in  region  (i) 

time 

flow  velocity  in  region  (i) 
incident  shock  wave  velocity 
velocity  vector  (i)  wrt  (j) 
coordinate 
coordinate 
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specific  heat  ratio 


time  for  incident  shock  wave  to  travel  distance  L 

s 

deflection  of  flow  from  its  original  direction  while  entering 
region  (i) 

maximum- deflect ion  angle 
sonic  angle 
actual  wedge  angle 
= effective  wedge  angle 

= sin  ^(l/M),  Mach  angle 
flow  density  in  region  (i) 

incident  angle  between  flow  in  region  (i)  and  oblique  shock  wave 
first  triple-point-traj ectory  angle 
second  triple-point- trajectory  angle 
angle  between  I and  R 


flow  properties  as  measured  from  the  second  triple  point 
see  Fig.  50 


Subscripts 

(o) 

flow  ahead  of  I for  RR,  SMR,  CMR  and  DMR 

(1) 

flow  behind  I for  RR,  SMR, 

CMR  and  DMR 

(2) 

flow  behind  R for  RR,  SMR, 

CMR  and  DMR 

(2’) 

flow  behind  I in  laboratory  frame  of  reference 

(3) 

flow  behind  M for  SMR,  CMR 

and  DMR 

00 

flow  behind  M^  for  DMR 

(5) 

flow  behind  R^  for  DMR 

Abbreviations 


complex-Mach  reflection 
double-Mach  reflection 


flow  deflection 

incident  shock  wave  in  RR,  SMR,  CMR  and  DMR 
kink 

left-hand  side 

Mach  stem  in  SMR,  CMR  and  DMR 
second  Mach  stem  in  DMR 

reflection  point  in  RR  or  location  of  M along  the  wedge  in 
SMR,  CMR  and  DMR 
wedge  corner 

reflected  shock  wave  in  RR,  SMF.,  C>'R  and  DMR  wrt  T 

reflected  shock  wave  in  CMR  and  DMR  with  respect  to  K or 

second  reflected  shock  wave  in  DMR 

right-hand  side 

regular  reflection 

slipstream  in  SMR,  CMR  arid  DMR 

second  slipstream  in  DMR 

shock  (wave)  diffraction 

single-Mach  reflection 

shock  (wave) reflection 

first  triple  point  in  SMR,  CMR  and  DMR 

second  triple  point  in  DMR 
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1 . INTRODUCTION 


When  a planar  moving  incident  shock  wave  encounters  a sharp  compressive 
corner  in  a shock  tube,  two  processes  take  place  simultaneously.  The  incident 
shock  wave  is  reflected  by  the  wedge  surface,  whereas  the  induced  nonstationary 
flow  behind  it  is  deflected  by  the  wedge  corner.  Throughout  this  report  the 
first  process  will  be  referred  to  as  shock-wave  reflection  (SR),  the  second  as 
flow  deflection  (FD)  and  the  overall  phenomena  as  shock-wave  diffraction  (SD). 

Shock-wave  diffraction,  i.e.,  the  combination  of  shock-wave  reflection 
and  flow  deflection,  for  a given  gas,  depends  on  three  factors: 

(1)  Mach  number  of  the  incident  shock  wave  M . 

s 

(2)  Corner  wedge  angle  e . 

(3)  Initial  thermodynamic  state  of  the  gas,  i.e.,  temperature  T0  and 
pressure  P0  (for  a perfect  gas  this  is  not  required) . 

Since  the  nonuniform  terminology  of  the  various  reflection  processes 
used  in  the  literature  tends  to  be  confusing  the  following  terms  are  suggested 
and  will  be  used  throughout  this  report: 

Regular  reflection  (RR),  Fig.  la 

Single-Mach  reflection  ( SMR) , Fig . lb 

Complex-Mach  reflection  (CMR),  Fig.  lc 

Double-Mach  reflection  (DMR),  Fig.  Id 

Although  RR  and  SMR  were  first  noticed  by  the  distinguished  physicist 
and  philosopher  Mach  (Ref.  1)  as  early  as  1878,  almost  no  work  was  done  in  this 
field  until  the  1940' s,  when  Neumann  reinitiated  the  problem  (Refs.  2 and  4) . 

An  intensive  investigation  at  Princeton  University  under  the  supervision  of 
Prof.  Bleakney  finally  led  to  the  discovery  of  CMR  by  Smith  (Ref.  6)  and  DMR 
by  White  (Refs.  17  and  18) . Once  these  four  different  types  of  reflection  were 
found  it  became  necessary  to  establish  the  transition  criteria  between  them. 

The  transition  problem  from  RR  to  SMR  was  first  studied  by  Neumann  (Refs.  2 and 
4)  who  assumed  the  following: 

(1)  A perfect  gas, 

(2)  Two-dimensional  inviscid  flows. 

(3)  When  there  are  two  possible  solutions  (the  so-called  weak  and  strong 
families  of  shock  waves  for  producing  a required  deflection)  the  weak- 
shock  solution  will  occur.  This  is  an  experimental  fact,  which  was 
verified  in  the  study  of  two-dimensional  supersonic  wedge  flows  (Ref.  28). 
Although  an  explanation  of  minimum  entropy  for  weaker  shocks  is  sometimes 
advanced  it  has  not  been  proven  analytically. 

(4)  The  flow  problem  is  pseudo- stationary.  This  arose  from  the  observation 
that  any  point  on  the  wave  configuration  having  a plane  radius  vector  r 


1 


■ v 


with  the  corner  as  origin  was  transformed  to  a new  point  cr,  where  c is 
a scalar  constant.  This  means  that  instead  of  three  independent  variables 
x,  y and  t the  phenomenon  was  describable  in  terms  of  x/t  and  y/t,  i.e., 
x and  y may  be  measured  relative  to  any  point  moving  with  constant  velocity 
with  respect  to  the  corner;  in  other  words,  the  flow  is  self-similar. 

Using  the  above  assumptions  Neumann  argued  (Refs.  2 and  4)  that  in 
an  RR  the  streamline  deflection  angle  02  through  the  reflected  shock  wave  is 
equal  in  magnitude  but  opposite  in  sign  to  the  deflection  angle  0p  through 
the  incident  shock  wave,  i.e.,  01  + 92  = 0 (Fig.  7).  This  is  violated  when 
the  wedge  angle  9W  decreases  to  the  point  where  it  forces  9q  to  exceed  in 
magnitude  the  maximum  deflection  angle  02m.  This  criterion  will  be  referred 
to  as  the  "detachment"  criterion.  The  detachment  criterion  can  best  be 
illustrated  by  using  the  pressure- deflection  (P,  9)- shock  polars.  Consider 
Fig.  2a  where  the  I and  R-polars  represent  the  incident  and  reflected  shock 
waves,  respectively.  Since  the  net  deflection  through  an  RR  is  zero  the 
solution  is  at  the  point  where  the  R-polar  intersects  the  P/P0-axis,  i.e., 
state  (2')  on  Rui-  As  the  wedge  angle  0W  decreases  the  R-polar  moves  away 
from  the  P/P0-axis  until  it  becomes  tangent  to  it  [Fig.  2a,  state  (2')  on 
Riv] • Upon  further  decrease  in  9>,  the  R-polar  will  not  intersect  the  P/P0- 
axis  any  more  and  an  RR  is  not  possible.  Consequently,  the  detachment 
criterion  is  represented  by  the  Rqv  shock  polar.  (Note  that  the  term  detach- 
ment comes  from  steady  flows  where  the  oblique  shock  detaches  at  this  angle.) 

Some  disagreement  between  the  detachment  criterion  and  experiments 
were  found  by  Smith  (Ref.  6) . In  his  experiments  RR  persisted  beyond  the  limit 
determined  by  the  detachment  criterion.  Bleakney  and  his  students  (Refs.  9,  10, 
13,  15,  17,  18,  20,  22,  23)  tried  to  resolve  the  disagreement  found  by  Smith 
(Ref.  6).  Unfortunately  however,  they  were  unable  to  do  so.  Kawamura  and 
Saito  (Ref.  24)  who  also  tried  to  resolve  this  problem  by  making  use  of  shock 
polars  discovered  that  the  point  of  tangency  between  the  R-polar  and  the  P/PQ- 
axis  (i.e.,  R^y,  Fig.  2a)  can  lie  outside  or  inside  the  I-polar  depending  on 
whether  or  not  the  value  of  Mo  is  greater  or  less  than  a certain  change-over 
value  Moc-  Unlike  the  case  of  a diatomic  gas  where  different  values  of  Moc 
are  reported  by  different  investigators  (Refs.  24,  60  and  8l),  only  one  value 
is  reported  for  a monatomic  gas  (Ref.  S'J)  . All  these  values  are  for  perfect 
gases  only. 

Henderson  and  Lozzi  (Ref.  60)  investigated  the  RR  -> SMR-transition 
problem  experimentally  in  a wind  tunnel  and  in  a shock  tube.  They  introduced 
an  alternative  criterion  which  has  the  property  that  the  system  always  remained 
in  mechanical  equilibrium  (i.e.,  no  pressure  discontinuities)  during  transition. 
Consider  Fig.  2a  and  note  that, once  RR  terminates  and  SMR  forms,  the  solution 
moves  from  the  point  where  the  R-polar  intersects  the  P-axis  (state  2'  on  Rqv) 
to  the  point  where  the  I and  R-polars  intersect  (states  2 and  3 on  Riv) • Con- 
sequently, a sharp  pressure  change  (from  P2'  to  P2)  is  associated  with  this 
transition,  if  the  detachment  criterion  is  accepted.  Henderson  and  Lozzi 
argue  (Ref.  60)  that  "a  system  which  develops  a pressure  discontinuity  during 
transition  cannot  be  in  mechanical  equilibrium".  Furthermore  they  say  that, 

"If  a pressure  discontinuity  occurs  during  transition  then  an  unsteady  wave  of 
finite  amplitude  or  a finite  amplitude  band  of  waves  will  be  generated  in  the 
flow.  These  would  be  expansion  [waves]  for  RR  -»SMR  and  compression  [waves^ 
for  SMR  ->  RR . " However,  since  these  waves  have  not  been  observed,  they 
discarded  (Ref.  60)  the  detachment  criterion  and  suggested  an  alternative 
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criterion  that  enables  the  system  to  be  in  mechanical  equilibrium  during 
transition.  In  order  to  maintain  the  system  in  mechanical  equilibrium,  the 
transition  should  take  place  at  the  point  where  the  R-polar  intersects  the 
I-polar  (SMR  solution)  on  the  P/P0-axis  (ER  solution) . This  is  illustrated 
by  states  (2)  and  (3)  on  the  Rn  polar.  The  formulation  of  this  criterion 
yields  ©i  + ©2  = e3  = 0 811(1  ^ be  referred  to  as  the  "me-.hanical-equili- 

brium"  criterion.  Consider  polar  Riii  (Fig.  2a)  and  note  that  according  to 
the  detachment  criterion  an  RR  takes  place  at  (2')  since  the  termination 
case  Riv  was  not  reached.  However,  according  to  the  mechanical-equilibrium 
criterion,  RR  cannot  occur,  since  the  termination  criterion  given  by  Rii  was 
exceeded.  It  is  worth  noting  that  the  area  of  disagreement  in  the  (Mg,  em- 
plane between  these  two  different  criteria  (Fig.  2b)  is  very  large.  The 
mechanical-equilibrium  criterion,  rather  than  reducing  the  previously- 
mentioned  disagreement  between  theory  and  experiment  found  by  Smith  (Ref.  6), 
where  RR  occurred  even  below  the  line  of  the  detachment  criterion  in  shock- 
tube  experiments  at  low  incident  shock  wave  Mach  numbers,  made  it  even  worse, 
for  their  line  lies  above  the  detachment- criterion  line  for  all  Ms  > 1.68 
(Mo  = 2.23,  0O  = 41.20).  Nevertheless,  Henderson  and  Lozzi  should  be  credited 
for  their  new  physical  approach  to  the  problem.  It  is  worth  noting  that  in 
an  experiment  where  the  wedge  angle  is  changed  gradually  perhaps  one  might 
obtain  two  different  criteria  for  RR  -»SMR  and  SMR  -»RR  transitions.  Consider 
Fig.  2a  and  note  that  if  one  starts  with  a given  SMR  at  states  (2)  and  (3)  on 
Riv  and  the  wedge  angle  is  increased  slowly,  it  is  possible  that  states  a (SMR), 
b (SMR),  c (SMR  -» RR)  and  d (RR)  might  be  encountered,  and  hence  the  transition 
would  follow  the  mechanical-equilibrium  criterion.  However,  if  one  started 
with  a given  RR  state  (2)  on  Ri  arid  then  decreased  fa  gradually,  it  is  possible 
that  the  sequence  of  events  might  be:  states  d (RR),  c (RR),  e (RR),  f (RR) 
and  a (SMR).  This  sequence  of  events  follows  the  detachment  criterion.  Such 
experiments  have  not  been  made  to  date  and  need  further  study. 

During  their  attempt  to  substantiate  their  mechanical-equilibrium 
criterion,  Henderson  and  Lozzi  (Ref.  60)  found  a "remarkable  anomaly"  between 
their  results  from  wind-tunnel  and  shock-tube  experiments.  In  their  shock- 
tube  results  RR  continued  to  exist  below  the  detachment  and  mechanical -equili- 
brium transition  boundary  lines  (Fig.  2b)  in  a region  where  the  perfect-gas 
theory  had  no  RR  solution.  They  resolved  the  anomaly  by  advancing  that  those 
RR- configurations  found  below  their  mechanical-equilibrium  transition  line 
were  really  undeveloped  DMR- configurations  in  which  all  shock  waves,  slip- 
streams and  triple  points  typical  of  a well-developed  DMR  were  too  close  to- 
gether to  be  observed. 

Hornung  and  Kychakoff  (Ref.  65)  initiated  another  criterion  for  the 
termination  of  RR.  They  argued  that  in  order  for  an  SMR  to  form,  a length 
scale  must  be  available  at  the  reflection  point,  i.e.,  pressure  signals  must 
be  communicated  to  the  reflection  point.  This  single  argument  eventually  led 
them  to  two  different  termination  lines  for  RR,  depending  on  whether  the 
flow  under  consideration  is  steady  or  nonstationary. 

Consider  the  nonstationary  RR  in  Fig.  2c  (i)  and  note  that  the  length 
fa  can  affect  the  reflection  point  P jnly  when  a subsonic  flow  is  established 
between  Q and  P (in  a frame  of  refer  nee  attached  to  P) . In  a steady  flow 
(Fig.  2c,  ii)  the  geometrical  length  fa  can  affect  the  reflection  point  P only 
if  a propagation  path  exists  between  point  Q and  point  P via  the  expansion 
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wave  at  Q' . This  is  possible  only  if  the  flow  between  P and  Q'  is  subsonic. 
According  to  Hornung  and  Kychakoff  (Ref.  65)  this  could  happen  if  an  SMR 
existed  since  the  flow  behind  the  Mach  stem  is  always  subsonic.  Consequently, 
they  argued  that  transition  takes  place  the  very  first  time  when  an  SMR  can 
occur.  Consider  Fig.  2a  and  note  that  this  corresponds  to  states  (2)  and  (3) 
on  Rii-  That  also  represents  the  mechanical-equilibrium  criterion  of  Henderson 
and  Lozzi  (Ref.  60) . 

Consequently,  Hornung  and  Kychakoff' s sonic  criterion  led  to  two 
different  transition  lines  in  steady  and  nonstationary  flows.  In  steady  flows 
the  transition  satisfies  the  mechanical-equilibrium  criterion,  while  for  non- 
stationary flows  their  analysis  led  to  a new  transition  line,  which  will  be 
referred  to  as  the  "sonic"  criterion.  Note  that  if  the  sonic-transition  line 
was  drawn  in  Fig.  2b,  it  would  coincide  with  the  detachment  criterion- transition 
line  since  it  lies  below  it  only  at  very  low  incident  shock  waves  (Ref.  74) ■ 
Consequently  it  could  not  be  drawn  on  Fig.  2b.  Note  that  Hornung  and  Kychakoff 
(Ref.  65)  claim  to  have  experimental  data  obtained  in  both  shock-tube  and  wind- 
tunnel  flows  which  verify  their  criterion. 

The  transition  boundary  lines  in  Fig.  2b  were  calculated  using  the 
steady-flow  theory  for  regular  and  single-Mach  reflection.  This  is  justified 
by  the  fourth  assumption  discussed  previously.  In  1951  Jones,  Martin  and 
Thornhill  (Ref  l6)  used  this  assumption  to  transform  the  equations  of  motion 
so  that  instead  of  three  independent  variables  x,  y,  t only  two  x/t  and  y/t 
were  sufficient  to  describe  the  phenomenon.  Their  analysis  was  later  verified 
experimentally  by  Parks  (Ref.  19).  However,  since  in  the  nonstationary  case 
the  shock  waves  are  usually  curved,  the  results  from  the  steady- state  analysis 
are  correct  only  in  the  vicinity  of  the  reflection  point  for  RE  and  the  triple 
point  in  SMR,  CMR  and  DMR . 

Skews  (Refs.  57  and  58)  tried  to  resolve  the  disagreement  between 
analysis  and  experiment  by  arguing  that  the  boundary  conditions  for  RR  and  SMR 
should  be  relaxed.  In  the  case  of  SMR  he  assumed  that  the  pressure  and  flow 
deflection  immediately  behind  the  reflected  shock  wave  R need  not  be  the  same 
as  those  behind  the  Mach  stem.  Consequently  he  defined  two  relaxed  variables 
F = (P2  - Po)/(P3  - ?o)  and  e = 02  - 03-  (Note  that  in  the  usual  SMR  theory 
F = 1 and  e = 0.)  He  supported  this  new  approach  by  the  fact  that  some 
investigators  claimed  that  the  slipstream  was  not  a surface  of  discontinuity 
but  a shear- flow  mixing  region.  However,  since  excellent  agreement  was  found 
between  the  classical  (without  relaxed  boundary  conditions)  theory  and 
expe  iment  for  the  angles  between  the  various  shock  waves  and  slipstream 
(Refs.  6,  10,  13,  15,  18  and  43),  one  cannot  justify  the  need  to  relax  the 
boundary  conditions  for  SMR. 

In  the  RR-analysis  (Ref.  58)  Skews  relaxed  the  usual  boundary 
condition  that  the  flow  immediately  behind  the  reflected  shock  is  uniform  and 
parallel  to  the  wall.  Instead,  he  defined  a relaxation  angle  e between  the 
flow  direction  and  the  wall.  (Note  e = 0 is  the  usual  boundary  condition.) 
Although  his  numerical  results  were  in  better  agreement  with  experiments  one 
must  bear  in  mind  that  his  analysis  was  based  mainly  on  mathematical  concepts. 
Unlike  the  foregoing  SMR-case  physical  considerations  impose  the  boundary 
condition  that  the  supersonic  flow  behind  the  reflected  shock  wave  is  parallel 
to  the  rigid  wall . 

Once  RE  terminates,  three  different  types  of  reflection  can  occur 
in  nonstationary  flows,  i.e.,  SMR,  CMR  and  DMR.  White  (Ref.  17)  was  the 
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first  to  notice  that  when  the  flow  Mach  number  behind  the  reflected  shock  wave 
R,  becomes  supersonic  in  a frame  of  reference  attached  to  the  triple  point  a 
"kink"  K forms  in  R and  the  transition  SMR  -»  CMR  occurs.  A mechanism  for  the 
transition  was  later  suggested  by  Gvozdeva  et  al  (Ref.  45)  and  Henderson  and 
Lozzi  (Ref.  60) . Henderson  and  Lozzi  suggested  that  a "band  of  compression 
waves"  must  exist  in  a CMR.  These  compression  waves  then  converge  to  a shock 
wave  to  form  DMR  when  M2  > 1.  Unfortunately,  their  suggestion  was  not  substan- 
tiated analytically  or  experimentally.  In  addition  the  precise  value  of  for 
the  termination  of  CMR  and  the  formation  of  DMR  was  not  established  either. 

The  correct  gasdynamic  criterion  for  the  termination  of  CMR  and  formation  of 
DMR  was  established  during  the  present  study.  It  will  be  shown  in  Section  2.4.2 
that  the  flow  behind  the  reflected  shock  wave  must  be  supersonic  with  respect  to 
the  kink  of  a CMR  in  order  for  a DMR  to  form.  Consequently,  the  CMR  ^ DMR 
transition  occurs  at  ^2  = 1.  It  was  found  further  that  ^ = 1 corresponds  to 
M2  « 1.3  (M2  measured  wrt  the  first  triple  point)  and  hence  Gvozdeva  et  al 
(Ref.  45)  and  Henderson  and  Lozzi  (Ref.  60)  were  correct  in  predicting  the 
transition  at  M2  > 1- 

An  attempt  to  establish  some  transition-boundary  lines  experimentally 
was  made  by  Bazhenova  et  al  (Ref.  62).  However,  their  experiments  did  not  cover 
a significant  range  of  incident  shock  wave  Mach  number  (Mg)  and  corner  wedge 
angle  (ew).  Their  experimental  boundaries  for  SMR,  CMR  and  DMR  are  limited 
to  Mg  < 4.5.  Unfortunately,  experiments  made  by  other  researchers  also  did 
not  cover  a wide  range  of  interest.  Smith  (Ref.  6),  Kawamura  and  Saito  (Ref. 

24)  and  White's  (Ref.  17,  18)  experimental  data  covered  only  the  range  1 < 

Mg  < 2.75,  while  Henderson  and  Lozzi 's  (Ref.6o)  experiments  were  centred 
around  the  RR- termination  criterion  line.  Law  and  Glass  (Ref.  51)  were  the 
first  to  extend  the  range  of  incident  shock  Mach  number  up  to  Ms  < 8,  but 
their  corner  wedge-angle  range  at  the  lower  end  was  limited  (25°  < 0w  < 60*). 

Figure  3 which  is  a reproduction  of  Fig.  5 from  the  paper  by  Bazhenova 
et  al  (Ref.  62)  summarizes  all  the  theoretical  and  experimental  knowledge  (ex- 
cluding the  mechanical  equilibrium  criterion  for  the  termination  of  RR,  discussed 
earlier)  concerning  the  regions  and  boundaries  of  regular,  single-Mach,  complex- 
Mach  and  double-Mach  reflections  available  when  the  present  study  started.  Only 
the  termination  criterion  of  RR  is  calculated  theoretically  for  both  perfect 
and  inperfect  gases  (lines  1 and  2,  respectively).  Although  the  SMR  boundary 
line  is  also  calculated,  line  5 (for  a perfect  gas  only)  it  does  not  start  or 
terminate  at  any  other  boundary  line  and  hence  does  not  enclose  any  region. 

The  imperfect  boundary  lines  between  SMR,  CMR  and  DMR  (lines  3 and.  4)  were  all 
obtained  experimentally  and  they  do  not  encompass  any  closed  region  either. 

There  is  no  information  about  the  types  of  reflections  for  wedge  angles  in  the 
range  0 < ew  < 20°. 

Although  most  of  the  investigators  were  interested  in  finding  the 
correct  transition  criteria  from  one  type  of  reflection  to  another  there  were 
some  attempts  (Refs.  5,  15,  6l,  63,  64)  to  solve  the  entire  flow  field  for  a 
given  reflection  utilizing  various  analytical  or  numerical  methods.  Unfor- 
tunately the  numerical  results  suffer  from  a lack  of  agreement  for  the  same 
initial  conditions  depending  on  the  technique  used  for  a solution  (Refs.  6l, 

63  and  73) • 

In  view  of  the  above  literature  survey,  the  present  research  was 
directed  towards: 

(l)  Resolving  the  disagreements  concerning  the  termination  criterion  of  RR. 
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(2)  Establishing  the  correct  termination  criterion  of  CMR  and  hence  the 
formation  criterion  of  DMR. 

(3)  Defining  domains  and  transition  boundaries  of  the  various  reflections 
in  order  that  they  can  be  predicted  a priori . 

(4)  Extending  the  experimental  data  over  a much  wider  range,  i.e.,  0 < Qy  < 60° 
and  1 < Mg  < 10. 

(5)  Comparing  the  existing  numerical  predictions  of  the  isopycnic  (constant- 
density  lines)  in  the  flow  field  with  those  obtained  interferametrically. 

(6)  Resolving  and  clarifying  areas  of  disagreement  existing  in  the  literature 
concerning  this  problem. 

To  achieve  these  goals  it  was  decided  to  approach  the  problem  using  the 
following  physical  concepts: 

(1)  It  was  assumed  that  the  physical  processes  are  governed  by  the  usual 
equations  of  conservation  of  mass,  momentum  and  energy,  as  written  for 
oblique  shock  waves.  The  pressure  and  flow  directions  were  taken  as 
invariant  across  shear  layers,  and  very  near  the  wall  flow  was  assumed 
to  be  parallel  to  the  wall. 

(2)  The  shock-reflection  process  depends  on  the  flow- deflection  process  behind 
it.  Consequently,  the  interaction  between  these  two  processes  should  be 
investigated  and  understood.  Note  that  this  approach  was  initiated  by  Law 
(Ref.  48)  and  Law  and  Glass  (Ref.  51)  • Unfortunately,  it  was  not  adopted 
by  other  investigators. 

(3)  Real-gas  effects  strongly  influence  the  location  of  the  transition  boundary 
lines  in  the  (Mg,  q^) -plane.  Consequently,  they  should  be  included  in 
solving  the  equations  of  motion.  This  is  substantiated  by  lines  1 and-'ig 
of  Fig.  3 that  illustrate  the  significance  of  real-gas  effects  in  shifting 
the  RR- transition  boundary  line. 

Due  to  the  difficulties  in  solving  the  nonlinear  equations  of  motion 
it  was  decided  to  restrict  the  inperfect-gas  model  to  thermodynamic  and  chemical 
equilibrium. 

The  inperfect-gas  model  in  this  report  takes  into  account  the  excita- 
tion of  the  internal  degrees  of  freedom  such  as  rotation-vibration  coupling, 
vibration  and  dissociation  in  nitrogen,  and  electronic  excitation  and  ionization 
in  argon  (see  Appendix  A in  Ref.  48  for  details). 

Note  that  the  flow  usually  passes  through  two  shock  waves  (incident 
and  reflected) . Consequently,  even  if  the  incident  shock  wave  is  very  weak 
the  temperature  behind  the  reflected  shock  wave  can  be  such  that  the  contri- 
bution of  the  vibrational  degree  of  freedom  in  a diatomic  gas  can  no  longer 
be  neglected. 

(4)  The  CMR  and  SMR  should  be  treated  from  a frame  of  reference  attached  to 
the  kink  of  the  CMR  or  the  second  triple  point  of  the  DMR.  Therefore, 
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the  correct  transformation  of  the  frame  of  reference  from  the  first 
triple  point  should  be  investigated,  understood  and  formulated. 
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During  the  present  study  all  the  criteria  for  the  formation  and 
termination  of  RR,  SMR.  CMR  and  DMR  have  been  established  analytically.  Con- 
sequently, the  (Mg,  0^;  and  the  (Ms,  Qy) -planes  were  divided  into  the  domains 
of  the  different  types  of  reflection  and  diffraction  processes,  respectively. 

Over  100  experiments  have  been  performed  in  the  10  cm  x l8  cm  IfTIAS 
Hypervelocity  Shock  Tube  in  nitrogen  and  argon,  at  an  initial  temperature  of 
nearly  300  K and  a pressure  of  15  torr.  The  shock-Mach-number  range  was 
2 < Mg  < 8 over  a series  of  wedge  angles  2°  < < 60°.  Dual-wavelength 

laser  interferograms  were  obtained  by  using  a 23- cm  dia  field  of  view  Mach- 
Zehnder  interferometer.  For  each  and  every  type  of  diffraction  the  shock 
shapes  and  density  field  (isopycnics)  as  well  as  the  density  distribution 
along  the  wedge  surface  were  deduced  from  the  corresponding  interferograms, 
and  comparisons  were  made  with  existing  numerical  analyses. 

The  experimental  results  from  the  present  experiments  and  other 
sources  substantiate  the  present  analysis.  In  addition  areas  of  disagreement 
which  existed  in  the  literature  have  been  clarified  and  resolved.  Finally, 
some  comparisons  between  steady  and  nonstationary  shock  reflections  were  made 
and  discussed. 


2.  OBLIQUE  SHOCK-WAVE  REFLECTION  IN  NONSTATIONARY  FLOWS 

Since  the  incident  shock  wave  moves  with  a constant  velocity,  the 
entire  problem  can  be  considered  from  a pseudo- stationary  point  of  view,  by 
attaching  a frame  of  reference  to  any  point  on  the  incident  shock  wave.  The 
above  suggests  that  instead  of  three  independent  variables  x,  y and  t the 
phenomenon  is  now  describable  in  terms  of  x/t  and  y/t  and  the  flow  problem 
is  self- similar . 


2.1  Reasons  for  Reflection 

Consider  now  a planar  incident  shock  wave  I,  having  Mach  number 
Mg  (Fig.  4a),  colliding  with  a sharp  compressive  corner  of  angle  ew.  Denote 
the  states  ahead  of  and  behind  it  by  (0)  and  (l),  respectively.  Attach  a frame 
of  reference  to  the  reflection  point  P where  I meets  the  wedge  surface.  The 
flow  in  state  (0)  moving  parallel  to  the  wedge  surface  approaches  I with  a 
velocity  UQ  = Us  cosec  or  a Mach  number  Mq  = Ms  cosec  0o,  where  0o  = 

90°  - ©w  is  the  incident  angle.  While  passing  through  I the  flow  is  deflected 
towards  the  wedge  surface  by  an  angle  from  its  original  direction,  and  its 
dynamic  and  thermodynamic  properties  are  changed.  The  deflection  causes  the 
flow  in  state  (l)  to  approach  the  wedge  surface  obliquely  at  an  angle  91 
(Fig.  4a). 
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From  single  gasdynamic  considerations  this  situation  is  analogous 
to  the  case  of  having  a steady  flow  (Pi,  Tp  and  Mp)  over  a wedge  9p  (Fig.  5). 

It  is  well  known  that  depending  on  Mp  this  flow  can  negotiate  the  corner 
either  through  a subsonic  turning  (Fig.  5a)  if  Mp  < 1,  or  through  a bow  shock 
wave,  straight  and  attached  (Fig.  5b),  curved  and  attached  (Fig.  5c)  or  curved 
and  detached  (Fig.  5d)  depending  on  Bp  (perfect  gas),  when  the  flow  is  super- 
sonic Mp  > 1 (see  Appendix  A for  further  details).  For  an  inperfect  gas  the 
initial  temperature  Tp  and  pressure  Pp  are  also  inportant . 

Applying  the  above  steady  flow  results,  the  pseudo-stationary 
situation  suggests  that  if  Mp  > 1 the  flow  can  negotiate  the  wedge  surface  either 
through  an  attached  (to  point  P)  shock  wave  (straight  or  curved,  Fig.  6b  and 
6c,  respectively)  which  will  result  in  an  RE,  or  through  a detached  shock 
wave  that  will  result  in  a reflection  through  a three- shock  confluence 
configuration  rather  than  two  (Fig.  6d) . In  the  case  of  Mp  < 1 the  analogy 
suggests  that  the  flow  adjusts  its  direction  through  a subsonic  turning 
(Fig.  6a) . However,  it  has  been  found  experimentally  throughout  this  study 
that  in  the  nonstationary  case  SMR  occurs  even  for  combinations  of  Ms  and  9w 
for  which  Mp  < 1 wrt  point  P.  We  suggest  that  the  reflection  (when  Mp  < 1 
arises  from  the  interaction  between  the  incident  shock  wave  and  the  bow  shock 
generated  by  the  wedge  (Fig.  4b).  Let  the  frame  of  reference  be  attached  to 
point  T,  where  these  two  shock  waves  meet.  State  (2)  results  from  state  (l) 
on  passing  through  a bow  shock  wave  R.  Consequently  P2  > Pi  and  P2/P0  > Pp/Po • 

To  satisfy  the  last  condition  the  portion  of  the  incident  shock  wave  I that 
lies  below  T (dashed  line  TP)  must  move  forward  from  P to  M in  order  to 
become  more  normal  to  the  oncoming  relative  flow,  forming  a Mach  stem  TM, 
triple  point  T and  slipstream  S (omitted  for  clarity) . 

The  triple  point  trajectory  path  angle  x plays  a significant  role 
when  it  makes  9w  = 8w  + X large  enough  such  that  Mp  will  always  be  greater 
than  unity  wrt  triple  point  T (a  necessary  condition  for  the  existence  of  R) . 

Note  that  when  the  induced  flow  behind  the  incident  shock  wave 
becomes  subsonic  (R£>>  < 1),  and  no  shock  waves  arise  from  the  corner,  since 
the  flow  can  turn  subsonically  (Fig.  5a),  neither  the  first  reason,  Fig.  4a 
(the  only  one  quoted  in  the  literature),  nor  the  one  added  here,  Fig.  4b, 
apply. 


Perhaps  the  most  significant  conclusion  from  the  above  discussion 
is  that  shock-wave  reflection  in  the  nonstationary  case  depends  additionally 
on  the  flow-deflection  process  over  the  wedge  corner,  and  cannot  be  treated 
independently  of  it . 

2.2  Regular  Reflection  (RR) 

The  regular  reflection  shock-wave  configuration  (two- shock  con- 
fluence) is  shown  in  Fig.  7-  The  frame  of  reference  is  attached  to  the 
reflection  point  P and  it  is  moving  parallel  to  the  wedge  surface  with  the 
constant  velocity  UQ  = Us  cosec  0O,  where  0o=9°o-  0w  is  the  angle  of 
incidence  between  the  incident  shock  wave  I and  the  oncoming  flow  UQ-  Upon 
passing  through  I the  flow  (Uo)  is  deflected  by  an  angle  0p  from  its  original 
direction.  This  deflected  supersonic  flow  Up  is  then  redeflected  by  the 
reflected  shock  wave  R,  by  an  angle  02  = -0p}  so  that  the  flow  is  again 
parallel  to  the  wedge  surface,  meeting  the  required  boundary  condition. 
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It  has  already  been  mentioned  (Section  2.1)  that  the  flow  behind  R 
(U2)  could  be  either  supersonic  (Mg  > l)  or  subsonic  (M2  < 1)  wrt  point  P. 
Consequently,  two  different  types  of  RR  are  possible  since  the  reflected 
shock  wave  can  be  either  straight  or  curved  (Figs.  6b  and  6c,  respectively) 
in  the  vicinity  of  the  reflection  point  P. 

2.2.1  Analytical  Formulation 

In  the  frame  of  reference  attached  to  the  reflection  point  P,  the 
incident  and  reflected  shock  waves  (I  and  R)  can  be  treated  using  steady-flow 
theory..:  ■ Considering  I and  R separately  and  making  use  of  the  oblique  shock 
wave  relations  the  equations  of  motion  for  the  regular  reflection  are; 

For  I: 


pota*0o  = pitan^o  ■ ®i)  (2.1) 

poUosin0o  = piUisin^o  " 9l)  ( 2- 2) 

Po  + poUo2sin\  = P1  + P1U12sin2(0o  - 0L)  (2.3) 

ho  + I Uo2sin\  = hl  + | U12sin2(0o  - 91)  (2.4) 

For  R; 


P1tan01 

= p2tan(01  - 

e2) 

(2.5) 

piuisin0i  = 

! P2U2sin(01 

- e2> 

(2.6) 

P1 

2 2 

+ p1U1  sin  01 

P2  + P2U2 

sin2(0L  - 

02) 

(2.7) 

hl 

■ 1 ..  2 . 2 . 

+ 5 Ui  sln 

- + 5 u22 

sin2(0^  - 

02) 

(2.8) 

The  boundary  condition 

is: 

°2  ' '61 

(2.9) 

If  equilibrium  is  assumed,  two  thermodynamic  properties  are  suffi- 
cient to  define  a state,  e.g.,  p = p(P,T),  h = h(P,T)*.  Consequently,  the 
above  nine  equations  have  thirteen  independent  variables,  namely:  PD,  Pj_, 
P2,  TOJ  Tlf  T2  > Uq,  U^,  Ug,  0q»  0i>  9j_  and  Q2.  Thus,  if  four  of  them  are 
known,  the  remaining  nine  can  be  calculated.  The  four  known  parameters  are 
usually  the  flow  pressure  and  temperature  ahead  of  the  incident  shock  wave 
Po  and  To>  respectively,  the  flow  velocity  Uo  = Us  sec Gw  and  its  angle  of 
incidence  with  the  primary  shock  wave  0O  =90*  - &w-  Note  that  in  shock 
tube  experiments  PQ,  T0,  Mg  and  Gw  are  the  known  or  measured  parameters. 
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2. 2. 1.1  Perfect-Gas  Solution 


When  the  gas  under  consideration  is  assumed  to  be  calorically  and 
thermally  perfect  the  above  nine  equations  can  be  reduced  to  the  following 
four  equations  (Ref.  31): 


tane 


1,2 


J^2 


- 1 


* <1  - V 


1 + Mo2 


^r^,i  +fri 


Ju2 

^i 


2-Z.l  + _±i2 


7 + 1 


0,1 


1/2 


(2.10,  11) 


zJlA  + _k 

7-1 


7 + 1 + ° 


- 1 


P 

_c 

p. 


(2.12) 


0i  + ^ “ 0 (2.13) 

The  above  set  of  equations  involves  seven  unknowns,  namely:  7,  Mq, 

Mq,  9q,  02,  Pq/P0  ^2/^1*  Consequently,  in  order  to  solve  this  set,  three 
of  the  variables  must  be  known.  If  7,  Mq,  Pl/PQ  are  selected,  the  rest  can 
be  easily  calculated,  since  the  solution  now  involves  solving  one  equation  at 
a time  with  one  unknown,  i.e.,  9q  is  calculated  from  Eq.  2.10,  then  Mi  and 
©2  from  Eqs.  2.12  and  2.13,  respectively,  and  finally  P2/P1  from  Eq.  2.11. 

This  also  suggests  that  Eqs.  2.10  to  2.13  are  reducible  to  one  equation  with 
one  unknown. 

It  should  be  mentioned  that  for  this  case  (perfect  gas),  the  equations 
of  motion  (Eqs.  2.1  to  2.9)  have  been  reduced  by  Henderson  (Ref.  31)  to  a single 
polynomial  of  order  six. 

2. 2. 1.2  Ingerfect-Gas  Solution 


When  real-gas  effects  such  as  vibration,  vibration-rotation  coupling 
and  dissociation  in  diatomic  gases  or  electronic  excitation  and  ionization  in 
monatomic  gases  are  considered,  the  equations  of  motion  (Eqs.  2.1  to  2.9) 
cannot  be  simplified.  Consequently,  one  has  to  solve  nine  nonlinear  algebraic 
equations  with  nine  unknowns.  However,  since  the  nine  equations  of  motion 
consist  of  two  similar  sets  of  four  equations  each  (2.1  to  2.4  and  2.5  to  2.8), 
a method  was  developed  in  which  each  set  is  treated  independently.  The  final 
solution  is  obtained  when  the  last  equation  (Eq.  2.9)  is  satisfied. 

Consider  the  first  set  (Eqs.  2.1  to  2.4)  which  involves  eight 
variables  Pq,  T0,  Uq  and  0O  on  the  LHS  and  Pq,  Tq,  Uq  and  9q  on  the  RHS . 


L 
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Since  the  variables  on  the  LHS  consist  of  the  initial  conditions  and  hence  are 
known,  the  RHS  can  be  calculated  and  pp,  Tp,  Up  and  0p  are  found. 

The  second  set  (Eqs.  2.5  to  2.8)  which  again  consist  of  eight 
variables  Pp,  Tp,  Up  and  0p  on  the  LHS  and  P2,  T 2 , l>2  and-  02  on  the  RHS  is 
slightly  different  since  not  all  the  parameters  on  the  LHS  are  known.  Only 
Pi,  Ti  and  Ui  which  were  obtained  from  the  solution  of  the  first  set  are 
known,  while  0i  is  still  unknown.  However,  if  we  guess  a value  of  0p,  the 
second  set  becomes  exactly  the  same  as  the  first  one  and  the  four  parameters 
on  the  RHS  (P2,  T2>  U2  and  02)  can  calculated.  The  correct  value  of  0p 
must  satisfy  the  boundary  condition,  i.e.,  Eq.  2.9.  The  Newton-Raphson  method 
ensuf«s  a rapid  convergence  to  the  correct  value  of  0p.  A computer  program 
(based  on  the  above  method),  as  well  as  a typical  output  is  given  in  Appendix  B. 

2.2.2  Graphical  Solution  (Shock-Polar  Presentation) 

A graphical  solution  of  an  RR  can  be  obtained  by  making  use  of 
(P,  e)-shock  polars  (Ref.  27) • Consider  Fig.  8 in  which  the  I-polar  corre- 
sponds to  POJ  T0  and  Mo  = Mg  sec^.  State  (l)  behind  I can  be  easily 
obtained  from  the  known  incident  angle  0o  = 90°  - 9w  Once  state  (1)  is 
determined  the  shock  polar  that  corresponds  to  Pp,  Tp  and  Mp,  i.e.,  the 
R-polar,  can  be  drawn.  The  solution,  state  (2),  is  at  the  point  where  the  net 
deflection  equals  zero  (-6  = 0),  i.e.,  where  the  R-polar  intersects  the  P/PQ- 
axis.  Figure  8 reveals  that  theoretically  two  solutions  are  possible,  (2)  and 
(2').  However  (2')  corresponds  to  the  so-called  strong  solution  and  is  dis- 
carded on  physical  grounds  that  it  is  not  observed  in  practice. 

It  is  worthwhile  mentioning  again  that  the  reflected  shock  wave  R 
is  straight  (in  the  vicinity  of  the  reflection  point  P)  if  the  solution  on 
the  R-polar  lies  in  the  interval  e - 0£»  (see  Fig.  9a)  whereas  for  the  small 
interval  0sjj  - 0mg  a curved  reflected  imock  wave  is  obtained  (Fig.  9b)  . An 
RR  with  a curved  reflected  shock  wave  can  be  seen  in  Fig.  20  of  Ref.  23. 

2.2.3  Termination  of  Regular  Reflection 

As  noted  earlier  (Section  1)  three  different  criteria  for  the 
termination  of  RR  exist  in  the  literature.  The  first  one  is  due  to  Neumann 
(Ref.  2)  and  is  known  as  the  "detachment"  criterion,  which  states  that  RR 
terminates  when  the  R-polar  becomes  tangent  to  the  P-axis  on  the  (P,  e)-plane 
(Fig.  2a,  Riv)  . Analytically  this  can  be  expressed  as: 


61  + ®2  = ° (2*1*0 

in 

The  second  criterion,  the  "mechanical  equilibrium"  criterion,  is 
due  to  Henderson  and  Lozzi  (Ref.  60) . It  states  that  RR  terminates  when  the 
I and  R-polars  intersect  on  the  P axis  (Fig.  2$,  Rpp)  . Analytically  it  assumes 
the  relation: 


01  + 02  = e3  = 0 (2.15) 
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The  third,  criterion  is  due  to  Hornung  and  Kychakoff  (Ref.  65)  which 
states  that  RR  terminates  when  the  R-polar  intersects  the  P-axis  at  its  sonic 
point , i . e . , 


ex  + e2  =0  (2.16) 

s 

Since  the  sonic  point  and  the  point  of  maximum  deflection  are  usually  very 
close,  it  is  impractical  to  distinguish  between  this  criterion  and  the 
detachment  criterion,  and  hence  only  the  first  two  criteria  will  be  considered 
here.  One  should  note  that  the  third  criterion,  the  so-called  sonic  criterion, 
lies  slightly  below  the  detachment  criterion  at  very  low  Mach  numbers. 

Hornung  and  Kychakoff  who  suggested  the  third  criterion  ( Eq.  2.16) 
found  that  two  different  criteria  for  the  terminations  of  RR  exist,  one  for 
steady  flow  and  one  for  nonstationary  flow.  They  claim  that  their  recent 
finding  was  substantiated  by  their  experiments  in  both  a wind  tunnel  and  a 
shock  tube . Their  data  for  nonstationary  flows  (Refs.  65  and  74)  do  not  agree 
with  the  "mechanical-equilibrium"  criterion.  One  should  also  keep  in  mind 
that  Henderson  and  Lozzi's  experimental  data  revealed  a "remarkable  anomaly" 
between  the  results  obtained  in  steady  and  nonstationary  (or  pseudo-stationary) 
flows,  where  RR  was  obtained  in  a region  where  the  "perfect-gas  theory"  (for 
both  t-, criteria)  has  no  RR  solution.  Henderson  and  Lozzi  (Ref.  60)  resolved 

this  disagreement  for  their  shock-tube  data  by  arguing  that  the  RR-like  con- 
figurations were  really  undeveloped  DMR- configurations,  where  both  the  triple 
points  as  well  as  the  slipstreams  and  the  shock  waves  typical  of  a well- 
developed  DMR  were  too  small  to  be  observed.  (They  estimate  the  primary  and 
secondary  confluences  as  being  0.1  mm  apart.)  A few  experiments  conducted  in 
the  present  work  in  order  to  verify  their  argument,  using  an  optical  magnifi- 
cation of  5-4  for  Mg  « 4.7,  9w  = 60°  and  a wedge  8.5  cm  long,  failed  to  show 
any  sign  of  a DMR. 

In  all  the  experiments  claimed  by  Henderson  and  Lozzi  to  be  undeveloped 
DMR-configurations,  the  reflected  shock  wave  R is  straight  near  the  reflection 
point.  This  indicates  that  the  flow  behind  R is  supersonic  (l^>  > 1)  and  hence 
parallel  to  the  wall.  However  a DMR  should  produce  subsonic  flow.  Consequently, 
all  the  available  evidence  favours  Hornung  and  Kychakoff' s (Ref.  65)  conclusion 
that  the  "mechanical  equilibrium"  criterion  established  by  Henderson  and 
Lozzi  (Ref.  60)  is  inapplicable  to  nonstationary  flows.  Therefore,  the  detach- 
ment criterion,  i.e.,  the  case  where  the  R-polar  becomes  tangent  to  the  P-axis, 
Eq.  2 .l4  is  the  correct  criterion  for  the  termination  of  RR. 

2. 2. 3.1  Change-Over  Mach  Number  - Mqc 

It  was  oted  that  Kawamura  and  Saito  (Ref.  24)  were  the  first  to  show 
that  the  point  of  ^angency  between  the  R-polar  and  the  P-axis  (the  termination 
criterion  of  RR)  can  lie  inside  or  outside  the  I-polar  (Figs.  10a  and  10b, 
respectively)  depending  on  whether  the  value  of  Mo  is  le  ss  or  greater  than  a 
certaJn  change-over  value  - Moc-  Unfortunately,  their  value  for  Mqc  = 3. 203 
(7=1.4)  does  not  match  their  other  variables  P0/Pi  = 0.433  and.  0o  = 4l.5°  which 
result  in  Moc  = 2.198.  Henderson  and  Lozzi  (Ref.  60)  pointed  out  that  a misprint 
occurred  and  that  it  should  read  Mqc  = 2.203,  which  is  in  good  agreement  with 
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the  above  values.  Unfortunately,  Henderson  and  Lozzi  (Ref.  60)  who  quoted  the 
value  Moc  = 2.23  in  the  text  of  their  paper  used  the  number  2.40  -in  the  caption 
of  their  Fig.  3>  thereby  adding  to  the  confusion  caused  by  Kawamura  and  Saito's 
misprint.  A different  value,  Moc  = 2.202,  was  recently  calculated  by  Molder 
(Ref.  8l)  for  a perfect  diatomic  gas. 

The  value  calculated  here  yields  Mqc  = 2.190  (Ms  = 1.450,  and  9^  = 48.55°) 
for  a perfect  diatomic  gas  (7  = 1.4)  and  Moc  = 2.185  (Ms  = 1.449  and  = 48.46°) 

for  imperfect  nitrogen  at  PQ  = 15  torr  and  T0  = 300  K.  The  significance  of  real- 
gas  effects  (vibration  and  vibration-rotation  coupling)  even  at  this  low  Mach 
number  (Ms  = 1.45)  is  clear,  albeit  small.  In  the  case  of  a perfect  monatomic 
gas  (7  = 5/3)  the  value  calculated  by  Molder  (Ref.  8l)  is  2.470  and  our  calculated 
value  of  Moc  is  2.453  (Ms  = 1.540  and  0^  = 51-13*)-  The  present  value  applies 
to  imperfect  argon  at  PQ  = 15  torr  and  T0  = 300  K. 

The  shock-polar  combinations  that  correspond  to  Moc  for  perfect  and 
imperfect  nitrogen  and  argon  (or  any  monatomic  gas)  are  shown  in  Figs.  11a,  lib 
and  11c,  respectively. 

When  RR  terminates  the  reflection  can  be  achieved  only  by  a three- 
shock  confluence,  i.e.,  a triple  point  (T) . In  nonstationary  flows  three  types 
of  reflection  can  occur  depending  on  the  flow  Mach  number  behind  the  reflected 
shock  wave  R.  One  more  difference  between  these  three  types  of  reflection 
and  RR  lies  in  the  fact  that  they  all  make  the  flow  near  the  wedge  surface 
(behind  the  Mach  stem)  subsonic  (M3  < l) , while  in  the  case  of  an  RR  the  flow 
near  the  wall  (behind  the  reflection  point)  is  usually  supersonic. 

2.3  Single-Mach  Reflection  (SMR) 

When  the  flow  behind  the  reflected  shock  wave  R is  subsonic,  i.e., 

M2  < 1,  in  a frame  of  reference  attached  to  the  triple  point  T,  it  can  negotiate 
the  wedge  surface  by  turning  subsonically . Consequently,  there  is  no  physical 
need  for  additional  shock  waves. 

The  wave  configuration  of  an  SMR  is  shown  in  Fig.  12.  The  frame  of 
reference  is  attached  to  the  triple  point  T and  it  moves  along  a straight  line 
(the  triple  point  trajectory  at  an  angle  *.)  with  a constant  velocity  U0  = Us 
cosec  0O  where  0O  = 90°  - 9w  ( 9w  = 0w  + X/  fc^e  an6le  of  incidence  between 
the  incident  shock  wave  I and  the  oncoming  flow  and  )(  is  the  angle  between  the 
triple  point  trajectory  and  the  wedge  surface.  Throughout  this  report  9w  will 
be  referred  to  as  the  effective  wedge  angle,  whereas  9^  is  the  actual  wedge 
angle.  The  Mach  stem  M lies  ahead  of  the  incident  shock  wave  I.  It  is  normal 
to  the  wedge  surface  but  not  necessarily  straight . 

The  flow  in  state  (0)  can  reach  the  region  bounded  by  R and  M by 
passing  through  two  shock  waves  I and  R or  only  one  shock  wave  M depending 
upon  whether  it  is  above  or  below  the  triple  point  trajectory  However,  gas- 
dynamic  considerations  imply  that  the  gas  must  be  compressed  to  the  same 
pressure,  and  must  move  in  the  same  direction  on  either  side  of  the  slipstream. 
Consequently,  all  other  thermodynamic  properties  such  as  entropy,  density, 
temperature , etc  , are  different  giving  rise  to  a slipstream  S.  It  divides 
the  two  regions  of  equal  pressure  and  flow  direction,  but  different  thermodynamic 
states . 
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2.3.1  Analytical  Formulation 


The  shock  waves  I,  R and  M are  again  treated  separately,  using  the 
oblique- shock- wave  relations  across  each  one  of  them.  Consequently,  the  equations 
of  motion  for  an  SMR  are: 

For  I: 


For  R: 


For  M: 


potan0o  = P1tan(0Q  - 0^ 

(2.17) 

PoV^o  = plUlsin((^o  ' 01} 

(2.18) 

Po  + poUo2sin^o  = P1  + P1u12sin2(0Q  - 

9l) 

(2.19) 

ho  + 2 Uo  sin  ^o  hl  2 U1  sin  (0o  - 

9l) 

(2.20) 

P1tan01  = p2tar.(*L  - ©2) 

(2.21) 

plUlSin^l  = p2Uosin^l  ' 02^ 

(2.22) 

P1  + piui2sin2<^i  = P2  + P2U22sin2^l  ’ 

02) 

(2.23) 

hl  + \ U12sin201  = h2  + | U22sin2(01  - 

e2) 

(2.24) 

pQtan03  = p3tan(03  - 03) 

(2.25) 

poU0&in^3  = P3^3s^n( $3  • e3) 

(2.26) 

Po  + p0u02sin2^3  = P3  + p3U32sin2(^3  - 

e3} 

(2.27) 

hQ  + \ UQ2sin203  = h3  + | U32sin2(03  - 

93) 

(2.28) 

boundary  conditions  are: 

P2  = P3 

(2.29) 

e3  = 0X  + e2 

(2.30) 

14. 


It  will  be  shown  in  Section  2-3-2,  using  the  (P,  e)-plane,  that 
two  combinations  (families)  of  SMR  are  possible,  the  usual  one,  in  which  the 
reflected  shock  wave  R redeflects  the  flow  so  that  03  = 0q  - ©2  311  cl.  one  where 
the  flow  is  being  further  deflected  while  passing  through  R,  i.e.,  63  = 0q  + ©2* 
For  Mo  < 1.55  in  nitrogen  (or  M0  < 1.7^  in  argon)  only  the  latter  occurs,  and 
hence  the  plus  sign  (+)  should  be  used  in  Eq.  2.30.  In  the  range  1.55  < Mo  < 
Moc  in  nitrogen  (or  1.71*  < Mo  < Moc  in  argon)  both  combinations  are  possible. 

For  Mo  > Moc  only  the  first  combination  occurs  and  hence  the  minus  sign  (-) 
should  be  used  in  Eq.  2.30.  (Note  Moc  was  introduced  in  Section  2.2.3.1>  where 
the  corresponding  value  for  argon  or  nitrogen  can  be  found.) 

The  above  fourteen  equations  have  eighteen  independent  variables, 
namely  PQ,  Pq,  P2,  P3,  T0,  Tq,  T2,  T3,  UG,  Uq,  U2,  U3,  0o,  0q,  03,  0q,  02  and 
03.  Consequently,  if  four  of  the  eighteen  parameters  are  known  the  remaining 
fourteen  can  be  calculated.  The  four  chosen  parameters  are  again  the  flow 
pressure  Po  and  temperature  T0  ahead  of  the  incident  shock  wave,  the  flow 
velocity  Uo  = Ussec©w  and  the  angle  of  incidence  between  it  and  the  incident 
shock  wave  I,  0o  = 90*-  0w»  Unfortunately,  unlike  the  case  of  RR,  the  initial 
conditions  in  shock  tube  experiments,  namely  Po,  To,  Ms  and  0w  are  not  suffi- 
cient to  define  these  four  chosen  parameters,  since  0w  = 0w  + X involves  one 
more  parameter  q,  which  is  not  known  before  a photograph  of  the  wave  configura- 
tion is  obtained. 

One  way  of  overcoming  the  problem  was  suggested  by  Law  and  Glass 
(Ref.  51)  who  developed  a graphical  method  to  obtain  x = x(po>  T0>  Ms,  By) . 
Their  graphical  method  was  in  good  agreement  with  experiments  (in  oxygen) 
only  in  the  range  25°  < 0w  < ^5°.  However,  the  present  analytical  formulation 
of  their  graphical  method  gives  better  agreement  (in  nitrogen) , over  a greater 
range,  i.e.,  5°  < 0w  < 1+5**  For  the  range  0W  < 5°  a new  method  is  suggested 
(to  be  discussed  subsequently) . 


2. 3. 1.1  Prediction  of  Triple  Point  Trajectory 


As  noted,  a fairly  good  method  for  predicting  the  values  of  x the 
range  25°<  0w  < ^5°  was  developed  by  Law  and  Glass  (Refs.  b8  and  51) , based  on 
the  experimental  observation  that,  except  for  very  strong  shock- wave  diffractions, 
the  Mach  stem  M is  only  slightly  curved.  Consequently,  they  assumed  a straight 
Mach  stem  normal  to  the  wedge  surface,  and  introduced  an  additional  independent 
geometrical  relation  (Fig.  12): 


03  = 90*  - x (2.31) 

This  additional  equation  together  with  the  fourteen  equations  of  motion  (Eqs. 
2.17  to  2.30)  can  now  be  used  to  find  expressions  of  the  form: 


* ■ x(Pe.  V “s’  8„>  <2-32> 

This  is  done  by  keeping  P0,  T0  and  Mo  constants  and  changing  0o  [0O  = 90*  - 
(0W  + x)]  in  small  steps,  and  then  finally  solving  graphically  for  each  set 
of  POJ  T0>  Mq  and  <pQ  the  equations  of  motion  (Eqs.  2.17  to  2.30).  Once  these 
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equations  are  solved  x is  calculated  from  Eq.  2.31  and  Ms  from  the  relation 
Mg  = Mosin0o.  Consequently,  for  a given  value  of  Mo>  each  value  of  9w  + x 
(90°  - 0O)  corresponds  to  a value  of  Ms  and  Repeating  the  above  with 
different  values  of  Mo  but  still  maintaining  the  same  pressure  and  temperature 
(P0  and  To)  will  produce  a set  of  points  (Ms,  x)  f°r  each  ( 0w  + x)*  Thus  x 
vs  Ms  can  be  plotted  for  constant  (&w  + x)-  Using  these  curves  x can 
determined  for  any  Ms  and  0w  at  a given  PQ  and  T0- 

Law  and  Glass  (Ref.  51)  found  that  this  method  is  in  good  agreement 
with  experimental  results  only  in  the  range  25*  < 9W  < 45*;  consequently  they 
suggested  that  an  alternative  method  should  be  developed  for  smaller  wedge 
angles . 

An  alternative  method  for  very  small  wedge  angles  9w  < 4°  is 
suggested  in  the  following.  This  method  is  again  based  on  an  experimental 
fact,  i.e.,  the  fact  that  at  small  wedge  angles  an  SMR  occurs  even  though 
the  flow  behind  the  incident  shock  wave  I is  subsonic  (Mq  < 1)  wrt  the  reflec- 
tion point  P (Fig.  4a).  It  has  been  discussed  previously  that  the  triple-point- 
trajectory  angle  plays  a significant  role  (at  these  small  wedge  angles)  when  it 
makes  the  effective  wedge  angle  0w  = ©w  + X large  enough  so  that  Mq  will  always 
be  greater  than  unity  wrt  the  triple  point  T (Fig.  4b).  Consequently,  one  can 
easily  find  a relation  of  the  form  0w  (Mq=l  = ew  (Ms>  po>  To)  Tor  which  the  flow 
behind  the  incident  shock  wave  is  exactly  sonic,  i.e.,  Mq  = 1.  Once  this 
relation  is  derived,  x can  be  found  from: 

x * 6” Iml»I  e”  (2'33) 

It  will  be  shown  in  Section  7-2  that  for  very  small  wedge  angles  the  value  of  x 
predicted  by  Eq.  2.33  is  in  better  agreement  with  experimental  results  than  the 
one  obtained  by  Law  and  Glass  (Eq.  2.32). 

Note  that  although  Po  and  T0  do  not  appear  explicitly  in  Eq.  2.33 
as  they  do  in  Eq.  2.32,  the  first  term  on  the  RHS  of  Eq.  2.33  depends  on  the 
initial  pressure  and  temperature  and  hence  both  methods  of  predicting  x have 
in  general  the  form  x = x(po,  T0>  Ms,  9w) • 


2. 3-1. 2 Perfect-Gas  Solution 


(2.37) 


1 + 


1 + 


e3  = 02  + 91 


P P P 

_3  =_1  _2 

P P„  PT 
o o 1 


(2.38) 

(2.39) 


This  set  of  six  equations  involves  nine  variables, namely,  7,  Mq,  Mp,  Sq,  02  > e3> 
pl/po,  P2/po  8X1(1  p3/po-  Consequently,  three  variables  must  be  defined  for  a 
solution.  If  7,  Mq  and  Pi/P0  are  selected,  9p  and  Mp  can  be  immediately  cal- 
culated from  Eqs . 2.34  and  2.37,  respectively,  and  the  number  of  equations  to 
be  solved  is  reduced  to  four. 


A further  simplification  was  made  by  Henderson  who  reduced  this 
system  of  equations  of  motion  to  a single  polynomial  of  degree  ten,  with  the 
pressure  ratio  P3/P0  as  the  polynomial  variable.  The  polynomial  coefficients 
were  taken  to  be  functions  of  7,  Mq  and  Pp/PQ.  A detailed  description  is  given 
in  Ref.  31.  Note  that  a polynomial  of  degree  ten  yields  ten  mathematical  roots. 
Henderson  shows,  however,  that  from  simple  physical  consideration  and  the 
possibility  of  double  roots,  seven  of  the  ten  roots  can  be  discarded.  A 
detailed  discussion  of  the  possible  roots  can  be  found  in  Refs.  31  and  48. 


2. 3. 1.3  Imperfect-Gas  Solution 

When  real-gas  effects  are  considered  the  equations  of  motion  (Eqs. 
2.17  to  2.30)  cannot  be  simplified.  Consequently,  the  fourteen  nonlinear 
algebraic  equations  must  be  solved.  A method  of  solution,  very  much  the  same 
to  the  one  utilized  for  RR  (Section  2. 2. 1.2)  was  developed  and  solved  for  the 
first  time. 

Note  that  Eqs.  2.17  to  2.30  consist  of  three  similar  _ sets 
of  four  equations  each  (2.17  to  2.20,  2.21  to  2.24  and  2.25  to  2.28)  which  are 
treated  separately.  The  final  solution  is  obtained  when  Eqs.  2.29  8X1(1  2.30 
are  satisfied.  The  first  set  (Eqs.  2.17  to  2.20)  involves  the  following  eight 
variables:  Po,  T0,  UQ  and  0O  on  the  LHS  and  Pp,  Tp,  U^_  and  0p  on  the  RHS. 

Since  all  variables  on  the  LHS  consist  of  initial  conditions  and  are  known, 
the  RHS  can  be  solved  and  Pp,  Tp,  Up  and  9^  obtained. 

The  second  set  (Eqs.  2.21  to  2.24)  consists  again  of  eight  variables: 
Pp,  Tp,  Up  and  0p  on  the  LHS  and  P2,  T2,  U2  and  02  on  the  RHS;  thus  by  assuming 
a value  of  0p  the  LHS  becomes  known  (P^_,  T^  and  Up  were  calculated  in  the  first 
set)  and  the  RHS  can  be  solved  to  give  P2,  T2,  U2  and  02*  The  third  set  (Eqs. 
2.25  to  2.28)  also  consists  of  eight  variables:  P0,  T0>  UQ  and  0 3 on  the  LHS 
and  P3,  T3,  U3  and  03  on  the  RHS.  It  is  solved  in  the  same  fashion  as  the 
second  set,  i.e.,  assume  a value  of  03  so  that  the  LHS  becomes  known  (PQ,  T0 
and  Uq  are  the  initial  conditions)  and  then  solve  for  P3,  T3,  U3  and  ©3.  The 
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correctly  assumed  values  of  and  <£3  are  those  which  will  result  in  values 
of  Pg,  P3,  01  > 02  63  that  will  satisfy  the  boundary  conditions,  i.e., 

Eqs.  2.29  and  2.30.  Fast  convergence  is  assured  by  using  the  Newton-Raphson 
method.  A congcuter  program  (based  on  the  above  method)  as  well  as  a typical 
output  is  given  in  Appendix  C. 

2.3.2  Graphical  Solution  (Shock- Polar  Presentation)  • 

An  SMP.  in  the  (P,  e)-plane  is  shown  in  Fig.  13 • The  R-polar  does 
not  intersect  the  p/P0-axis,  and  hence  an  RP.  is  not  possible.  State  (2) 
behind  the  reflected  shock  wave  R,  and  state  (3)  behind  the  Mach  stem  M,  in 
the  vicinity  of  the  triple  point  T,  are  at  the  point  where  the  R-polar  inter- 
sects the  I-polar  [state  (2)  on  the  R-polar  and  state  (3)  on  the  I-polar],  since  it 
satisfies  the  equations  of  motion  (Eqs.  2.17  to  2.28)  as  well  as  the  boundary 
conditions  (Eqs.  2.29  and  2-30),  i.e.,  03  = 01  - 02  and  P3  = P2.  Note  that 
since  the  Mach  stem  is  not  straight  it  corresponds  to  a small  portion  of  I 
rather  than  single  point  (3)  , see  Law  and  Glass  (Ref.  51)  for  details. 

It  was  stated  earlier  that  a necessary  condition  for  an  SMR  to 
exist  (once  RR  is  terminated)  is  that  the  flow  behind  the  reflected  wave  R 
(Fig.  12)  is  subsonic  (M2  < 1)  wrt  the  triple  point  T.  This  means  that  the  sub- 
sonic. portion  of  R-polar  (i.e.,  the  one  above  the  sonic  point  Sr)  intersects 
with  the  I-polar.  One  should  note  that  if  the  intersection  was  on  the  supersonic 
portion  of  R (i.e.,  below  Sr)  the  flow  behind  it  wculd  be  supersonic  (M2  > 1) 
and  since  02  / 0 a n additional  deflecting  mechanism  would  be  needed  to  turn 
away  this  supersonic  flow  from  the  wedge  surface,  or  make  it  subsonic  before 
colliding  with  it. 

When  Mo  < MoCj  i.e.,  when  the  R-polar  becomes  tangent  to  the  P/PQ- 
axis  inside  the  I-polar  (Fig.  10a),  the  termination  of  RR  will  always  be 
followed  -by  an  SMR.  In  this  case,  the  R-polar  will  intersect  the  I-polar 
at  a point  above  its  point  of  maximum  deflection,  i.e.,  on  che  strong- shock 
family  portion.  In  this  range  (Mq  < Moc)  the  I and  R-polars  can  take  on  two 
basically  different  combinations.  In  one  of  them  the  solution  in  the  (P,  0)- 
plane  (Fig.  l4a)  indicates  that  the  flow  is  re  deflected  by  R so  that  03  = 02  < ©1> 
while  the  other  combination  (Fig.  l4b)  gives  63  = 02  > 01  and  hence  the  flow 
is  being  further  deflected  by  R.  The  present  modification  of  Eq.  2. 30  (the 
boundary  condition  for  the  equations  of  motion)  that  was  presented  earlier  was 
made  for  this  reason. 

The  limiting  condition  between  these  two  cases  is  the  combination 
where  03  = 82  = 01  (Pig*  l4c) , for  which  there  is  no  deflection  through  R (in 
the  vicinity  of  the  triple  point  T) , i.e.,  R is  normal  to  the  streamlines.  The  analy- 
sis suggests,  furthermore,  that  for  Mo  < 1.55  in  nitrogen  and  Mg  < 1.74  in  argon, 

'only  the  second  combination  of  the  I and  R-polars  (the  one  in  which  03  = 02  > ©l) 
occurs,  while  for  1.55  < Mo  < Mqc  in  nitrogen,  and  1.74  < Mq  < Moc  in  argon, 
both  combinations  are  possible. 

When  Mq  > Moc  (i.e.,  the  R-polar  becomes  tangent  to  the  P-axis  outside 
the  I-polar,  Fig.  10b)  two  interesting  cases  exist.  In  one  of  them  Mqc  < Mo  < 

2 .40,  2 .41  and  2.69  for  perfect  nitrogen,  imperfect  nitrogen  and  argon,  respec- 
tively, the  R-polar  intersects  the  I-polar  in  the  portion  between  its  sonic  point 
and  its  point  of  maximum  deflection,  jdr  (Fig.  15),  which  corresponds  to  the 
stb sonic  portion,  and  hence  11?  < 1 and  an  RMF.  will  still  occur. 
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For  higher  values  of  Mo,  the  termination  of  RE  will  result  in  M2  > 1 
and  hence  an  SMR  is  not  possible.  However,  if  9^,  (or  Q^)  is  decreased  and 
Mq  kept  constant  the  R-polar  moves  away  from  the  P-axis,  and  a situation  is 
reached  in  which  once  again  the  R-polar  intersects  the  I-polar  at  its  sonic 
point  s (Fig.  l6,  Rn)  . Upon  a further  decrease  of  9w,  the  R-polar  intersects 
the  I-polar  baloV  its  sonic  point  (Fig.  16,  Riii)  and  hence  again  an  SMR  is 
obtained.  Thus  one  can  conclude  that  any  value  of  Mo  can  be  matched  with  a 
corresponding  value  of  0O  [0o  =90°  - (9w  + x)]  for  which  Ife  = 1.  The  signi- 
ficance of  this  family  of  SMR- solutions  lies  in  the  experimental  fact  that  it 
corresponds  to  an  SMR  in  which  the  slipstream  is  curled  back,  perhaps  owing 
to  the  extremely  high  flow  deflections  involved.  (For  example,  on  Figs.  10a 
and  10b,  9p  » 1*  and  5°,  respectively,  while  in  Fig.  l6,  9p  > 2?°). 

2.3.3  Termination  of  Single-Mach  Reflection 

As  described  earlier,  the  SMR  can  be  viewed  as  a reflection  process 
through  a three- shock- confluence  (I,  R and  M) , by  which  the  flow  behind  R is 
made  subsonic  (M2  < l) , wrt  the  triple  point  T,  and  hence  can  negotiate  the 
solid  surface.  However,  when  the  Mach  number  of  the  incident  shock  wave  is 
increased,  the  flow  Mach  number  behind  the  reflected  shock  wave  R is  also 
increased,  and  finally  it  exceeds  the  value  of  one.  Consequently,  an  SMR  is 
no  longer  sufficient,  since  again  we  have  a supersonic  flow  (M2  > l)  directed 
towards  a solid  surface  (92  / 0 always).  On  the  (P,  G)-plane  this  means  that 
the  intersection  point  between  the  I and  R-polars  lies  on  the  supersonic  part  of 
the  R-polar  (Fig.  16,  Rp) . Physically,  this  supersonic  flow  will  have  to  be 
made  either  parallel  to  the  solid  surface  or  become  subsonic  before  it  reaches 
the  wall  by  means  of  a shock  wave  (or  a compression  wave)  so  that  it  can  nego- 
tiate the  solid  surface  subsonically.  In  practice  the  second  process  occurs 
and  a CMR  forms,  where  the  flow  Mach  number  behind  the  reflected  shock  wave 
is  reduced  by  passing  through  a compression  wave.  At  higher  values  of  M2 
this  compression  wave  becomes  a ah  ods  wave  and  CMR  terminates.  Henderson  and 
Lozzi  (Ref.  60)  were  the  first  to  suggest  this  sequence  of  events,  i.e.,  a 
band  of  compression  waves  which  can  converge  to  form  a shock  wave.  The 
existence  of  this  compression  wave  was  first  verified  experimentally  during 
this  study.  The  reflection  associated  with  the  conpression  wave  is  a CMR 
while  the  one  associated  with  the  shock  wave  is  a DMR. 

2.4  Complex-Mach  Reflection  (CMR) 

The  CMR  shock-wave  configuration  is  shown  in  Fig.  17.  Unlike  the  case 
of  an  SMR  where  the  reflected  shock  wave  R has  continuous  ( smooth)  curvature 
until  it  finally  terminates  at  the  surface,  here  the  curvature  reverses  dis- 
continuously.  (in  some  CMR  only  a smooth  reversal  of  curvature  takes  place.) 
As  a result  a kink  forms  in  the  reflected  shock  wave.  R is  straight  between 
the  triple  point  T and  the  kink  K,  but  the  remaining  portion  curves  continu- 
ously until  it  terminates  at  the  wall.  It  will  be  shown  that  a compression 
wave  is  located  at  the  kink.  When  it  converges  to  form  a shock  wave  in  DMR, 
the  kink  becomes  the  second  triple  point. 

2.4.1  Graphical  Solution  (shock- Polar  Presentation) 

As  mentioned  previously,  the  existence  of  a band  of  compression 
waves  suggested  by  Henderson  and  Lozzi  (Ref.  60)  was  verified  experimentally 
during  the  present  study.  In  order  to  solve  the  problem  analytically  or 
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graphically  it  is  assumed  that  this  compression  wave  exists  between  the  kink 
K (Fig.  17)  of  the  reflected  wave  R and  the  slij stream  S.  Consequently,  the 
flow  from  region  (2)  is  compressed  (by  passing  through  this  compression  wave) 
to  a new  higher-pressure  state  (region  5).  Furthermore,  we  impose  the  boundary 
condition  that  the  pressure  and  flow  direction  in  states  (4)  and  (5)  are  equal, 
and  hence  they  correspond  to  a single  point  in  the  (P,  e)-plane. 

The  solution  is  no  longer  as  simple  as  in  the  RR  or  SMR  case,  since 
the  kink  and  its  vicinity  should  be  treated  from  a frame  of  reference  attached 
to  the  kink.  Thus  the  solution  must  be  done  in  two  steps:  first,  when  the 
frame  of  reference  is  attached  to  the  triple  point  T,  and  then  when  it  is 
attached  to  the  kink.  Consequently,  the  relative  motion  of  the  kink  wrt  the 
triple  point  should  be  known  (see  Appendix  D) . 

Figure  l8  represents  the  graphical  solution  of  a CMR.  States  (0), 

(l),  (2)  and  (3)  are  found  in  the  same  way  as  for  an  SMR  (see  Section  2.3.2), 
and  states  (l1),  (2')  and  (3')  are  found  in  the  same  way  as  for  a DMR  (see 
subsequent  discussion  in  Section  2.5.2). 

Once  state  (2M  is  found,  the  (p,  e)-curve  of  the  compression  wave 
should  be  drawn  from  (2').  Unfortunately,  real-gas  effects  prevent  derivation 
of  an  explicit  algebraic  formula  of  the  (P,  e)-curve  of  a compression  wave.  In 
order  to  overcome  this  difficulty,  assume  that  the  compression  process  is  perfect 
(reasonable  for  such  a weak  compression)  and  make  use  of  the  Prandtl-Meyer  func- 
tion, -0  + const  = v(M) , where  e is  the  deflection  angle  measured  wrt  the 
direction  of  a given  flow  with  Mach  number  M,  and  v(M)  is  the  Prandtl-Meyer 
function  -(Ref . 27). 

Let  any  point  inside  the  compression  be  denoted  by  i so  that: 


-9^  + const  = v(M^) 


where: 

v(M)  tan  1 Va(r^  - 1)  - tan  "Vm2  - 1 

Equation  2.40  holds  also  in  state  (2'),  Fig.  l8,  and  hence: 

-0^  + const  = v(M£) 

Defining  9g  = 0 results  in: 


(2.40) 


(2.41) 


const  = v(Mp 


(2.42) 


Inserting  Eq.  2.42  into  Eq.  2.40  gives: 


0.  = v(Mp  - v(M.) 


(2.43) 
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The  isentrqpic  relations 


1 


1 7-1 


are  employed  to  obtain: 


i + -2Lf-i  (np2 

1 + ^ Mt"  ;J 


7-1 


(Note  that  P2  = Pg-)  Rearrangement  of  the  above  equations  results  in: 


-1 


M. 


2 + (Hp2 


r 


7 - 1 


7 - 1 


(2.44) 


Inserting  Eq.  2.44  into  Eq.  2.43  yields  a relation  of  the  form: 


where 


9i  = v(N£)  - f(M£,  P2,  Pi)  (2.45) 

/(*£)  =V«  tan_17a[(^)2  - 1]  - tan'^Hp2  - 1 


Since  and  P2  are  known  [state  (2')  was  already  found],  Eq.  2.45  represents 
the  (P,  0) -relation  across  a perfect  compression  wave. 

Once  this  curve  is  drawn  from  state  (2'),  states  (4)  and  (5)  are 
obtained  at  the  point  where  this  curve  intersects  the  R' -polar.  State  (4) 
is  on  the  R' -polar  while  state  (5)  is  on  the  compression  wave  polar. 

Note  that  since  the  compression  wave  under  consideration  is  weak 
it  can  be  replaced  by  a weak  shock  wave.  For  this  case,  the  CMR  graphical 
solution  becomes  the  same  as  the  DMR  graphical  solution  that  is  given  in 
Section  2.5*2. 
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It  has  been  found  in  the  present  study  that  when  the  flow  Mach 
number,  Mg,  behind  the  reflected  shock  wave  R becomes  supersonic  in  a frame 
of  reference  attached  to  the  kink  K (Fig.  17),  then  CMR  terminates  resulting 
in  DMR  (Fig.  19) . The  kink  K becomes  the  second  triple  point  Ti*  The  portion 
of  the  reflected  shock  wave  R,  between  T and  Ti,  is  straight.  The  two  triple 
points  and  their  flow  and  wave  systems  are  shown  in  Fig.  20.  In  the  following 
the  appropriate  expression  for  Mg',  i.e.,  the  flow  Mach  number  behind  the 
reflected  shock  wave  in  a frame  of  reference  attached  to  the  kink  is  developed. 

Consider  Fig.  21  in  which  the  kink  K (or  the  second  triple  point 
Ti)  moves  with  the  velocity  VTiT  (or  Vkt)  wrt  the  first  triple  point  T.  The 
flow  velocity  behind  the  reflected  shock  wave  R,  V2T  s U2  is  directed  along 
a line  that  is  inclined  with  an  angle  <t>i  - 92  to  the  direction  of  motion  of 
the  kink.  Consequently,  using  simple  vector  analysis  and  the  (x,  y) -coordinate 
system  shown  in  Fig.  21, 


1 
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'lZI1  V2TCOS^l  " 62^  ’ VTlT 
■ V2TSln(0l  * 92> 

where  V^-,  and  are  the  x and  y components  of  the  velocity  of  the  flow 

behind  R wrt  the  kink  K or  the  second  triple  point  T^_  • Thus  the  absolute  flow 
velocity  wrt  K becomes: 

lV2T1l  = tV2T  + VT^T  " 2V2TVT1TCOS^l  “ 02^ 

Rearranging  the  above  expression  results  in: 

/ VTqT  ^ 


V = V 
V2T1  V2T 


1 + 


VT-|T  \2 


2T 


- 2 


V J - 92> 


1/2 


and  in  terms  of  Mach  numbers  one  obtains: 


^ = M2 


ni/2 


(2.116) 


The  velocity  of  the  second  triple  point  Tq  wrt  the  first  triple 
point  T,  developed  in  Appendix  D,  is  given  by: 


T-jT  = ^ cosec^i  + ’ ei)U£ 


(2. U7) 


Note  that  the  above  expression  for  Vt^T  (see  Appendix  D for  details)  is  based 
on  an  analytical  method  for  predicting  the  location  of  the  kink  originally 
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suggested  by  Law  and  Glass  (Refs.  48,  51).  Their  method  was  found  to  be  in 
very  good  agreement  with  experiments  only  in  the  range  9w  < 40°,  while  for 
the  range  9w  > 40#  the  agreement  became  progressively  worse  (Ref.  46) . 

2.5  Double-Mach  Reflection  (DMR) 

When  the  flow  Mach  number  Mg  behind  the  reflected  shock  wave  R (in 
a frame  of  reference  attached  to  the  second  triple  point  Ti)  exceeds  the  value 
of  one,  a shock  wave  is  formed  at  the  kink  of  the  CMR,  and  consequently  a DMR 
is  obtained.  The  shock  wave  configuration  of  a DMR  is  shown  in  Fig.  19.  The 
additional  shock  wave  Rp  results  in  a second  triple  point  Ti,  and  a second 
slipstream  S]_,  since  the  flow  behind  Tp  is  obtained  from  state  (l)  by  passing 
through  either  one  shock  wave  (Mp)  or  two  shock  waves  (R  and  Rp) . The  first 
triple  point  trajectory  angle  is  again  x,  while  the  second  triple  point  trajec- 
tory angle  is  x' • 

States  (l) , (2)  and  (5)  can  be  found  by  solving  Eqs.  2.17  to  2.30. 

It  is  worth  mentioning  that  the  major  assumption  behind  Eqs.  2.17  to  2.30  is 
that  states  (0),  (l),  (2)  and  (3)  are  uniform.  Consequently,  one  can  conclude 
that  they  describe  the  first  triple  point  of  a DMR  better  than  that  of  an  SMR 
since  in  a DMR  state  (2)  is  supersonic  and  it  is  obtained  from  state  (l)  through 
a straight  shock  wave  (R) . Recall  that  in  the  case  of  an  SMR,  R is  curved  and 
hence  the  flow  behind  it  is  not  uniform.  Note  that  Eqs.  2.17  to  2.30  are 
applicable  only  in  the  vicinity  of  the  triple  point  with  a confluence  of  curved 
shock  waves,  and  further  away  from  it  when  the  shock  waves  are  straight. 

2.5-1  Analytical  Formulation 


The  analytical  formulation  is  based  on  the  similarity  between  the 
flow  fields  associated  with  the  first  and  second  triple  points.  Consider 
Figs.  20a  and  20b,  which  describe  the  flow  fields  in  the  vicinity  of  T and 
Tp  respectively.  Let  the  notation  a s b mean  that  symbol  'b'  in  Fig.  20b 
is  equivalent  (analogous)  to  symbol  'a'  in  Fig.  20a;  thus  one  can  write: 

state  (l)  s state  (0) , state  (2)  s state  (l) , state  (5)  3 state  (2), 

state  (4)  s state  (3); 


and  finally, 


R=I,  R-i  = R,  M-.  = M,  S,  = S; 


u-  H Uo,  U*  3 U1S  U'  3 U2,  U£  3 U3; 


02  S 9i»  05  H 02’  04  S e3; 


0 P = 0O5  <t>2  = 01  and  0£  = <t>y 


Applying  this  analogy  to  Eqs.  2.17  to  2.30  will  result  in  the  equations  of 
motion  for  the  second  triple  point. 
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For  R: 


p^tan0^  = p2tan(0^  - 9£) 

(2.48) 

P1U1  sin<t>l  = P2U2  sin(^L  " 0P 

(2.49) 

P1  + = P2  + P2U22sin2^i  - 

*2> 

(2.50) 

1^+i  U^2sin20^  = h2  + | U^2sin2(0[  - 

ep 

(2.51) 

For  R^: 

P2tan02  = p^tan(  02  - Qp) 

(2.52) 

P2U^sin0^  = P^U^sin(02  - 0£) 

(2.53) 

P2  + P2Ui2sin2^2  = P5  + P5U52sin2(^  - 

e5} 

(2.54) 

h2  + \ Ug2sin202  = h5  + | U'2sin2(0^  - 

05} 

(2.55) 

For 

P1tan0]^  = p4tan(0^  - e£) 

(2.56) 

P1U|sin01^  = Pl+U^sin(0^  - 04) 

(2.57) 

P1  + Piui2sin2<^  = p4  + P^sin2^  - 

94) 

(2.58) 

\ + \ Ui2sin204  = h4  + | U^sin2^  - 

94) 

(2-59) 

And  finally  the  boundary  conditions  are: 

P5  =P4 

(2.60) 

04  = 0'  + 0' 

(2.61) 

Since  thermodynamic  properties  do  not  depend  on  the  frame  of  reference 
of  the  solution  pi  = Pi»  P2  = P2>  etc,  811(1  hence  the  prime  was  omitted  from  the 
thermodynamic  variables  in  the  above  equations. 
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It  has  been  shown  in  Section  2. 3 .1.3  that  once  POJ  TOJ  U0  and.  0O 
are  known,  Eqs.  2.17  to  2.30  can  be  solved.  Thus,  in  order  to  solve  Eqs. 

2.48  to  2.6l,  Pi,  Tq,  Uq  and  0q  should  be  defined.  In  the  following  expressions 
for  Uq  and  0q  are  developed.  (Note  that  Pq  and  Tq  are  known  from  the  solution 
of  Eqs.  2.17  to  2-30) . 

Consider  Fig.  21  in  which  the  second  triple  point  Tq  moves  along  the 
x-axis  (defined  by  the  direction  of  the  reflected  shock  R)  with  the  velocity 
VTqT  wrt  the  first  triple  point  T.  Simple  vector  analysis  yields  that: 

V1T1  = VlTCOS^l  " VT^T 
V1T1  ' VlTSln,fl 

where  Vq^  and  Vqq-^  are  the  x and  y components  of  the  velocity  of  the  flow  in 
state  (l)  wrt  the  second  triple  point  Tq,  i.e.,  the  x and  y components  of 
Uq.  The  value  of  Uq  then  becomes: 

“i  " <°12  + \ i - 2uivt1toos^)1/2 

whereas  the  value  of  0^  is: 


-1 


0|  = tan  ( - 


U^sin0^ 


I C°S0, 


Inserting  the  value  of  VTqT  (Eq.  D.3  in  Appendix  D)  into  the  above  two  equations 
results  in: 


U1  = (Ui2  + l2us2  ' 21  uqUscos^q)l/2 


(2.62) 


, U sin0 

^i  = tan"  Cu^cTs^-T 


(2.63) 


where  L = Lq/Ls  and  Lq  and  Ls,  defined  in  Fig.  D.la,  should  be  measured  from 
an  actual  interferogram. 


It  is  mentioned  in  Appendix  D that  Law  and  Glass  (Ref.  51)  developed 
an  analytical  method  to  predict  L in  the  range  < 40°  (Eq.  D.12)  . Using  their 
method  one  obtains: 


L 


cosec^  + 0Q 


91> 


(2.64) 
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Thus  for  the  range  0w 40°  Eqs.  2.62  and  2.63  reduce  to; 

2 / Pn  ? 2 

U1  + ) Us  cosec  (^1  + 0o  • 01} 


Ui  = 


- 2U^Ug  f — j cosec(0 ^ + 0o  - 01)cos01 


1/2 


0|  = tan 


U^sin0^ 


U1cos01  - -2  Uscosec(01  + 0Q  - 01) 


A further  singplifi  cation  can  be  made  using  Eq.  2.18  and  the  relation  Us  = Uosin0o 
that  together  result  in 


-2  Us  - UlSin(0o  - 9,) 


Inserting  this  relation  into  the  above  expressions  for  and  0^  results  in 


Ui  =U1 


1 + 


sin  (0Q  - e1) 


- 2 


sin (0O  - 91)cos01  -|  1/2 


sin  (0X  + 0Q  - 01)  sin(01  + 0Q  - e1)  J 


(2.65) 


0^  = tan 


-1 


sin0 


COS01  - sin(0o  - e1)cosec(01  + 0Q  - 01) 
Using  trigonometric  functions  Eq.  2.66  can  be  finally  reduced  to 


(2.66) 


1 


= 01  + 0O  " 91  (2-67) 

Since  all  the  variables  in  Eqs.  2.66  and  2.67  are  known  once  Eqs. 

2.17  to  2.30  are  solved,  the  4 required  parameters,  Pq,  Tq,  Uq  and  0q, for 
solving  Eqs.  2.48  to  2.6l  are  analytically  determined.  It  will  be  shown  later 
that  in  spite  of  the  limitation  on  the  wedge  angle  (9w  < 40°)  that  was  intro- 
duced by  the  use  of  Eq.  2.64  in  the  foregoing  development,  Eqs.  2.66  and  2.67 
can  almost  always  be  used  since  the  DMR  domain  is  centred  in  this  angular 
range. . 

2.5*2  Graphical  Solution  (Shock-Polar  Presentation) 

The  graphical  solution  of  DMR  is  akin  to  CMR  described  in  Section 
2.4.1.  The  solution  is  again  done  in  two  steps,  first  with  the  frame  of 
reference  attached  to  the  triple  point  T and  then  when  it  is  attached  to  the 
second  triple  point  Tq.  The  relative-motion  relations  of  the  second  triple 
point  Tq  wrt  the  first  triple  point  T are  given  in  Appendix  D.  Once  they  are 
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known  the  solution  is  as  follows . First  attach  the  frame  of  reference  to 
the  triple  point  T and  draw  the  I and  R-polars  (Fig. 22a)  in  the  same  way  as 
described  in  Section  2.3-2.  States  (2)  and  (3)  are  at  the  point  where  the  I 
and  R-polars  intersect.  Note  that  the  flow  directions  in  states  (0),  (l), 

(2)  and  (3)  are  all  measured  wrt  the  first  triple  point  trajectory,  i.e., 

0O  = 0,  since  the  flow  in  state  (0)  moves  parallel  to  the  first  triple  point 
trajectory.  Consider  now  Fig.  22b  in  which  e and  0'  are  the  angles  corresponding 
to  the  same  direction  measured  ir.  two  different  frames  of  reference  T and  Tq, 
respectively.  Simple  geometrical  consideration  yields: 


0 = 0'  + (x'  ~ X)  (2.68) 


or  alternatively: 

0 = 6'  + (0O  - <t>'Q)  (2-69) 

where  <po  = 90°  - 9w  - X and  0o  = 90*  - 9w  - x'  (0o  is  given  by  Eq.  D.l4)  . 
Rearranging  Eq.  2.69  results  in: 


0'  = e - (0O  - 0;)  (2.70) 

Consequently,  state  (l)  in  a frame  of  reference  attached  to  Ti  is  at  0^  = 0q  - 
(0o  “ 0o)  and  pq  = P1  (n0  change  in  thermodynamic  properties  by  changing  a 
frame  of  reference).  Once  state  (l')  is  known,  the  R'-polar  (i.e.,  the  R-polar 
in  a frame  of  reference  attached  to  Tq)  that  corresponds  to  Mq  (Eq.  2.65), 

Pq  = Pq  and  T-j_  = Tq  is  constructed.  States  (2')  and  (3')  are  at  the  point 
where  a constant  pressure  line  drawn  from  states  (2)  and  (3)  intersects  the 
R'-polar.  Once  state  (2')  is  known  the  Ri-polar  (l^g*  , Pg>  = P2  and  Tg>  ~ Tg) 
is  drawn.  States  (4)  and  (5)  (Fig.  19)  are  at  the  point  where  the  Ri-polar 
intersects  the  R'-polar.  The  initial  conditions  for  the  case  illustrated  in 
Fig .22a  are  Ms  = 3-681,  ew  =.40°,  x = 4.8°,  PQ  = 15-3  torr  and  Tc  = 297.4  K. 
Consequently  the  corresponding  values  for  the  I polar  are  Mo  = 5-l89>  po  = p5 -3 
torr  and  T0  = 297-4  K.  The  first  triple-point  solution  results  in  Mq  = 1.984, 
Pi  = 239.6  torr,  Ti  = 1062.2  K,  P2  = 474.2  torr,  T2  = 1302-9  K,  01  = 32.26* 
and  0q  = 43.10°.  Therefore  the  R-polar  that  corresponds  to  Mq,  Pq  and  Tq  is 
drawn  from  a point  0 = 0q  on  the  I polar.  Using  the  above  listed  values  and 
Eq.  D.l4  results  in  <p'Q  = 42.6°  and  hence  (l')  is  shifted  by  2.6°  (0q  - 0q) 
from  state  (l)  along  a constant  pressure  line.  With  the  help  of  Eq.  2.65  and 
the  above  values  one  finds  that  Mq  = 0.8237  x Mq  = 1.634  and  the  R'-polar 
that  corresponds  to  Mq,  Pq  and  Tq  is  drawn  from  state  (l').  States  (2')  and 
(3')  are  at  the  point  where  a constant  pressure  line  drawn  from  (2)  and  (3) 
intersects  the  R'-polar,  i.e.,  Po>  = Pg  = 4.74.2  torr.  The  present  solution 
from  the  R'-polar  yields  Pg  = 474.2  torr,  T2  = 1269.6°,  m£  = 1.125.  Conse- 
quently a polar  corresponding  to  these  values  was  drawn  from  states  (2')  and 
(3')  to  obtain  states  (4')  and  (5')  or  (4)  and  (5)  at  the  intersection  of  this 
polar  with  the  R'-polar.  Note  that  although  T2’  f Tg  probably  due  to  the  use 
of  Eq.  2.64,  the  relative  error  is  only  2.55$. 
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2.5-3  Prediction  of  the  Second  Triple-Point-Trajectory  Angle  x' 

The  second  triple-point- trajectory  angle  x'  can  he  expressed  by 
using  Fig.  D.la  as  follows: 


¥'  = 90°  - e - 0' 
A w *o 


Inserting  the  value  of  (f>'0  (Eq.  D.l4)  into  Eq..  2.71  results  in: 


(2-71) 


1 - ^2 


X’  = 90°  - 


0 - tan 

w 


-1 


pl 


c°t0o  " c°t  (0j_  + 0O  - ex) 


(2.72) 


Thus  for  a set  of  given  initial  conditions  P0j  To>  Ms  and.  9w  the  value  of  x 
(the  first  triple-point-tra j ectory  angle)  can  be  calculated  using  the  appropriate 
method  (given  in  Section  2.3 .1.1).  Once  x is  known,  $q,  0q  and  pq  can  be  found 
by  solving  Eqs.  2.17  to  2.30.  (Note  that  0O  = 90°  - 0w  - X 8114  Mo  = Msc°sec0o) , 
and  consequently  x’  can  be  calculated  from  Eq.  2.72. 


It  should  be  mentioned  again  that  although  this  method  of  predicting 
X’  is  limited  to  0W  < 40°  (see  Appendix  D for  details)  it  can  be  used  here 
since  the  DMR-domain  exists  mainly  in  this  range  (see  subsequent  discussions). 
Recall  that  ^pendix  Dwas  developed  in  general  for  the  kink  K of  a CMR  as  well 
as  the  second  triple  point  Tq  of  a DMR.  Consequently  this  method  of  predicting 
X1  is  applicable  also  for  a CMR  for  which  x’  is  defined  as  the  kink  trajectory 
angle. 


2.5.4  Termination  of  Double-Mach  Reflection 


Following  the  gasdynamic  reasons  for  the  formation  of  the  kink  in  a 
CMR  (i.e.,  termination  of  SMR)  where  the  flow  in  state  (2)  behind  R (Fig.  12) 
became  supersonic  wrt  T,  and  then  the  formation  of  DMR  (i.e.,  termination  of 
CMR)  when  the  flow  in  state  (2)  became  supersonic  wrt  K (Fig.  17),  one  might 
assume  that  if  the  flow  Mach  number  in  state  (4)  behind  Mq  exceeds  the  value  of 
unity  wrt  the  second  triple  point  Tq  (Fig.  19)  a new  kink  will  form  on  Mq  and 
DMR  will  terminate.  If  the  flow  in  state  (4)  becomes  supersonic  wrt  this  new 
kink  a triple-Mach  reflection  (TMR)  can  form. 

In  order  to  verify  the  above  hypothesis,  reflections  must  be  obtained 
using  very  strong  incident  shock  waves  as  well  as  long  compression  models  in 
order  to  allow  the  shock  wave  configuration  to  develop  to  a significant  size. 

Unfortunately,  experiments  with  high  Mach  numbers  (Mg  > 8)  involve 
the  risk  of  damaging  the  interferometric-quality  windows.  Consequently  it  was 
not  possible  to  verify  the  existence  of  TMR  in  the  present  facility. 

2.6  Nonstationary-Oblique  Shock-Wave  Reflection  in  (Mg,  0^)-Plane 

The  formation  and  termination  criteria  for  RR,  SMR,  CMR  and  DMR 
having  been  established  can  be  shown  on  the  (MS)  0^) -plane  to  give  the  domains 
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of  the  different  reflection  processes  and  their  boundaries.  A necessary 
condition  for  a reflection  to  occur  is  that  the  flow  behind  the  incident  shock 
wave  I must  be  supersonic  (Mi  > l)  wrt  the  reflection  point  P (in  the  case  of 
RR)  or  the  triple  point  T (for  the  other  reflections)  (Figs.  4a  and  4b,  respec- 
tively) . Thus  the  (Ms,  0y)-plane  [(Mo,  0o)-plane]  can  be  divided  into  two  major 
domains ■ One  of  them  corresponds  to  Mi  < 1,  where  shock-wave  reflection  is  not 
possible  and  the  other  corresponds  to  Mi  > 1 where  shock  reflections  occur. 

This  region  of  shock-wave  reflection  is  then  divided  into  four  domains  corre- 
sponding to  the  four  reflection  processes,  i.e.,  RR,  SMR,  CMR  and  DMR. 

The  criterion  for  predicting  the  termination  of  RR  makes  use  of  the 
boundary  condition  that  the  flow  downstream  of  the  reflection  point  (Fig.  7) 
must  be  parallel  to  the  wall,  i.e.,  02  + 01  = 0 (Eq.  2.9).  When  this  condition 
is  violated  (i.e.,  0 0 increases  to  a point  where  0i  exceeds  the  detachment  angle 
(92m)  in  state  (l)  (Fig.  7),  RR  terminates.  Therefore,  the  termination  criterion 
is: 

el  + 02m  = 0 ^-*73) 

When  RR  terminates,  3 different  types  of  reflection,  i.e.,  SMR,  CMR 
and  DMR,  can  occur  depending  on  the  Mach  number  of  the  flow  behind  the  reflected 
shock  wave  R.  As  long  as  the  flow  behind  R is  subsonic  wrt  the  first  triple 
point  T,  an  SMR  occurs.  When  this  flow  becomes  supersonic  wrt  T,  SMR  terminates 
and  a CMR  forms.  Consequently,  the  termination  criterion  for  SMR  and  the  forma- 
tion criterion  for  CMR  is: 

Mg  = 1 (2.74) 

CMR  terminates  when  the  flow  behind  R becomes  supersonic  wrt  the  kink  K.  There- 
fore the  termination  criterion  for  CMR  and  the  formation  criterion  for  DMR  is: 

m;  = 1 (2.75) 

It  is  worthwhile  mentioning  that  the  line  = 1 corresponds  approximately  to 
M2  = 1.30  in  both  nitrogen  and  argon.  Alternatively  one  may  use  the  following 
empirical  criteria  for  the  existence  of  SMR,  CMR  and  DMR  in  nitrogen  and  argon. 

SMR  occurs  only  if: 

Mg  < 1 (2.76) 

A CMR  takes  place  when: 

1<H><1.3  (2.77) 

A DMR  results  for  all 

Hp  > 1-3  (2.78) 

In  the  following  the  domains  of  RR,  SMR,  CMR  and  DMR  in  the  (Ms,  0^,)- 
plane  are  given.  Note  that  0^  = By  + \ in  the  domains  of  SMR,  CMR  and  DMR 
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while  in  the  RR  domain  9^,  = ew  since  7 = 0 by  definition.  Note  also  that 
since  Mo  = Mg  sec0^  and  0O  = 9°*  - 6^  one  might  assume  that  the  reflection 
domains  amply  to  steady  flows  as  well.  However,  it  will  be  shown  later  that 
in  a steady  flow  CMR  and  DMR  cannot  occur  due  to  geometrical  limitation,  and 
that  the  termination  criterion  of  RR  is  different  to  that  described  by  Eq. 

2.73.  Therefore  the  following  domains  are  for  : or s t at i on ar y flows  only. 

2.6.1  Regions  and  Transition  Boundaries  in  a Diatomic  Gas-Nitrogen 

The  domains  of  RR,  SMR,  CMR,  BMP.  and  NR  and  the  boundaries  between 
them  in  the  (Mg,  9^) -plane  for  nitrogen  are  shown  in  Fig.  23.  The  dashed 
boundary  lines  are  for  a perfect  gas  (7  = 1.4)  while  the  solid  lines  account 
for  real-gas  effects  (vibration,  rotation- vibration  coupling  and  dissociation) 
with  four  different  initial  pressures:  P0  = 1,  10,  100  and  1000  torr  and  a 
constant  initial  temperature  T0  = 300.  It  can  be  seen  that  the  imperfect-gas 
boundary  lines  start  to  diverge  from  those  of  a perfect  gas  at  very  low 
incident  shock  wave  Mach  numbers  owing  to  temperature- dependent  vibrational 
excitation.  At  higher  values  of  Mg  each  boundary  line  splits  and  diverges 
according  to  its  initial  pressure  as  a result  of  dissociation.  At  still  higher 
shock  Mach  numbers  electronic  excitation  and  ionisation  would  play  a similar 
role.  Under  such  conditions  the  sharp  boundary  lines  which  exist  for  a perfect 
gas  between  the  domains  of  the  different  types  of  reflection  are  replaced  by 
a multiplicity  of  lines  depending  on  the  initial  pressure.  For  example , one 
should  expect  an  RR  for  Ms  = 10  and  0^  =45°  when  the  initial  pressure  is  1000 
torr  or  more  and  a DMR  when  P0  = 100  torr  or  less.  As  Mg  approaches  unity  the 
line  0p  + §2  approaches  the  line  Mi  = 1.00.  These  two  lines  are  coincident 
at  the  origin  Ms  = 1,  0^  =0. 

It  can  be  seen  from  Fig.  23  that  if  the  shock  Mach  number  Mg  is 
fixed  while  varying  the  effective  wedge  angle  9^  (i.e.,  the  actual  wedge  angle 
0W) , then  different  domains  of  reflection  are  encountered.  This  figure  clears 
up  same  problems  and  disagreements  between  various  investigators  who  reported 
different  sequences  of  events  as  one  passed  through  a range  of  wedge  angles 
9W  for  a fixed  Mg.  They  were  unaware  of  the  domains  and  boundaries  presented 
in  Fig.  23.  The  different  sequences  of  events  (and  experimental  and  numerical 
investigations)  are  summarised  as  follows:  1.00  < Ms  < 1.60  (see  point  a) 

RR  ->  SMR  (Bleakney  and  Taub,  Ref.  10),  1.60  < Ms  < 2.69  (see  point  b) , RR 
CMR  -»SMR  (Smith,  Ref.  6),  for  Ms  > 2.69,  RR  -+DMR  -♦CMR  -♦SMR  (Kutler  et  al. 
Ref.  63,  and  Shankar  et  al,  Ref.  64) . RR  and  SMR  can  occur  with  a proper 
choice  of  Qy  (or  0W)  for  any  incident  shock  wave  Mach  number  Mg,  while  the 
other  reflections  are  limited  to  defined  values  of  Ms*  e.g.,  CMR  occurs  only 
for  Mg  > 1.60  (see  point  a)  and  DMR  only  for  Mg  > 2.69  (see  point  b) . 

2.6.2  Regions  of  Transition  Boundaries  in  a Monatomic  Gas-Argon 

The  domains  of  RR,  SMR,  CMR,  BMR  and  NR,  and  the  boundaries  between 
them  in  the  (Ms , 0w) -plane  for  argon  are  shown  in  Fig . 24 . The  dashed  boundary 
lines  are  for  a perfect  gas  (7  = 5/3)  while  the  solid  lines  account  for  real-gas 
effects  (electronic  excitation  and  ionization)  with  four  different  initial 
pressures:  PQ  = 1,  10,  100  and  1000  torr  and  a constant  initial  temperature 
To  = 300  K. 

Unlike  a diatomic  gas  (Fig.  23)  the  imperfect-gas  boundary  lines 
start  to  diverge  from  those  of  a perfect  gas  only  at  relatively  high  incident 
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shock  wave  Mach  numbers.  Therefore,  the  perfect-gas  theory  is  adequate  in  the 
range  Mg  < 5-36  for  P0  > 1 torr,  T0  = 300  K (the  point  where  the  RE  imperfect- 
gas  boundary  lines  start  to  diverge) . It  should  be  noted,  however,  that  as  the 
initial  pressure  (PQ)  decreases  the  corresponding  imperfect-gas  boundary  line 
splits  at  a lower  value  of  Mg.  For  example,  the  imperfect- gas  boundary  lines 
for  the  termination  of  SMR  (tfe  = 1)  split  at  Ms  = 7-62,  8.15,  8.50  and  8.77 
for  P0  = 1,  10,  100  and  1000  torr,  respectively.  Therefore,  the  sharp  boundary 
lines  which  exist  for  a perfect  gas  between  domains  of  different  reflections 
are  replaced  by  a multiplicity  of  lines  depending  on  the  initial  pressure  (when 
the  temperature  is  fixed) . Consequently,  one  can  expect  an  SMR  for  Ms  = 10  and 
9w  = 30*  when  PQ  > 100  torr  and  a CMP.  when  P0  < 100  torr. 

It  can  again  be  seen  from  Fig.  2k  that  if  the  incident  shock  wave 
Mach  number  is  fixed  while  varying  the  effective  wedge  angle  0^  (or  ew) , different 
domains  of  reflections  are  encountered.  The  different  sequences  of  events  are: 
1.00  < Mg  < 1.85,  RR  -» SMR  (point  a),  1.85  < Ms  < 3-17,  RR  -+CMR  -* SMR  (point  b) , 
for  Mg  > 3-17,  RR  -+DMR  -»CMR  -+SMR.  RR  and  SMR  can  occur  with  a proper  choice 
of  0W  (or  ew)  for  any  incident  shock  wave  Mach  number  Ms,  while  the  other 
reflections  are  limited  to  definite  values  of  Ms,  e.g.,  CMR  can  occur  only  for 
Ms  > 1.85  (point  a)  and  DMR  only  for  Ms  > 3*17  (point  b) . 


3-  FLOW  DEFLECTION  PROCESS  OVER  A WEDGE  CORNER 

It  has  been  mentioned  earlier  (Section  2.1)  that  the  diffraction 
process  in  nonstationary  flows  also  depends  on  the  flow-deflection  process  of 
the  flow  induced  by  the  incident  shock  wave  over  the  corner. 

Consider  a planar  shock  wave  propagating  in  a shock  tube  and  denote 
the  state  behind  it  as  (2').  For  any  given  set  of  initial  conditions  (PQ,  TQ) 
and  incident  shock  wave  Mach  nuntoer  M. , the  induced  flow  Mach  number  • as 
well  as  the  pressure  and  temperature  (Pp1  and  T2O  can  be  calculated.  Conse- 
quently, the  corresponding  sonic  deflection  angle  eS2<  and  the  angle  of  maximum 
deflection  (detachment  angle)  9^,  can  be  determined. 

Thus  the  (Ms,  ew) -plane  is  now  divided  into  two  main  regions;  one 
corresponds  to  Mg'  < 1,  where  the  induced  flow  is  subsonic  and  hence  turns 
over  the  comer  subsonically  (Fig.  5a),  and  the  other  corresponding  to  Mg'  >1 
where  the  flow  is  supersonic.  The  latter  region  is  subdivided  into  three  regions 
of  different  flow  deflection  processes  (Ref.  36); 

(1)  0 < 9u  < eS2,  for  deflection  through  a straight  and  attached  oblique  shock 
wave  (Fig.  5*0 . 

(2)  6s2'  < 9w  < 9mp,  for  deflection  through  a curved  and  attached  shock  wave 
(Fig.  5c). 

(3)  9w  > 9mp' , where  the  deflection  is  through  a curved  and  detached  shock 
wave  (Fig.  5d) . 

Since  the  maximum  separation  between  ©mp'  and  eS2<  is  usually  very  small  (for 
nitrogen  at  PQ  = 15  torr  and  T0  = 300  K,  - 0S2'  = 0.63*)  only  two  regions 

0 < ew  < 00^,1,  where  the  shock  wave  is  attached,  and  > 9^ ,,  where  the  shock 
wave  is  detached,  need  be  considered  for  practical  purposes. 
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3 -1  Flow  Deflection  Process  in  a Diatomic  Gas  - Nitrogen 

The  above-mentioned  regions  of  different  types  of  deflection  processes 
for  nitrogen  are  shown  in  the  (Ms,  ew)-plane  in  Fig.  25.  The  dashed  line  is 
for  a perfect  gas  ( y = 1.4),  while  the  solid  lines  are  for  imperfect  nitrogen 
with  different  initial  pressures  (F0  = 1,  10,  100  and  1000  torr)  and  fixed 
initial  temperature  (T0  = 300  K) . One  should  note  that  the  line  Mg'  =1  corre- 
sponds to  Ms  = 2.068  (for  a perfect  gas)  and  2.055  when  vibrational  energy  is 
included.  Since  the  difference  between  these  two  values  of  Ms  is  too  small  to 
be  plotted,  only  one  line  is  shown.  It  can  be  3aer.  from  Fig.  25  that  the 
deflection  process  depends  strongly  on  real-gas  effects  (i.e.,  initial  pressure 
Po  for  a given  temperature  T0) • For  Ms  = 10  and  9W  = 35®  the  flow  will  nego- 
tiate the  corner  through  a detached  shock  wave  if  PQ  < 1 torr  or  through  an 
attached  shock  wave  if  PQ  > 10  torr. 

Since  the  imperfect-gas  boundary  line  diverges  from  the  perfect 
gas  even  at  the  lowest  value  of  Ms  one  must  conclude  that  the  perfect-gas 
theory  is  not  adequate  and  imperfect-gas  theory  should  be  used. 

3 -2  Flow  Deflection  Process  in  a Monatomic  Gas  - Argon 

The  regions  of  different  types  of  reflection  processes  for  argon 
in  the  (Ms,  0W) -plane  are  shown  in  Fig.  26.  .Again  the  dashed  line  is  for  a 
perfect  gas  (7  = 5/3)  while  the  solid  lines  are  for  imperfect  argon  with 
different  initial  pressures  (P0  = 1,  10,  100  and  1000  torr)  and  a constant 
initial  temperature  (T0  = 300  K) . 

Unlike  the  case  of  a diatomic  gas  (see  Fig.  25)  the  imperfect-gas 
boundary  lines  start  to  diverge  from  the  perfect-gas  lines  at  quite  high  values 
of  Mg  (i.e.,  Ms  = 7.36)  and  consequently  the  perfect-gas  theory  is  adequate 
in  the  range  Ms  < 7-36.  It  is  worthwhile  mentioning  that  the  line  1^2'  =1 
corresponds  to  Ms  = 2.758  for  both  perfect  and  imperfect  gases. 

The  dependence  of  the  deflection  process  on  real-gas  effects  (at 
high  values  of  Ms)  is  again  clearly  seen  (Fig.  26).  For  Ms  = 10  and  9w  = 10® 
the  flow  will  negotiate  the  corner  through  an  attached  shock  wave  if  Po  < 10 
torr  or  through  a detached  shock  wave  if  PQ  > 100  torr.  Note  that  unlike 
a diatomic  gas  (Fig.  25)  the  region  that  corresponds  to  a deflection  with  the 
aid  of  an  attached  shock  wave  is  very  small  compared  to  the  regions  of  detached 
shock  wave  and  subsonic  turning.  Therefore,  in  general  (unless  the  wedge 
angles  are  very  small  and  the  initial  pressures  very  low) , only  two  flow 
deflection  processes  will  be  observed,  i.e.,  a subsonic  turning  over  the  corner 
or  a turn  with  the  aid  of  a detached  shock  wave. 


4.  SHOCK-WAVE  DIFFRACTION  PROCESS  IN  NONSTAIICHARY  FLOWS 

The  two  independent  phenomena  discussed  in  the  previous  chapters  , 
i.e.,  shock  reflection  (Section  2)  and  fxow  deflection  (Section  3)  interact 
and  give  rise  to  the  overall  shock-wave  diffraction  phenomenon.  To  show  this 
process  the  figures  corresponding  to  shock-wave  reflection  (Figs.  23  and  24 
for  nitrogen  and  argon,  respectively)  and  the  figures  for  flew  deflection 
(Fig.  25  for  nitrogen  and  Fig.  26  for  argon)  are  superimposed.  However,  the 
vertical  axis  of  Figs.  23  and  24  and  that  of  Figs.  25  and  26  are  different. 
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In  Figs.  23  and  24  the  vertical  axis  is  9^  while  in  Figs.  25  and  26  it  is 
9W.  The  effective  wedge  angle  9w  equals  (the  actual  wedge  angle)  in  the 
RR-regime  and  9w  + X elsewhere.  Thus,  the  triple  point  trajectory  angle  % 
should  be  subtracted  from  the  appropriate  curves  of  Figs . 23  and  24  prior  to 
any  superposition. 

The  interaction  between  the  shock-wave  reflection  phenomenon  and  the 
induced-flow  deflection  process  causes  the  reflected  shock  wave  R to  curl  back 
towards  the  model  and  terminates  at  the  wedge  corner  on  the  shock- tube  wall. 

Since  the  shock-wave  configuration  is  growing  with  time,  the  point  where  R 
terminates  at  the  shock- tube  wall  moves  towards  the  oncoming  shock-induced 
flow,  therefore  increasing  the  oncoming  relative  flow  Mach  number.  Consequently, 
the  subsonic  turning  region  shown  in  Figs.  25  and  26  cannot  be  established  in 
nonstationary  flows . At  the  limiting  case  of  a degenerated  incident  shock  wave 
(Ms  = 1)  the  reflected  shock  wave  becomes  a Mach  wave.  Therefore,  as  there  are 
four  reflection  processes  (RR,  SMR,  CMR  and  DMR)  and  two  deflection  processes 
(an  attached  or  detached  shock  wave),  a maximum  of  eight  different  shock- 
diffraction  systems  are  possible. 

In  Section  2.3 .1.1  two  methods  of  predicting  x were  suggested.  The 
first  one  is  due  to  Law  and  Glass  (Ref.  51),  where  straight  Mach  stem  (normal 
to  the  wedge  surface)  is  assumed,  and  the  second,  developed  here,  makes  use 
of  the  experimental  fact  that  shock-wave  reflections  occur  even  at  very  low 
wedge  angles  when  the  flow  Mach  number  behind  the  incident  shock  wave  I (Fig. 

4b)  is  much  smaller  than  unity  (Mq  < 1)  wrt  the  reflection  point.  The  methods 
were  discussed  in  detail  in  Section  2. 3. 1.1.  Consequently  only  the  results, 
i.e.,  x = x(0w>  Ms,  P0,  To)  for  nitrogen  and  argon  will  be  presented  in  the 
following. 

4 .1  Triple-Point-Trajectory  Angle  in  a Diatomic  Gas  - Nitrogen 

The  present  analytical  extension  of  the  graphical  method  by  Law  and 
Glass  (Ref.  51)  for  predicting  the  triple-point-trajectory  angle  x was  applied 
to  obtain  Figs.  27  to  31- 

The  triple  point  trajectory  angles  x as  a function  of  the  incident 
shock  wave  Mach  number  Ms  with  the  effective  wedge  angle  9^  as  a parameter  for 
a constant  P0  and  To  shown  in  Fig.  27.  The  dashed  lines  are  for  a perfect 
gas  (7  = 1.4)  while  the  solid  lines  are  for  an  imperfect  gas.  Since  the 
perfect-gas  model  results  in  no  reflection  for  q(!,  = 22.2*  (see  Fig.  23.)  only 
a solid  line  (i.e.,  inperfect  gas)  is  shown  for  this  effective  wedge  angle. 

It  can  be  seen  that  the  perfect  gas  result  yields  values  of  x greater  than 
for  an  inperfect  gas.  At  low  Mach  numbers  the  inperfect-gas  lines  approach 
those  of  a perfect  gas.  At  higher  Mach  numbers  and  lower  9^  the  dependence  of 
X upon  Ms  becomes  stronger  for  an  inperfect  gas  due  to  dissociation.  However, 
for  a perfect  gas  x is  almost  independent  of  Ms*  (Note  that  at  higher  Mg  the 
dashed  lines  level  out.) 

At  9w  = *+0°,  the  perfect  gas  gives  an  almost  constant  value  of  x through- 
out the  entire  region  of  Ms.  In  other  words,  one  might  conclude  that  a wedge 
angle  exists  for  which  x is  independent  of  Mg  (a  constant).  This  fact  is  more 
clearly  illustrated  in  Fig.  28  showing  the  variation  of  x with  Qy  for  a given  Ms. 
Only  3 lines  (Ms  =4,7  and  10)  are  drawn  for  clarity.  Note  that  all  3 lines 
intersect  practically  at  one  point  (9w  » 35.05°,  x 6.45  and  9w  = 41.50°). 

A reproduction  of  Fig.  28  for  inperfect  nitrogen  is  shown  in  Fig.  29;  again  only 
3 curves  are  given  (Ms  =4,7  and  10) . The  lines  corresponding  to  Ms  = 4 and  7 
intersect  at  9w  = 45-38°  and  x = 3.30°  (i.e.,  9 w = 48.68°).  A comparison  between 
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inperfect  and  perfect  gas  results  is  shown  in  Fig.  30.  Only  the  lines  that 
correspond  to  Ms  = 10  are  shown.  It  can  be  clearly  seen  that  the  value  of  x 
predicted  using  a perfect  gas  is  greater  than  for  an  imperfect  gas. 


The  dependence  of  x on  Ms  and  By  is  shown  in  Fig.  31 • The  dashed 
lines  are  again  for  a perfect  gas  {y  = 1.4)  while  the  solid  lines  are  for 
imperfect  nitrogen.  It  can  again  be  seen  that  the  perfect  gas  values  of  x 
are  higher  than  those  predicted  by  the  inperfect  gas  theory.  The  discrepancy 
between  the  two  values  (i  .e . , perfect  and  inperfect  nitrogen)  increases  as 
Mg  increases  and  ew  decreases.  Again  the  perfect  gas  lines  level  out  at 
high  Mach  numbers  and  x is  almost  independent  of  Mg  or  x ~ x(9w)  only- 

Law  and  Glass  developed  this  method  of  predicting  v and  found  that 
it  agrees  with  experiments  only  in  the  range  25® 9w  ^ 45®  (Ref*  51)*  It 
will  be  shown  later  that  during  the  present  research  their  method  was  found 
to  be  in  good  agreement  with  experiments  in  a much  larger  range  of  wedge 
angles  5®  < 0W  < 45°.  This,  we  believe,  is  due  to  the  fact  that  they  solved 
Eqs.  2.17  to  2.31  graphically,  while  during  the  present  study  an  analytical 
solution  was  obtained.  Their  graphical  solution  was  rechecked  and  verified. 
Since  their  method  was  found  to  be  good  in  the  range  5*  9w  ^ 45*  an  alter- 
native method  was  developed  for  the  range  < 5*  which  was  described  in  detail 
in  Section  2. 3 .1.1.  The  results  are  as  follows. 

Figure  32  represents  the  variation  of  the  triple  point  trajectory 
angle  ^ with  the  incident  shock  Mach  number  Ms  and  the  actual  wedge  angle 
Ow,  as  a parameter,  for  very  small  wedge  angles  ( 0w  < 5*)-  Note  that 
although  Fig.  32  is  based  only  on  the  experimental  fact  that  SMR  occurs  even 
for  very  low  wedge  angles  (i.e.,  x should  be  large  enough  to  keep  9^  above 
the  line  Mi  = 1 of  Fig.  23;  recall  that  % = 9^  - 9^),  it  resembles  the 
characteristics  of  Fig.  31,  for  small  wedge  angles  5°<  0w  < 20°  that  results 
from  physical  considerations.  In  both  figures  ^ is  a decreasing  function  with 
increasing  Ms  and  9w  and  while,  for  the  inperfect  gas  model,  the  dependence 
of  x upon  Ms  is  prominent  as  Ms  increases  due  to  dissociation,  the  perfect 
gas  lines  (dashed)  level  out  so  that  x becomes  nearly  independent  of  Ms,  i.e., 

X = x(®w)  only  (Fies • 31  and  32) . 

4 .2  Triple-Point-Trajectory  Angle  in  a Monatomic  Gas  - Argon 

The  dependence  of  x on  Ms  for  fixed  9^  is  shown  in  Fig.  33*  The 
solid  lines  for  imperfect  argon  in  ionisational  equilibrium  start  to  diverge 
from  the  perfect-gas  lines  y = 5/3  (dashed)  slightly  above  Ms  = 7.  Conse- 
quently, in  the  range  1 < Ms  < 7 the  perfect-gas  solution  is  adequate. 

The  perfect-gas  lines  level  out  as  Mg  increases  making  v independent 
of  Ms  on  x = x(0w)  only*  This  can  be  more  clearly  seen  in  Fig.  34,  in  which 
X is  plotted  against  9^-  for  fixed  values  of  Ms.  For  clarity  only  3 lines 
axe  drawn  (Ms  =4,5  and  10) . All  the  lines  intersect  practically  at  one 
point,  9w  = 29.87®,  x = 10.22®  (i.e.,  9w  = 40.09®) . This  value  differs  very 
slightly  from  the  corresponding  value  of  a diatomic  gas  found  earlier  (9w  = 
41.50®) . 


A similar  plot  for  imperfect  argon  is  given  in  Fig.  35*  Again  only 
3 lines  (Mg  =4,7  and  10)  are  shown.  Only  two  lines,  those  corresponding 
to  Ms  = 4 and  7 intersect  (at  9w  = 32.72®,  x ~ 8-83®,  9^  = 41.55®) • Therefore 
a wedge  angle  9w  does  not  exist  for  which  x Is  completely  independent  of  Mg. 

A comparison  between  perfect  and  inperfect  argon  is  shown  in  Fig.  36.  It  can 
be  clearly  seen  that  the  value  of  x predicted  by  perfect-gas  theory  is  larger 
than  the  one  for  inperfect  argon . 
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The  depender.ee  of  ^ on  M£  for  giver.  ^ is  shown  in  Fig.  37-  Again 
the  inperfect  argon  (solid  line)  starts  to  diverge  from  the  perfect  gas  values 
7 = 5/3  (dashed  lines)  only  for  Ms  > 7-  Therefore  the  perfect  gas  values  are 
adequate  in  the  range  1 < Ms  < 7-  The  perfect-gas  lines  level  out  as  Ms 
increases,  giving  j_  = on4y . However,  electronic  excitation  and  ior.ir,a- 

tion  start  when  Ms  > 7 and  the  depender.ee  of  upon  Ms  increases  for  imperfect 
argon. 


The  previous  method  for  very  low  wedge  angles  was  also  used  for  argon 
to  obtain  Fig.  38-  The  solid  lines  are  for  imperfect  argon  in  ioni nation 
equilibrium  while  the  dashed  lines  are  for  a perfect  monatomic  gas.  It  can 
again  be  seen  that  for  Mg  > 2,  ^ is  a decreasing  function  of  Mg  and  9w>  while 
for  imperfect  argon  the  dependence  on  Ms  becomes  strong  at  Ms  > 7 when  ionisa- 
tion starts.  The  perfect-gas  lines  (Fig.  38,  dashed  lines)  level  out  and  x 
becomes  independent  of  Ms  or  ^ = x(0w)  • This  type  of  behaviour  is  similar 
to  the  results  that  were  obtained  using  the  method  of  Law  and  Glass  (Ref.  51) 
as  shown  in  Fig.  37* 

4 .3  Domains  of  Shock-Wave  Diffractions  in  a Diatomic  Gas  - Nitrogen 

Figures  23  (with  ^ subtracted)  and  25  were  superimposed  to  obtain 
Fig.  39*  Only  the  boundary  line  corresponding  to  P0  = 15  torr  is  reproduced. 

Out  of  a maximum  of  eight  possible  shock-wave  diffractions  only  seven  are 
obtained  in  the  range  1 < Ms  < 10.  The  unobtainable  diffraction  is  an  RR 
with  an  attached  shock  wave  at  the  wedge  corner.  Note  that  if  the  lines 
9l  + S2m  = 0 and  the  attached/detached  lines  are  extrapolated  beyond  Ms  = 10 
they  might  intersect,  and  hence  an  RR  with  an  attached  shock  wave  could  be 
obtained.  The  seven  different  shock  wave  diffractions  in  the  range  1 < Mg  < 10 
appear  in  Fig.  39  and.  are  listed  in  Table  1.  They  consist  of  an  RR  with  a 
detached  shock  wave  (region  1) , SMR  having  a detached  or  attached  shock  wave 
(regions  2 and  3>  .respectively) ; OMR  with  deflection  through  a detached  or 
attached  shock  wave  (regions  4 and  5 5 respectively);  and  finally  DMR  with  a 
detached  or  an  attached  shock  wave  (regions  6 and  7,  respectively) . 

4 .4  Domains  of  Shock-Wave  Diffractions  in  a Monatomic  Gas  - Argon 

Figures  24  (with  ^ subtracted)  and  26  were  superimposed  to  obtain 
Fig.  40.  Only  the  lines  corresponding  to  P0  = 15  torr  are  reproduced.  In 
the  range  1 < Mg  < 10  this  time  only  six  diffractions  out  of  a maximum  of  eight 
are  possible . The  two  unobtainable  diffractions  are  an  RR  and  a DMR  with  an 
attached  shock  wave  at  the  corner.  However,  if  the  lines  9i  4-  02E  = 0,  = 

1.00  and  the  attached/detached  line  (dw  = dmo<)  £ne  extrapolated  beyond  Ms  = 10, 
the  two  former  lines  might  intersect  the  latter  one,  resulting  in  the  missing 
two  diffractions.  The  six  different  shock- wave  diffractions  in  the  range 
1 < Mg  < 10  are  RR  with  a detached  shock  wave  ( region  1) ; SMR  with  a detached 
or  an  attached  shock  wave  (regions  2'  and  3>  respectively);  CMR  with  a 
detached  or  an  attached  shock  wave  (regions  4 and  5j  respectively);  and 
DMR  with  a detached  shock  wave  (region  6)  The  six  different  diffractions 
are  shown  in  Fig.  40  and  listed  in  Table  2. 
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. OBLIQUE  SHOCK-WAVE  REFLECTION  IN  STEADY  FLOWS 


Prior  to  the  present  study  it  was  accepted  among  most  researchers 
that  the  shock  wave  reflection  phenomenon  in  the  nonstationary  case  is  equi- 
valent to  the  reflection  process  in  steady  flows,  once  a frame  of  reference 
is  attached  to  the  triple  point  in  SMR  or  the  reflection  point  in  RR  (Fig.  4l)  . 
However,  since  several  differences  have  already  been  pointed  out  throughout 
this  report  it  is  of  interest  to  summarize  the  phenomenon  in  steady  flow  as 
well  as  to  compare  the  steady  and  nonstationary  shock-wave  reflections . Note 
that  one  cannot  use  the  term  "diffraction"  in  steady  flows  as  only  reflection 
takes  place. 

When  a steady  supersonic  flow  encounters  a wedge  with  an  angle  9p 
(Fig.  42)  it  will  be  deflected  by  an  oblique  shock  wave  I through  an  angle  9q 
to  become  parallel  to  the  deflecting  surface.  The  angle  of  incidence  as 
well  as  all  the  thermodynamic  and  dynamic  properties  of  state  (l)  behind  the 
oblique  shock  wave  I are  defined  by  P0,  T0,  Mq  and  9i-  If  the  flow  behind  the 
oblique  shock  wave  I is  supersonic  (Mq  > 1)  a reflected  shock  wave  occurs  to 
turn  it  away  from  and  make  it  parallel  to  the  wall.  It  is  an  experimental 
fact  that  the  weak  solution  is  the  one  that  occurs.  Consequently  Mq  will  be 
greater  than  unity. 

If  the  value  of  9i  (see  Fig.  42)  is  smaller  than  the  maximum  deflec- 
tion angle  92m  of  the  flow  in  state  (l)  an  RR  occurs  (Fig.  43a) , while  for 
92  > ®2m  an  RR  is  not  possible  and  an  SMR  occurs  (Fig.  43b)  . It  is  not  assumed 
that  this  is  the  transition  criterion  from  RR  to  SMR.  Note  that  an  RR  or  an 
SMR  can  also  be  obtained  other  than  from  a rigid  wall  by  the  intersection  of 
two  oblique  shock  waves  (generated  by  two  wedges  in  a steady  supersonic  flow) 
as  shown  in  Figs . 44a  and  44b . Actual  interferograms  of  RR  and  SMR  in  a 
supersonic  wind  tunnel  are  shown  in  Figs.  45a  and  45b  (standard  interferogram) 
and  Figs.  46a  and  46b  (infinite  fringe  interferogram) . (We  are  indebted  to 
Dr.  H.  G.  Hornung  of  the  Australian  National  University  for  supplying  us  with 
these  interferograms;  see  Ref.  77*) 

It  was  noted  in  Section  2.2.3  that  the  termination  criterion  for  RR 
in  steady  flows  is  different  from  a nonstationary  flow  (Refs.  65  and  74) . 

While  in  nonstationary  flows,  the  "detachment"  criterion  or  the  "sonic" 
criterion  (Eqs.  2.l4  and  2.l6,  respectively)  is  used  (Fig.  2a,  Riv) • In  steady 
flow,  Hornung  and  Kychakoff  (Ref.  65)  verified  Henderson  and  Lczzi's  (Ref.  60) 
"mechanical- equilibrium"  criterion,  described  by  Eq.  2.15  (Fig.  2a,  Rqq) . It 
should  be  noted,  however,  that  Hornung  and  Kychakoff' s general  criterion  for 
the  transition  RR  ->  SMR  includes  both  specific  criteria  of  steady  and  non- 
stationary flows. 

5 .1  Types  of  Reflection  Obtainable 

It  is  an  experimental  fact  that  only  SMR  and  RR  have  been  obtained 
in  steady- flow  wind-tunnel  experiments.  It  is  suggested  that  the  reason  lies 
in  the  experimental  set-up  used  to  obtain  oblique  shock-wave  reflections  in 
steady- flow  experiments.  Recall  that  the  reason  for  the  termination  of  SMR 
and  formation  of  CMR  and  later  DMR  in  the  nonstationary  case  is  that  the  flow 
behind  the  reflected  shock  wave  R becomes  supersonic  wrt  the  triple  point  T 
(Mg  > 1)  while  it  is  directed  towards  the  wedge  surface  (see  discussion  in 
Section  2.4).  In  the  steady  case,  the  flow  in  state  (2)  is  always  interrupted 
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by  either  an  expansion  wave  or  a reflected  shock  wave  (Figs.  47a  and  47b, 
respectively) , which  arises  from  the  wedge  used  to  generate  the  incident 
oblique  shock  wave  I.  Thus  the  flow  in  state  (2)  behind  the  reflected 
shock  wave  is  prevented  from  reaching  the  w’edge  surface  again.  Therefore 
the  need  to  turn  this  flow  away  from  the  wall  is  eliminated.  Consequently, 
in  steady  flows  unlike  nonstatiopary  flows  it  is  possible  to  have  SMR  with 
supersonic  flow  behind  the  reflected  wave. 

5*2  Steady  Oblique  Shock  Wave  Reflection  in  (M_,  0q) -Plane 

In  steady  flows  the  parameters  Ms  and  Qy,  which  were  used  for 
nonstationary  flows  cease  to  be  significant.  Instead  the  flow  Mach  number 
in  state  (0),  Mo  and  the  flow  incidence  angle  0O  to  the  oblique  shock  wave  I 
are  used.  Consequently  the  reflection  phenomenon  is  described  in  the 
(Mo,  0O) -plane. 

5*2.1  Steady  Flow  Regions  and  Transition  Boundaries  in  a Diatomic  Gas  - Nitrogen 

The  domains  of  different  types  of  oblique  shock  wave  reflection  in 
nitrogen  are  shown  in  Fig.  48.  The  solid  boundary  lines  account  for  real-gas 
effects,  while  the  dashed  lines  are  for  a perfect  diatomic  gas  (7  = 1.4) . 

The  (Mo,  0o)-plane  is  divided  into  3 domains  of  RR,  SMR  and  NR.  The  signi- 
ficance of  real-gas  effects  on  shifting  the  boundary  lines  can  be  clearly 
seen.  The  perfect-gas  theory  can  be  used  only  in  the  narrow  range  1 < Mo  < 2.22. 

At  = 2.22  the  imperfect-gas  lines  start  to  diverge  from  the  perfect-gas  line 
due  to  vibrational  excitation.  At  higher  Mach  numbers  (Mo  > 8.6)  the  imperfect- 
gas  line  splits  into  a multiplicity  of  lines  for  different  initial  pressures  at 
constant  temperature  due  to  dissociation.  At  even  higher  Mach  numbers  electronic- 
excitation  and  ionization  would  play  a similar  role.  When  Mq  approaches  unity 
the  line  01  4-  82  = 03  = 0 approaches  the  line  Mq  = 1.00.  These  two  lines  are 
coincident  at  Mo  = 1 and  0o  = 90*. 

5*2.2  Steady  Flow  Regions  and  Transition  Boundaries  in  a Monatomic  Gas  - Argon 

The  domains  of  SMR,  RR  and  NR  for  argon  are  shown  in  Fig.  49.  The 
solid  lines  are  again  for  imperfect  argon  while  the  dashed  lines  are  for  a perfect 
monatomic  gas  (7  = 5/3) • 


Unlike  the  case  of  a diatomic  gas  where  the  perfect-gas  theory  was 
adequate  only  for  very  small  values  of  Mq  (Mo  < 2.22,  Fig.  48),  here  (Fig.  49) 
it  can  be  used  over  a much  wider  range,  1.00  < Mo  < 8-95*  At  Mo  = 8.95  the 
imperfect-gas  boundary  line  splits  with  initial  pressure  (at  constant  tempera- 
ture) due  to  electronic  excitation  and  ionization.  Consequently  one  can  expect 
an  SMR  when  Mo  = l6,  0O  = 30°  and  Po  > 10  torr,  while  an  RR  will  be  obtained  if 
Po  < 1 torr. 

5 .3  Comparison  Between  Steady  and  Nonstationary  Shock  Wave  Reflection 

It  has  already  been  mentioned  that  in  spite  of  the  similarity  between 
shock-wave  reflection  in  steady  and  nonstationary  flows,  several  differences 
exist  between  them.  Therefore  it  is  worthwhile  to  summarize  these  differences. 
The  first  and  most  important  difference  lies  in  the  fact  that  the  shock-wave- 
reflection  process  in  the  nonstationary  case  always  interacts  with  the  flow 
deflection  process  over  the  comer,  while  in  steady  flows  only  reflection  occurs. 
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Consequently,  unlike  that  in  steady  flow,  an  NR  domain  exists  in  the  (Mo,  0o)- 
plane  (Figs.  48  and  49);  in  the  nonstationary  case,  there  is  no  actual  NR 
domain  (Figs.  39  and.  40)  in  the  (Ms,  Qy) -plane.  Recall  that  in  the  nonstationary 
case,  a reflection  occurs  even  when  the  flow  behind  the  incident  shock  wave  I 
becomes  subsonic  wrt  point  P (Fig.  4a)  due  to  the  interaction  with  the  flow- 
deflection  process  (Fig.  4b).  In  the  steady  flow  case  only  one  reason  exists 
for  shock-wave  reflection,  i.e.,  a need  to  turn  the  supersonic  flow  moving 
towards  a solid  surface  (Fig.  42).  In  the  nonstationary  case  there  are  two 
reasons . The  first  is  analogous  to  the  steady  case  and  holds  for  Mq  > 1 
(Fig.  4a),  while  the  other  applies  to  Mq  < 1 (Fig.  4b).  Note  that  in  steady 
flows  Mq  < 1 yields  no  reflection  (NR).  A more  detailed  discussion  of  these 
two  reasons  is  given  in  Section  2.1. 

In  Section  2.2.3  it  was  noted  that  the  termination  criteria  of  RR 
are  different  for  steady  and  nonstationary  flows.  In  nonstationary  flows  either 
the  "detachment"  criterion  (Eq.  2.14)  established  by  Neumann  (Refs.  2,  4)  or  the 
"sonic"  criterion  (Eq.  2.16)  established  by  Hornung  and  Kychakoff  (Ref.  65)  holds, 
while  in  steady  flow  the  "mechanical-equilibrium"  criterion  (Eq.  2.15),  established 
by  Henderson  and  Lozzi  (Ref.  60)  and  Hornung  and  Kychakoff  (Ref.  65),  applies. 

Once  RR  is  terminated  three  different  types  of  reflection  can  occur 
in  nonstationary  flows,  i.e.,  SMR,  CMR  and  DMR.  In  steady  flows,  however,  CMR 
and  DMR  cannot  occur  due  to  the  geometrical  limitations  used  to  obtain  the 
reflection  (see  discussion  in  Section  5-1)  and  hence  only  an  SMR  is  possible. 

It  should  be  noted,  however,  that  while  in  steady  flow  one  can  have  an  SMR 
with  supersonic  flow  behind  the  reflected  shock  wave,  in  nonstationary  flows 
SMR  is  terminated  when  the  flow  behind  the  reflected  wave  becomes  supersonic. 

The  interaction  between  the  shock-wave- reflection  process  and  the 
flow-deflection  process  in  nonstationary  flow  results  in  a three- shock  con- 
fluence where  the  reflected  wave  R,  Mach  stem  M and  slipstream  S are  usually 
curved  (not  in  the  vicinity  of  the  triple  point  where  they  are  all  straight, 
see  Ref.  53).  In  steady  flow,  these  shock  waves  and  slipstream  are  usually 
straight,  and  hence  the  flow  regions  bounded  by  them  are  much  more  uniform 
than  those  in  the  nonstationary  case.  One  should  note  that  in  steady  flow 
the  slipstream  never  curls  back,  a typical  property  of  nonstationary  reflec- 
tions . 

Due  to  the  fact  that  the  incident  shock  wave  in  nonstationary  reflec- 
tions always  moves  into  a quiescent  flow,  the  reflection  point  (Fig.  4a)  is 
always  independent  of  the  boundary-layer  growing  behind  it.  However,  in  steady 
flow  the  boundary-layer  effects  depend  on  the  geometrical  set-up  used  to  obtain 
the  reflection.  A severe  boundary- layer  interaction  is  eliminated  if  the  inter- 
section of  shock  waves  is  used  to  obtain  the  reflection  (Figs.  44a  and  44b). 

The  triple-point  trajectory  angle  x has  no  physical  meaning  in  steady 
flows,  where  the  configuration  is  fixed  wrt  time,  i.e.,  does  not  grow  linearly 
with  time.  However,  in  steady  flow  the  length  of  the  Mach  stem  can  be  regarded 
as  a similarity  parameter.  So  far  a theory  for  predicting  the  length  of  the 
Mach  stem  has  not  been  published.  The  equations  of  motion  (Eqs.  2.17  to 
2.30)  are  all  satisfied  if  one  moves  the  triple  point  to  any  point  along  the 
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incident  shock  wave  (Fig.  50)  and  draws  the  corresponding  reflected  shock 
wave  R*,  Mach  stem  M*  and  slipstream  S*.  However,  repeated  experiments  for 
the  same  initial  conditions  and  geometrical  configurations  provide  identical 
shock-wave  configurations.  This  suggests  that  the  factors  governing  the 
reflection  process  in  steady  flow  are  related  to  the  waves  emanating  from  the 
wedge  (Figs.  47a  and  47b).  It  should  be  noted  here  that  Molder  (Ref.  76) 
developed  a theory  for  the  reflection  of  curved  shock  waves  in  steady  flows, 
using  the  assumption  of  constant  curvature,  by  which  the  length  of  the  Mach 
stem  can  be  determined  analytically.  However,  his  theory  fails  to  predict  the 
length  of  the  Mach  stem  in  the  case  of  the  reflection  of  straight  shock  waves. 
(Owing  to  the  lack  of  experimental  data  on  the  reflection  of  curved  shock  waves, 
Molder’ s analysis  remains  unsubstantiated.) 

The  incident  shock  wave  was  observed  to  be  accelerating  somewhat  when 
it  passed  over  the  compression  model  in  the  test  section.  A constant-velocity 
shock  wave  was  measured  when  the  model  was  removed.  The  reason  for  this  accelera- 
tion is  not  known.  Although  this  acceleration  is  usually  very  small  it  might  have 
an  influence  on  the  reflection  process.  Consequently  Eqs.  2.1  to  2.9  for  an  RR 
and  Eqs.  2.17  to  2.30  for  the  triple  point  T describe  the  reflection  phenomenon 
in  steady  flows  better  than  they  do  in  nonstationary  flows,  where  the  flow  Mach 
nuntoer  Mq  ahead  of  the  incident  shock  wave  [state  (0) ] changes  slightly  as  the 
incident  shock  wave  accelerates,  i.e.,  Ms  increases. 


6.  EXPERIMENTAL  TECHNIQUES  AND  INSTRUMENTATION 

The  UCIAS  10  cm  x 18  cm  Hypervelocity  Shock  Tube  (Ref.  32)  was  used 
to  generate  the  incident  shock  waves.  An  updated  and  detailed  description  of 
the  entire  facility,  its  performance  maintenance  and  operational  techniques  is 
reported  in  Ref.  71.  Therefore  in  the  following  only  a "brief  discussion  is  given. 

The  cold  runs  were  used  to  generate  the  shock  waves  in  the  range  2 < 

Ms  < 8 of  the  present  study.  The  initial  conditions,  i e.,  the  pressure  ratio 
across  the  mylar  diaphragm  and  the  driver  gas  combinations  used  to  obtain 
different  incident  shock  waves  into  nitrogen  and  argon  are  given  in  Tables  3 
and  4,  respectively  The  shock  wave  Mach  numbers,  Ms,  given  in  Tables  3 and 
4 represent  the  average  values  obtained  throughout  the  present  study.  Conse- 
quently, the  initial  conditions  listed  in  Tables  3 and  4 will  produce,  in 
general,  shock  waves  in  the  range  Ms  ± 0.1  of  the  predicted  value.  The  initial 
pressures  P0  in  the  range  5-00  < PQ  < 65. 00  torr  were  monitored  with  two  Wallace 
and  Tiernan  type  FA  160  (0-50  torr  and  0-200  toir)  dial  gauges.  A more  accurate 
measurement  of  P0  was  made  using  an  oil  manometer  (0.2-40  torr)  . 

The  nonstationary  process  was  recorded  using  a 23-cm  dia  field  of 
view  Mach-Zehnder  interferometer  (Refs  21,  71)  The  light  source  was  a 
giant-pulse  ruby  laser  (Ref  71)  Simultaneous  dual- frequency  interferograms 
were  taken  at  6943A  and  3471. 5A,  respectively.  Kodak,  Royal  X Pan  (1250  ASA), 

10  x 13  cm  (4  x 5 inch)  plate  film  was  used  to  record  the  interferograms. 

A schematic  instrumentation  diagram  of  the  equipment  used  throughout 
the  present  study  is  shown  in  Fig.  51-  When  the  shock-wave  arrives  at  station 
D,  time-counters  1,  2,  3>  4 and  5 are  triggered  simultaneously.  (Note  that 
these  5 counters  are  connected  to  a common  external  1-megacycle  oscillator, 
and  hence  have  a uniform  time  base  of  1-microsecond  resolution.)  These  5 
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counters  then  stopped  in  succession  as  the  shock  wave  passed  subsequent 
detectors  at  stations  F,  G,  H,  I and  J.  An  additional  counter  (No.  7)  was 
used  between  stations  I and  J to  indicate  the  time  taken  for  the  shock  wave 
to  travel  past  the  test  section.  Since  the  shock  wave  velocity  changes 
somewhat  due  to  the  convergence  of  the  cross  section  caused  by  the  model 
(Refs.  39  and.  46),  the  measurement  of  counter  No.  7 was  used  only  as  a monitor 
of  the  shock  location. 

The  control  of  the  laser  light  source  operation  is  also  included  in 
Fig.  51.  The  laser  flashlamp  capacitor  bank  was  usually  triggered  by  the  shock 
arrival  at  station  F,  suitably  delayed  in  a Tetronix  type  555  oscilloscope  to 
get  900  ^sec  of  energy  punning  by  the  time  the  Pockels-cell  shutter  opened. 

The  Pockels  cell  Q,- switch  was  always  triggered  from  station  I,  delayed  correctly 
in  a pulser  delay  unit  to  take  the  interferograms  at  a desired  time  (or  location 
of  the  shock  wave  wrt  station  I) . 

To  monitor  the  actual  flashlamp  pumping  time  in  each  experiment, 
synchronization  pulses  from  both  the  flashlamp  capacitor  bank  and  the  Pockels- 
cell  pulser  unit  were  used  to  start  and  stop,  respectively,  a UTIAS  microsecond 
counter  (Fig.  51 > Counter  No.  6) . 

As  an  additional  check  on  shock- tube  performance  in  each  experiment, 
a Kistler  type  6oi-B,  piezo-electric  pressure  transducer  was  used  to  monitor 
the  pressure  variation  with  time  behind  the  incident  shock  wave.  The  transducer 
output  was  taken  directly  to  a charge  amplifier  (Kistler,  Model  504)  with  a 
final  signal  displacement  on  the  designated  oscilloscope  (Fig.  51) • 

6.1  Accuracy  of  Measurements 

A detailed  description  of  the  absolute  errors  associated  with  the 
various  measurements  involved  in  the  operation  of  the  shock  tube  facility 
is  given  in  Ref.  71-  Consequently,  only  a brief  summary  follows. 

6.1.1  Incident  Shock  Wave  Mach  Number  - M 

s 

The  maximum  possible  relative  error  in  calculating  the  incident 
shock  wave  Mach  number  was  (Ref.  71): 

E(M  ) 

— „ S"  = (1.15  M + 10 .18)  x 10  for  nitrogen  (6.1) 

M S 

S 

E(MS}  , -3 

-jj = (1.05  x Ms  + 10.18)  x 10  J for  argon  (6.2) 


Consequently,  in  the  range  of  incident  shock  wave  Mach  numbers  used  in  the 
present  study,  2 < Mg  < 8,  E(MS)  = 0.02  for  both  nitrogen  and  argon  at  the 
lower  end  (Ms  = 2j  while  for  Mg  = 8,  E(Mg)  = 0.l6  and  0.15  for  nitrogen  and 
argon,  respectively.  Equations  6.1  and  6.2  are  shown  in  Fig.  52a.  It  can  be 
seen  that  the  relative  error  E(Ms)/Ms  increases  linearly  as  Ms  increases. 
Consequently  the  absolute  error  E(Mg)  is  proportional  to  Mg2. 


6.1.2  Initial  Pressure  PQ  and.  Temperature  Tq 

Initial  pressures  in  the  range  5 < Po  < 40  were  measured  with  an  oil 
manometer.  The  pressure  PQ  is  calculated  from  (Ref.  71): 


Pen  H(nm) 
13-5951 


(6.3) 


where  H is  the  oil-height  difference  in  mm  as  measured  in  the  manometer,  and 
Poii  is  the  density  of  the  oil  and  it  is  given  by: 


Poii  = 1-0690  + 9-5  x IQ"  (25-T) 


(6.4) 


I, the  .oil  tape rature,.  was  assumed  to  be  equal  to  the  laboratory  temperature 
measured  near  the  manometer  (T  = Tp) • The  absolute  error  associated  with  the 
initial  pressures  is  then  (Ref.  71): 


E(P  ) = [1.4  x 10~^  H(mm)  + 7.86  x 10  2]  torr 


(6.5) 


Consequently,  the  maximum  possible  errors  associated  with  lowest  (5  torr)  and 
highest  (37  torr)  pressures  that  were  measured  with  the  oil  manometer  are  0.08 
torr  and  0.09  torr,  respectively. 

Pressures  in  the  range  PQ  > 40  torr  were  measured  with  a Wallace 
and  Tiernan  type  FA  160  (0-200  torr)  within  an  accuracy  of  ±0.2  torr  and 
hence  the  maximum  possible  error  was  E(P0)  =0.4  torr.  A plot  of  the  absolute 
error  involved  in  the  pressure  measurements  in  the  present  study  5 < Po  < 65 
is  given  in  Fig.  52b. 

The  initial  terperature  To  was  measured  to  an  accuracy  of  ± 0.1°  with 
a standard  mercury  bulb  thermometer  inserted  in  a 5 cm  deep,  6 mm  dia  oil-filled 
port  in  the  upper  wall  of  the  shock  tube.  Consequently,  an  error  E(Sk)  =0.2° 
was  associated  with  all  temperature  readings. 


7-  EXPERIMENTAL  RESULTS  AM)  DISCUSSION 

The  above  discussed  facility  was  used  to  obtain  and  record  non- 
stationary oblique  shock-wave  diffraction  in  nitrogen  and  argon.  Fifty-eight 
successful  experiments  were  made  in  nitrogen,  with  9 different  compression 
corners  (wedges)  having  0w  = 2*,  5*}  10*,  20°,  26.56",  30",  40",  50*  and  60". 
In  the  case  of  argon  48  successful  experiments  were  performed  with  7 different 
wedge  angles, 0w  = 2",  10",  20",  30",  40",  50*  and  60*.  All  angles  had  an 
accuracy  of  ± 0 .0167*. 

Unfortunately,  due  to  mechanical  limitations  wedge  angles  lower  than 
2°  could  not  be  machined.  Wedge  angles  greater  than  60*  were  not  of  interest 
as  they  would  have  resulted  in  RR,  thereby  not  adding  any  inportant  data.  The 
incident  shock  wave  Mach  numbers  Mg,  wedge  angle  and  initial  pressures  and 


and  temperatures  (P0  and  To)  associated  with  each  experiment  are  listed  in 
Table  5 for  nitrogen  and  Table  6 for  argon.  These  tables  also  include  the 
experimental  values  of  the  primary  and  secondary  triple  point  trajectory 
angles  x and-  x'  as  well  as  the  type  of  reflection  obtained  with  the  initial 
conditions.  Note  that  x'  is  listed  also  for  CMS;  it  is  omitted,  however, 
for  CMR  where  the  location  of  the  kink  is  not  sufficiently  clear.  Since 
the  location  of  the  kink  in  a CMR  is  not  as  well  defined  as  the  location  of 
the  second  triple  point  in  a DMR , the  accuracy  of  measuring  x'  in  a CMR  is 
smaller  than  in  a DMR.  Consequently,  the  absolute  errors  are  E(xO  = 1*  for 
CMR  and  E(x')  = 0.5®  for  DMR.  (The  errors  associated  in  measuring  Mg,  Po 
and  To  are  discussed  in  Chapter  6.) 

Owing  to  the  size  of  the  present  report,  most  of  the  interferograms 
have  been  omitted.  However,  as  they  are  inportant  as  a data  bank  that  might 
be  helpful  in  the  future  to  various  researchers,  all  the  actual  interfero- 
grams of  the  nonstationary  oblique  shock  wave  reflections  that  were  obtained 
throughout  the  present  study  (see  Tables  5 and  6 for  details)  are  given  in 
Ref.  72.  In  all,  106  experiments  were  conducted. 

7.1  Verification  of  Nonstationary  Oblique  Shock-Wave  Reflection  in  (M_,  e')- 
Plane  S W 


In  the  present  analysis  the  two  different  phenomena  of:  shock- wave 
reflection  by  the  wedge  surface  and  the  induced-flow  deflection  over  the  wedge 
corner  were  treated  independently  (in  Chapters  2 and  3,  respectively) . In 
this  manner  Figs.  23  (nitrogen)  and  24  (argon)  were  produced  for  reflection 
and  Figs.  25  (nitrogen)  and  26  (argon)  for  flow  deflection.  A superposition 
of  these  figures  (Chapter  4)  produced  the  overall  phenomena  of  oblique  shock- 
wave  diffractions  shown  in  Figs.  39  and-  40  for  nitrogen  and  argon,  respectively. 
Recall  that  the  appropriate  value  of  x was  subtracted  from  Figs . 23  and  24  in 
order  to  superimpose  them  with  Figs.  25  and  26. 

In  the  following,  the  present  analysis  of  nonstationary  reflection  of 
oblique  shock  waves  is  compared  with  experiments.  Figures  53  and  54  are 
reproductions  of  Figs.  23  and  24  with  Po  = 15  torr  only,  for  clarity.  Note 
that  although  most  of  the  present  experiments  were  made  with  PQ  ~ 15  torr, 
those  corresponding  to  low  Mach  numbers  (Mg  « 2 in  nitrogen  and  argon)  and 
high  Mach  numbers  (Ms  ^ 6.79  in  nitrogen  and  Ms  ^ 7 >53  in  argon)  were  obtained 
with  slightly  different  initial  pressures  (see  Tables  5 and  6 for  details) . 
However,  Figs.  23  and  24  indicate  that  the  initial,  pressure  becomes  signi- 
ficant for  the  boundary  lines  Ml  = 1.00,  Vo  = 1.00  and  Mg'  = 1.00  (Fig.  23) 
at  Ms  = 7.56,  7-51  and  7*47,  respectively,  and  for  the  boundary  lines  Mq  = 1.00, 
= 1.00  and  ^2'  = 1.00  (Fig.  24)  at  Ms  = 7.72,  7.65  and  7-32,  respectively. 
These  high  values  of  Ms  were  about  the  highest  used  throughout  the  present 
study.  For  the  boundary  line  9i  + Qg™  = 0 of  Figs.  23  and  24  the  initial 
pressure  becomes  significant  at  Ms  = 6.05  and  5-56,  respectively.  However, 
in  the  neighbourhood  of  these  boundary  lines  the  experiments  did  not  go  up 
to  Ms  ~ 8,  since  in  similar  experiments  that  were  conducted  in  the  past 
with  the  same  facility  the  high-qjiality  optical  windows  were  burned  (Ref.  43) . 

7.1.1  Diatomic  Gas 


The  experimental  data  of  Smith  (Ref.  6)  for  Ms  = 1.25,  1.51,  2.10  and 
2.40,  and  White  (Ref.  17)  in  air,  Law  and  Glass  (Ref.  51)  in  oxygen,  and  our 
data  in  nitrogen  are  all  added  to  Fig.  53  (a  reproduction  of  Fig.  23  for 
Po  = 15  torr)  in  order  to  check  the  present  analysis  for  nonstationary 
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oblique  shock-wave  reflection  in  diatomic  gases.  [Note  that  the  term  reflec- 
tion is  used  here  as  the  results  are  plotted  in  the  pseudo-steady  plane 
(Ms,  04).]  Two  points  from  Smith's  experimental  data  that  are  reported 
to  be  RR  lie  slightly  below  their  analytical  terminating  boundary  line. 

We  believe  that  this  slight  disagreement  is  probably  due  to  Smith's  use  of 
air  (the  lines  in  Fig.  53  are  for  nitrogen)  . The  presence  of  20 Jo  oxygen  in 
air  would  cause  the  boundary  line  to  be  lower  than  that  of  nitrogen,  since 
the  oxygen  becomes  vibrationally  excited  at  a lower  tenperature,  and  hence 
the  line  would  start  shifting  downwards  at  a lower  Mach  number. 

Shock-wave  reflection  configurations  reported  by  Smith  (Ref.  6) 
to  be  SMR  in  the  range  2.10  < Ms  < 2.42  appear  to  lie  in  the  region  that 
corresponds  to  CMR  (Fig.  53)-  However,  a careful  check  of  Smith's  report 
reveals  that  he  observed  that,  "for  strong  shocks  (his  strongest  shock  was 
Ms  = 2.42)  a reversal  of  curvature  (in  the  reflected  shock  wave)  develops", 
and  furthermore  "the  portion  of  the  reflected  shock  near  the  triple  point 
that  appears  to  be  straight".  This  we  believe  corresponds  to  CMR.  It  is 
clear  that  although  he  had  noticed  a CMR- configuration , he  referred  to 
it  as  an  SMR  rather  than  propose  a new  type  of  reflection,  since  these  two 
configurations,  except  for  the  kink  in  R,  are  quite  similar  in  appearance. 

When  White  (Ref.  17)  discovered  DMR  in  1951,  the  importance  of  CMR  was 
recognized  as  a different  type  of  reflection.  However,  in  White's  report 
CMR  is  still  considered  as  SMR.  Note  that  one  experimental  point  of  White 
(Ref.  17),  Mg  = 2.3,  0w  = 39*2*  in  the  region  where  Smith  (Ref.  6)  reported 
SMR,  resulted  in  a CMR;  this  verifies  the  previous  remark  that  those  experi- 
ments by  Smith  are  CMR  and  not  SMR.  Excellent  agreement  can  be  seen  with 
all  the  other  boundary  lines  since  all  the  experimental  points  lie  in  their 
predicted  regions. 

7.1.2  Monatomic  Gas 

The  present  experimental  results  (Table  6)  in  argon  as  well  as 
some  data  from  Law  and  Glass  (Ref.  51)  in  argon  and  helium  are  all  shown  in 
Fig.  54  (a  reproduction  of  Fig.  24,  for  PQ  = 15  torr) . 

One  experiment  reported  by  Law  and  Glass  as  CMR  lies  inside  the 
DMR  regime  (Ms  = 7-66,  9^  r»45°).  However  this  experimental  point  is  sur- 
rounded by  DMR  points  that  were  obtained  during  the  present  study.  It  might 
be  that  the  initial  directions  of  the  fringes  chosen  by  Law  and  Glass  (differ- 
ent from  ours)  was  such  that  the  Rq-shock  wave  (Fig.  19)  of  a DMR  could  not 
be  seen  clearly.  Note  that  Bazhenova  et  al  (Ref.  62)  also  reported  DMR  in 
argon,  in  the  region  where  Law  and  Glass  did  not  notice  it. 

Out  of  all  our  experimental  points,  only  one  (Ms  = 4.44,  9^  = 45.5*) 
that  corresponds  to  a DMR  lies  outside  its  predicted  region,  in  the  CMR 
domain.  This  we  believe  is  due  to  the  use  of  Eq.  2.47  in  calculating  the 
CMR/DMR  boundary  line.  Recall  that  Eq.  2.47  was  found  by  Law  and  Glass 
(Refs.  48  and  51)  and  Bazhenova  et  al  (Ref.  46)  to  be  in  good  agreement  with 
experiments  only  in  the  range  9w  < 40*  while  for  the  range  9w  > 4o*  the 
agreement  becomes  progressively  worse.  Note  that  aside  from  this  slight 
disagreement  all  the  other  experimental  points  lie  inside  their  predicted 
regions . 

7.2  Triple  Point  Trajectory  Angle  - % Comparison  With  Experiments 

In  order  to  superimpose  the  shock- reflection  process  (Figs.  23 
for  nitrogen  and  24  for  argon)  with  the  flow- deflection  process  (Figs.  25 
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and  26  for  nitrogen  and  argon,  respectively)  the  value  of  x was  subtracted 
from  the  corresponding  boundary  lines  of  Figs.  23  and  24.  Two  different 
methods  for  the  prediction  of  x were  discussed  in  Section  2. 3. 1.1.  The 
first,  due  to  Law  and  Glass  (Ref.  51) , was  found  by  them  to  be  in  good 
agreement  with  experiments  (in  oxygen)  only  in  the  range  25*  < 9w  < 45®. 

The  other,  developed  here,  is  for  very  small  wedge  angles, 

9w  < 5°*  These  two  methods  were  applied  to  nitrogen  and  argon  (Sections 
4.1  and  4.2,  respectively)  to  obtain  Figs.  31  and  32  for  nitrogen,  and 
Figs.  37  and  38  for  argon.  In  the  following,  these  two  methods  are  checked 
against  the  experimental  results  in  nitrogen  and  oxygen. 

7-2.1  Diatomic  Gas 

The  present  data  (Table  5)  were  added  to  Fig.  55  (a  reproduction 
of  Fig.  31),  in  order  to  test  the  present  analytical  version  of  the  graphical 
method  by  Law  and  Glass  (Ref.  51)  for  predicting  x ln  nitrogen.  It  is  seen 
that  the  agreement  with  experiments  for  wedge  angles  in  the  range  5*  < 9w  < 4o® 
is  reasonably  good.  Most  of  the  experimental  points  in  this  range  lie  within 
± 0.5*  from  their  predicted  values.  However,  for  Gy  < 5*  the  agreement  becomes 
progressively  worse  as  decreases.  Note  that  while  Law  and  Glass  (Ref.  51), 

who  developed  this  method  of  predicting  x>  found  it  to  be  good  only  in  the 
range  25°  < 0W  < 45°  (in  oxygen) , it  is  seen  that  good  agreement  actually 
exists  over  a much  wider  range  5*  < 9w  < 45®.  This  is  probably  due  to  the 
fact  that  they  solved  Eqs.  2.17  to  2.31  graphically,  whereas  the  present 
solution  is  analytical.  During  the  present  study  their  graphical  solution 
was  checked  for  nitrogen.  It  was  found  to  exhibit  the  same  characteristics 
as  their  graphical  solution  for  oxygen.  Note  that  for  < 40*  the  actual 
values  of  x are  greater  than  those  predicted.  In  this  range  all  the  data 
points  fall  between  the  predicted  perfect  and  imperfect  gas  models. 

The  experimental  data  for  2®  and  5®  wedge  angles  (Table  5)  are 
shown  in  Fig.  56  (a  reproduction  of  Fig.  32)  in  order  to  check  the  present 
method  for  predicting  x for  Tow  wedge  angles.  It  can  be  seen  that  for  a 5" 
wedge  the  present  analytical  version  of  the  graphical  method  of  Law  and  Glass 
(Ref.  51)  is  actually  in  better  agreement  with  experiments.  results  than  the 
present  method  for  small  wedge  angles.  For  9w  < 4*  the  present  method  is  in 
good  agreement  with  experiments. 

It  should  be  noted  that  it  is  very  difficult  to  predict  the  value 
of  Xj  using  the  analytical  version  of  the  method  of  Law  and  Glass  for  very 
small  wedge  angles  (ow  < 5°)  owing  to  severe  convergence  problems  as  Bw 
becomes  small.  The  graphical  solution  is  possible,  however,  accuracy  is  very 
low  and  the  error  at  even  greater  wedge  angles  (5*  < 0^  < 20°)  becomes  signi- 
ficant (Ref.  51).  Consequently  this  is  the  only  existing  method  of  predicting 
X at  low  wedge  angles  (gw  < 5°) • 

7.2.2  Monatomic  Gas 

The  data  from  Table  6 for  0w  = 10®,  20®,  30° , 40®  and  50®  is  shown 
in  Fig.  57  (a  reproduction  of  Fig.  37)  in  order  to  test  the  present  analytical 
method  (based  on  Law  and  Glass,  Ref.  51)  for  predicting  x against  experiments 
for  argon.  Very  good  agreement  can  be  seen  with  the  wedge  angles  20®,  30®  and 
40°  while  for  10°  and  50®  the  actual  value  of  is  smaller  than  the  predicted 
one  by  more  than  1®  (the  error  bar) . For  these  two  wedge  angles  the  agreement 
is  fairly  good  for  small  values  of  Ms  (Ms  < 3 for  0„  = 50°  and  Ms  < 5 for 
0W  = 10° ) and  becomes  progressively  worse  as  Ms  increases . 
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The  data  points  for  ew  = 2*  (Table  6)  are  shown  in  Fig.  58  (a  reproduc- 
tion of  Fig.  38  for  9W  = 2*).  The  agreement  is  fairly  good  for  low  and  high 
Mach  numbers  while  in  the  range  3<Mg<6itis  not  so  good.  Here  as  well, 
owing  to  convergence  problems,  it  was  not  possible  to  get  a prediction  of  x 
using  our  analytical  version  of  the  method  of  Law  and  Glass  (Ref.  51)  for 
9w  = 2°. 


7.3  Verification  of  Nonstationary  Oblique  Shock  Wave  Reflections  in  (M  , 0 )- 
Plane  s w 

Figures  53  (nitrogen)  and  5*4  (argon)  correspond  to  nonstationary 
oblique  shock  wave  reflection  in  the  (Ms,  0w)-plane.  However,  while  the 
exact  value  of  the  effective  wedge  angle  9^,  can  be  measured  only  after  an 
actual  experiment  is  recorded,  the  actual  wedge  angle  Qy  is  a given  parameter. 
Consequently,  it  is  of  interest  to  obtain  the  reflection  phenomena  in  the 
(Ms,  6w) -plane  rather  than  the  (Mg,  9w) -plane.  This  is  done  by  subtracting 
the  appropriate  value  of  x (Figs.  31  and  32  for  nitrogen,  and  Figs.  37  and  38 
for  argon)  from  the  vertical  axis  of  Figs.  53  and  5*+  (nitrogen  and  argon, 
respectively) . 

The  line  0p  + 0gm  (Figs.  53  and  5*0  represents  the  boundary  of  the 
RR  process.  Consequently,  the  actual  value  of  x above  (and  on)  this  line  is 
zero.  Therefore  it  does  not  change  throughout  the  transformation.  Recall 
that  the  present  method  for  predicting  x for  very  small  wedge  angles  (Section 

2. 3-  1.1)  is  based  on  the  experimental  fact  that  in  the  nonstationary  case, 
the  incident  shock  wave  always  reflects.  Following  this  method  (Eq.  2.33) 
one  can  immediately  see  that  the  transformation  0w  -*  9w  eliminate  the  lines 
Ml  = 1.00  of  Figs.  53  and  5k.  In  other  words,  on  the  lines  Ml  = 1.00  (Figs. 

5^  and  5*+) , 0w  = X 811(1  since  0w  = 0w  + X by  definition  the  transformation 

0W  0w  yields  that  the  lines  Mj  = 1.00  coincide  with  the  lines  0w  = 0.  Con- 
sequently only  the  lines  M2  = 1.00  and  = 1.00  of  Figs.  53  and  5*4  need  be 
transformed. 

7.3- 1  Diatomic  Gas 

The  data  from  Smith  (Ref.  6),  White  (Ref.  17)  and  Bazhenova  et  al 
(Ref.  55)  in  air,  Law  and  Glass  (Ref.  51)  in  oxygen,  Bazhenova  et  al  (Ref.  55) 
and  the  present  in  nitrogen  are  all  shown  in  a combined  plot  on  Fig.  59  (ob- 
tained from  Fig.  53  by  subtracting  x)  in  order  to  verify  the  present  analysis 
of  nonstationary  oblique  shock-wave  reflection  in  the  (Mg,  0w) -plane  for 
diatomic  gases. 

Two  experiments  (present  one  and  one  fwm  Law  and  Glass)  with  0w  = *40* 
and  Mg  = 2.00  and  2.56,  respectively,  in  the  region  where  Smith  (Ref.  6)  re- 
ported SMR  results  in  CMR,  and  verifies  our  previous  remarks  that  those 
experiments  by  Smith  are  CMR  and  not  SMR.  Many  more  data  points  from  White 
(Ref.  17)  could  have  been  used.  However,  since  they  all  fall  in  their 
predicted  regions,  they  were  omitted  for  clarity.  It  can  be  concluded  from 
this  figure  that  our  analysis  of  the  nonstationary  oblique  shock-wave  reflec- 
tion in  the  (Mg,  0w) -plane  for  a diatomic  gas  is  substantiated. 

7-3-2  Monatomic  Gas 

The  data  of  Law  and  Glass  (Ref.  51)  in  helium,  Bazhenova  et  al 
(Ref.  55) , Law  and  Glass  (Ref.  51)  and  the  present  experiments  in  argon 
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appear  in  Fig.  60  (obtained  from  Fig.  54  by  subtracting  in  order  to  check  the 
present  analysis  for  nonstationary  oblique  shock-wave  reflection  in  a monatomic 
gas  on  the  (Mg,  ew)-plane.  Very  good  agreement  was  obtained  with  the  calculated 
regions  and  their  boundary  lines.  All  the  experimental  points  lying  outside 
their  predicted  regions  are  discussed  in  Section  7 .1.2.  It  can  be  concluded 
that  the  present  analysis  for  nonstationary  reflection  of  oblique  shock-waves 
in  a monatomic  gas  is  substantiated. 

7.4  Verification  of  Nonstationary  Oblique  Shock-Wave  Diffraction  in  (M_,  9;J) -Plane 

Figures  6l  and  62  are  reproductions  of  Figs.  39  and  40  for  nitrogen  and 
argon  respectively.  Unfortunately,  not  all  the  data  shown  in  Figs.  53 , 54,  59 
and  60  could  be  used  since  information  concerning  the  flow  deflection  process 
is  usual  1 y not  reported  [e.g.,  Bazhenova  et  al  (Ref.  55) , where  in  all  the 
recorded  pictures  the  wedge  corner  is  not  in  the  field  of  view].  Consequently, 
the  present  data  in  nitrogen  and  argon  as  well  as  a few  data  points  from  different 
sources  are  shown  in  Figs.  6l  and  62,  that  test  the  present  analysis  of  nonstationary 
oblique  shock- wave  diffraction  (shock  reflection  and  flow  deflection)  in  the  (Mg,  9w) ■ 
plane . 

7.4.1  Diatomic  Gas 

Our  experimental  data  points  from  Table  5 in  nitrogen  and  White's  (Ref. 

17)  in  air  are  shown  in  Fig.  6l  (a  reproduction  of  Fig.  39) • It  is  seen  that 
all  the  data  points  fall  inside  their  predicted  region.  Therefore,  our  analysis 
of  the  shock-wave  diffraction  in  nonstationary  diatomic  gases  is  substantiated, 
and  the  seven  different  types  of  diffractions  listed  in  Table  1 are  all  valid. 

It  is  worthwhile  to  refer  the  reader  to  White's  report  (Ref.  17)  where 
excellent  examples  of  shock-wave  diffractions  with  degenerated  reflected  shock 
waves  are  given  (Figs.  34,  38  and  42  of  Ref.  17).  Note  that  White  was  able  to 
obtain  very  weak  incident  shock  waves  (Mg  = 1.010,  1.047  and  1.022  for  these 
figures)  for  which  the  reflected  shock  wave  vanishes  as  it  degenerates  into  a 
Mach  wave.  For  further  details  see  Section  8. 

7.4.2  Monatomic  Gas 

Our  data  from  Table  6 in  argon  are  all  shown  in  Fig.  62  (a  reproduction 
of  Fig.  40) . All  experimental  points  (except  the  one  at  Ms  = 4.44,  Qy  = 4o* 
which  have  already  been  discussed)  lie  inside  these  predicted  regions.  Conse- 
quently, one  can  conclude  that  the  present  analysis  for  nonstationary  oblique 
shock-wave  diffraction  in  a monatomic  gas  is  substantiated.  Unfortunately,  out 
of  the  six  different  diffractions  predicted  by  the  present  analysis  (Table  2) 
only  five  were  observed  experimentally.  The  remainder,  a CMR  with  an  attached  shock 
wave  at  the  wedge  corner  (region  5>  Figs.  40  and  62),  was  not  observed  since  its 
domain  starts  at  Ms  =9.4  and  lies  beyond  the  practical  possibilities  of  incident 
shock  waves  (Ms  » 8)  without  damaging  the  interferometric  quality  windows  of  the 
test  section.  However,  in  light  of  Fig.  62  and  the  verification  of  5 regions  out 
of  the  6 listed  in  Table  2,  one  can  conclude  that  the  present  analysis  for  the 
diffraction  of  oblique  shock  waves  in  nonstationary  monatomic  flows  is  substantiated 
and  that  six  different  types  of  shock-wave  diffraction  exist  in  the  range  1 < Ms  < 

10. 
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7 .5  Summary  Discussions 

The  present  experiments  and  analyses  have  verified  the  regions  and 
boundaries  of  nonstationary  reflection  of  oblique  shock  waves  in  the  ( Ms,  0w) 
and  (Ms,  9W) -planes  as  well  as  those  of  ronstati onar y diffraction  of  oblique 
shock  waves  in  the  (Ms,  ewJ-plane  for  diatomic  and  monatomic  gases. 

It  is  worthwhile  to  refer  again  to  Fig.  3 which  shows  all  the  informa- 
tion concerning  these  domains  and  transition  boundaries  of  nons tationary  oblique 
shock- wave  reflection  in  the  (Ms,  q^) -plane  that  was  known  prior  to  the  present 
research.  A detailed  description  of  Fig.  3 is  given  in  Section  1,  and  hence 
only  a brief  comparison  follows.  While  Fig  59  covers  a wide  range  of  wedge 
angles  0 < 9W  < 90°  and.  shock  Mach  numbers  1 < Mg  < 10,  Fig  3 is  limited  to 
20°  < 9W  < 60°  (the  region  60°  < 0W  < 90°  corresponding  to  RR  was  probably  not 
drawn  since  it  does  not  contain  any  significant  information)  and  1 < Mg  < 9- 
However,  although  the  upper  limit  in  Fig.  3 is  Ms  = 9>  only  line  1 extends  all 
over  the  entire  range.  Lines  2,  3,  4 and  5 terminate  at  Ms  = 8.00,  4.75,  3*50, 
and  2.75,  respectively.  In  Fig.  59  all  lines  extend  over  the  entire  range,  up 
to  Mg  = 10.  However,  although  Fig  59  was  limited  to  Mg  = 10  (due  to  the  fact 
that  it  was  decided  to  make  experiments  only  up  to  Ms  = 8,  in  order  not  to 
damage  the  high  quality  test  section  windows)  the  present  analyses  apply  for 
Ms  > 10  as  well.  Unlike  the  case  of  Fig.  3 the  analytically- obtained  boundary 
lines  are  for  a perfect  gas  while  those  accounting  for  real-gas  effects  are  all 
experimental,  all  the  boundary  lines  in  Fig.  59  were  obtained  analytically, 
either  by  solving  Eqs.  2.1  to  2.9  for  an  RR  or  Eqs . 2.17  to  2.30  for  the  triple 
point  with  and  without  real-gas  effects.  It  is  probably  the  first  time  that  an 
analytical  solution  was  obtained  for  Eqs.  2.17  to  2.30  with  real-gas  effects. 

The  contributions  of  the  present  study  to  knowledge  concerning  the 
reflection  of  oblique  shock  waves  in  the  (Mg,  9w)  or  (Ms,  0o) -planes  are 
summarised  in  Tables  7 and  8.  While  RR  and  SMR  were  discovered  one  hundred 
years  ago  (1878)  by  Ernst  Mach,  it  is  only  about  35  years  since  GMR  was  first 
noticed  by  Smith  (1945)  and  DMR  discovered  by  White  (1951)  • The  formation 
criterion  of  RR  in  both  steady  and  nonstationary  flows  was  established  by 
Neumann  (1943),  as  well  as  its  termination  criterion  in  nonstationary  flows 
(1943)  More  recently  Henderson  and  Lossi  (1975)  suggested  a new  criterion  for 
the  termination  of  RR  that  was  found  (Ref.  65)  to  be  correct  only  in  steady 
flows.  However,  Hornung  and  Kychakoff  (1975)  established  a more  general 
criterion  for  the  termination  of  RR  resulting  in  the  sonic- criterion  (slightly 
below  Neumann's  criterion)  for  nonstationary  flows,  and  Henderson  and  Lossi 's 
criterion  for  steady  flows.  Henderson  and  Lossi  (1975)  were  the  first  to 
suggest  that  the  termination  of  SMR  is  followed  by  a CMR,  and  later  when  Mg 
increases  a DMR  is  formed.  However,  the  exact  criterion  for  the  termination 
of  CMR,  and  hence  the  formation  of  DMR,  was  established  only  in  the  present 
analytical  study  and  verified  by  experiments. 

Table  8 summarises  by  whom  and  when  the  boundary  lines  between  domains 
of  different  types  of  oblique  shock-wave  reflection  were  calculated.  For  a 
perfect  monatomic  gas  the  boundary  line  between  RR  and  SMR,  CMR  or  DMR  was 
calculated  by  Law  (1970)  and  all  the  other  boundary  lines  for  both  perfect  and 
imperfect  gases  were  calculated  in  the  present  study.  In  the  case  of  a diatomic 
gas,  the  perfect  gas  boundary  lines  RR/SMR,  CMR  or  DMR,  SMR/ CMR  and  CMR /DMR 
were  calculated  by  Neumann  (1943) , Bazhenova,  et  al  (1976)  and  Ben-Dor  (present), 
respectively.  (Note  that  since  Neumann  worked  in  low  Mach  numbers,  in  his  case 
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the  termination  of  RR  was  followed  by  an  SMR,  however,  his  criterion  of  the 
termination  of  RR  is  good  also  for  higher  values  of  Ms  where  RR  is  followed  by 
a CMR  or  DMR.)  All  the  boundaries  corresponding  to  an  imperfect  gas,  i.e., 
RR/SMR,  CMR  or  DMR,  SMR/CMR  and  CMR/DMB  were  calculated  in  the  present  study. 

8.  DENSITY-FIELD  INVESTIGATION 

The  experimental  facility  and  instrumentation  discussed  in  ..Chapter  6 
were  used  to  obtain  interferograms  of  nonstationary  oblique  shock- wave  reflec- 
tions. Lines  of  constant  density  (isopycnics)  were  obtained  from  them  using 
a method  developed  by  Whitten  (Ref.  75) • A detailed  description  of  this  very 
precise  data- evaluation  technique  is  given  in  Refs.  71  and  75-  Therefore,  only 
a brief  description  will  be  given  for  obtaining  the  isopycnics  (see  Appendix  E 
for  details) . 

Very  good  agreement  with  sophisticated  numerical  analyses  were  obtained 
previously,  using  the  identical  equipment  and  evaluation  techniques,  in  studies 
such  as  ionizing  shock-wave  structure  (Refs.  68  and  69)>  flat-plate  and  sidewall 
boundary  layers  (Refs.  70  and  75)  where  boundary  layer  refraction  errors  are 
more  significant  than  in  the  present  study.  Consequently,  one  can  conclude 
that  the  refraction  errors  from  sidewall  boundary  layers  had  a negligible 
effect . 

'The  errors  associated  with  the  present  experimental  and  evaluation 
techniques  are  discussed  in  detail  in  Ref.  71.  However,  since  the  number  of 
fringes  inside  the  shock-diffraction  region  was  small,  a more  conservative 
estimation  was  made  of  the  absolute  error  in  measuring  fringe  shifts.  The 
error  is  taken  as  E(s)  =0.1  rather  than  the  already  conservative  estimation 
used  by  Whitten  and  Ben-Dor  (Ref.  71),  i.e.,  0.05.  Consequently,  the  absolute 
errors  associated  with  the  densities  as  measured  from  the  interferograms  are: 
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= 4.292  X 

io-6 

for  argon 

(8.1) 

cur 
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for  nitrogen 
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Since  the  values  of  the  isopycnics  were  ail  normalized  wrt  the  flow  density 
ahead  of  the  incident  shock  wave  pQ,  the  relative  error  Ap/p0  is  used  through- 
out the  present  report. 

During  the  last  few  years  various  investigators  (Refs.  6l,  63  and  64) 
tried  to  develop  coaputer  codes  for  predicting  the  density  field  associated 
with  RR  and  SMR.  No  such  data  exist  for  CMR  and  DMR.  A detailed  comparison 
between  their  numerically  predicted  density  fields  and  the  present  results 
is  given  in  Ref.  73  and  briefly  .discussed  in  Section  8.2.  The  conparison 
reveals  that  all  numerical  methods  provide  reasonable  prediction  of  the  wave 
systems  and  their  shock  shapes  for  the  cases  that  were  compared.  However, 
the>  predict  rather  poor  values  and  locations  of  the  more  sensitive  indicators 
of  the  flow  isopycnics.  Consequently  the  various  codes  (Refs.  6l,  63  and  64) 
require  a reassessment  and  perhaps  a new  approach  in  the  light  of  the  disagree- 
ment with  the  present  detailed  and  very  accurate  interferometric  data.  Undoubtedly, 
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numerical  codes  will  evolve  in  the  future  that  will  reliably  predict  such 
complex  flow  in  inqperfect  gases.  No  doubt  the  present  data  will  provide  a 
solid  base  for  future  comparison. 

8 .1  Nitrogen 

The  seven  shock  wave  diffraction  domains  corresponding  to  regions  1 
to  7 of  Fig.  39  are  shown  in  the  interferograms  63a  to  63g,  respectively.  The 
density  distribution  (p/pD)  in  the  flow  fields,  in  terms  of  isopycnics  (n) 
associated  with  each  diffraction  process,  are  shown  in  Figs.  64a  to  64g.  The 
density  profiles  along  the  wedge  and  the  shock- tube  wall  appear  in  Figs.  65a 
to  65g.  In  the  following,  a general  description  is  given  of  each  diffraction 
as  well  as  their  similarities  and  differences. 


Figures  63a,  63b,  63c  and  63f  which  correspond  to  flow  deflections 
through  detached  shock  waves  (regions  1,  2,  4 and  6 of  Fig.  39) > shock  wave 
bifurcation  and  boundary  layer  roll-up  are  clearly  seen.  In  order  to  have  a 
closer  look  at  the  shock  wave  bifurcation  four  experiments  with  Ms  = 4.64, 

4.59>  4.60  and  4.72  and  By  = 40°  were  made.  The  process  was  recorded  26,  126, 
226  and  326  p.sec  after  the  incident  shock  wave  passed  the  centre  of  the  test 
section.  The  four  corresponding  interferograms  are  shown  in  Figs.  66a  to  66d, 
respectively.  Figure  66a  reveals,  as  expected,  a DMR  with  a detached  bifurcated 
shock  wave.  In  Figs.  66b  to  66d  the  incident  shock  wave  is  already  outside  of 
the  field  of  view;  however,  the  bifurcation  process  is  seen  to  be  growing  with 
time.  Figures  66b  and  66c,  for  example,  show  very  clearly  a bifurcated  shock 
and  a slipstream,  i.e.,  a three- shock  confluence  or  a triple  point.  Note  the 
clear  compression  at  the  wedge  corner  and  expansion  wave  as  the  wedge  becomes 
flat  giving  rise  to  a corner-expansion  flow.  In  Fig.  66d  the  reflected  shock 
wave  R is  seen  reflecting  from  the  upper  wall  of  the  shock  tube.  Unfortunately, 
the  upper  wall  is  not  in  the  field  of  view;  however,  it  looks  like  the  reflected 
wave  R reflected  through  a single-Mach  reflection  and  that  a Mach  stem  joins  it 
to  the  upper  wall . 


It  is  worthwhile  mentioning  that  White  in  his  pioneering  work  (Ref. 
17)  was  able  to  produce  some  very  weak  incident  shock  waves . He  evaluated  hi s 
interferograms  on  shock-wave  diffractions  at  very  low  Ms-  Figures  67a  and  67b 
are  reproductions  of  two  of  his  original  figures  (Figs.  34  and  46  of  Ref.  17). 
The  initial  conditions  for  these  two  pictures  are  Ms  = 1.010,  By  = 5-7°  (SMR) 
and  Ms  = 1.018,  0w  = 30°  (RR),  respectively.  For  these  low  values  of  Ms  the 
reflected  shock  wave  degenerates  to  a Mach  wave  and  the  flow  turns  over  the 
wedge  subsonically.  Note  that  for  these  two  examples  from  White  (Ref.  17)  the 
induced  flow  Mach  number  is  0.017  and  0.029,  respectively. 


When  the  flow  deflection  over  the  corner  is  achieved  through  an 
attached  or  detached  shock  wave  there  is  a sharp  density  jump  at  t point 
(Fig.  65) . In  the  case  of  an  attached  shock  wave  at  the  corner  (Figs.  63c, 

63e  and  63g)  the  highest  density  along  the  wall  is  measured  immediately 
behind  the  attached  shock  wave,  i.e.,  at  the  corner  (point  b,  Figs.  65c,  65e 
and  65g)  while  in  the  case  of  a detached  shock  wave  (Figs.  63b,  63d,  and  63g) 
the  highest  density  may  or  may  not  be  (Figs.  65a  and  65d)  behind  this  shock  wave 
or  at  the  corner.  It  can  occur  near  the  slipsteadm  (Fig.  65b)  or  near  the 
corner  (Fig.  65d) . 

The  existence  of  a compression  wave  at  the  kink  K of  a CMR  can  be 
cleanly  seen  in  Figs.  65d  and  64e  where  the  isopycnics  converge.  The  corre- 
sponding compression  strengths  (density  ratios)  are  approximately  6.880/6.447 
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and  8.778/8.127  or  1.067  and  1.080,  respectively.  For  the  compressions  of 
I.O67  and  1.080  an  increasingly  clearer  kink  can  be  seen  (Figs.  63d  and  63e) . 

The  equivalent  perfect  shock  wave  Mach  numbers  that  would  give  the  same 
compressions  are  1.040  and  1.048,  respectively.  The  calculated  flow  Mach 
numbers  behind  the  reflected  wave  ^ in  the  vicinity  of  the  triple  point  T 
are  1.251  and  1*319j  respectively.  Once  a DMR  is  formed  (Mg  > 1)  the  iso- 
pycnics  do  not  converge  any  more  (Figs . 64f  and  64g)  and  the  compression  wave 
is  replaced  by  a shock  wave. 

In  the  cases  of  SMR  (Figs.  64b  and  64c)  the  convergence  of  the 
isopycnics  corresponds  to  a weak  expansion  wave  rather  than  a compression  wave 
(follow  the  isopycnic  numbers) . The  strengths  of  those  expansion  waves  are: 
3-747/3-868  and  5-53/6.11  or  O.969  and  O.905,  respectively.  The  case  with  an 
attached  shock  wave  has  the  strongest  expansion.  This  is  also  true  for  CMR 
(Figs.  64c  and  64e) . Figures  65a  to  65g  show  that  the  density  along  the  wedge 
surface  always  increases  as  one  moves  from  the  Mach  stem  towards  the  point  where 
the  slipstream  disappears  into  the  boundary  layer.  Consequently  the  flow  behind 
the  Mach  stem  is  being  further  compressed. 

Although  the  density  flow  field  associated  with  the  various  shock 
wave  diffractions  differ  greatly,  nevertheless  they  do  have  some  similarities. 

For  example,  in  the  case  of  an  attached  shock  wave,  Figs.  64c  - SMR  and  64e  - CMR, 
the  isopycnics  tend  to  run  perpendicular  to  the  reflected  shock  wave.  In  the  case 
of  DMR  (Figs.  64 f and  64g)  a "corridor"  is  formed  for  the  second  slipstream.  In 
Fig.  63g  the  second  slipstream  is  not  visible,  probably  because  the  change  in  the 
density  is  not  sufficiently  large  to  establish  a noticeable  fringe  shift.  However, 
the  "corridor"  is  clearly  seen  in  Fig.  64g  and  consequently  a thin  dashed  line 
has  been  drawn  in  Fig.  65g  to  indicate  the  possible  location  of  the  second 
slipstream. 

The  density  at  any  point  (x,y)  of  Figs.  64a  to  64g  can  be  calculated 
either  by  interpolating  between  or  by  extrapolating  beyond  the  isopycnic  in 
the  vicinity  of  that  point.  However,  since  the  density  difference  between  the 
isopycnics  is  quite  small,  any  region  between  them  can  be  assumed  to  have  a 
uniform  average  value.  For  a region  where  the  density  change  was  not  suffi- 
ciently large  to  plot  isopycnics  it  can  be  assumed  as  uniform  with  the  indicated 
density  number.  For  example,  region  n = 6 bounded  by  R,  Rq  and  S in  Fig.  64g 
is  assumed  with  p = 8.87p0-  ™he  relative  error  hp/po  given  in  each  figure 
(Figs.  64a  to  64g)  was  obtained  from  Eq.  8.2  simply  by  dividing  with  the  corre- 
sponding value  of  pQ.  The  error  is  fixed  for  a particular  experiment.  It  can 
be  as  high  as  37 ■ 6%  for  PQ  « 5 torr  (Fig.  64g)  and  as  low  as  5-2$  for  PQ  « 37 
torr  (Fig.  64b) . Note  that  the  position  of  a point  on  any  isopycnic  was  drawn 
within  an  accuracy  of  ± 1 mm. 

8.2  Comparison  With  Some  Numerical  Data  and  Discussions 

The  different  conpared  cases  and  their  initial  conditions  are  shown 
in  Table  9.  Although  we  have  tried  to  perform  our  new  experiments  using  exactly 
the  same  initial  conditions  as  those  chosen  by  Schneyer  and  Kutler  and  Shankar 
for  their  numerical  analyses,  it  was  experimentally  convenient  to  use  nitrogen 
at  different  initial  pressures.  However,  owing  to  the  fact  that  nitrogen  and 
oxygen  can  be  treated  as  a perfect  gas  at  a shock  Mach  number  Mg  = 2.00,  the 
change  in  the  initial  pressure  was  not  significant.  Furthermore,  since  the 
value  of  the  flow  isopycnics  were  normalized  by  the  initial  density  ahead  of 
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the  incident  shock  wave  any  correctly- computed  and  actual  isopycnic  shapes  and 
values  must  be  the  same.  However,  for  strong  shock  waves  real  gas  effects  are 
important  and  must  be  taken  into  account,  for  HR  at  Ms  = 4.58  and  PQ  = 15-31 
torr . 


It  can  be  seen  from  Table  9 that  three  different  cases  were  compared. 
Case  1 results  in  an  RR  with  a weak  incident  shock  wave  (Ms  = 2.0)  and  it  was 
numerically  solved  by  Schneyer  and  Kutler  and  Shankar.  Case  2 is  again  an  KR 
but  with  a stronger  incident  shock  wave  (Ms  = 4.71)  and  it  was  analysed  by 
Kutler  and  Shankar.  Case  3 , an  SMR,  was  solved  by  Schneyer  using  two  different 
computer  codes,  the  two-dimensional  Euierian  code  THOR,  a revised  version  of 
the  HELP  (Ref.  78)  code,  and  the  two-dimensional  Lagrangian  code  CRAM,  based  on 
Wilkins'  formulation  (Ref.  79) • In  Kutler  and  Shankar  the  two-dimensional  time- 
dependent  Euler  equations  were  solved.  The  hyperbolic  partial-differential 
equations  were  transformed  to  introduce  self- similarity  and  the  distance  between 
the  corner  and  the  incident  shock  wave  was  used  for  normalization.  The  self- 
similar  transformation  reduces  these  equations  from  an  unsteady  to  a quasi- 
steady set  of  mixed  elliptic-hyperbolic  equations.  Then  the  equations  were 
made  totally  hyperbolic  by  reintroducing  a time-like  term.  The  final  set  of 
equations  were  written  in  a "strong  conservation  law  form"  and  solved  using 
MacCormack's  (Ref.  80)  second-order  finite-difference  algorithm. 

Comparison  With  Case  1 

The  shape  of  the  isopycnics  obtained  numerically  by  Schneyer,  Kutler 
and  Shankar  and  in  the  present  experiments  are  shown  in  Figs.  68a,  68b  and  68c, 
respectively.  The  contour  number  and  the  corresponding  density  ratio  are  tabu- 
lated in  Figs.  68a  and  68c.  Unfortunately  we  were  unable  to  obtain  the  contour 
numbers  from  Kutler  and  Shankar  (Ref.  63)  for  Fig.  68b.  It  is  worth  noting  that 
the  measured  density  ratios  immediately  behind  the  incident  and  reflected  shock 
waves  always  agree  well  with  theory.  It  can  be  seen  immediately  that  the  actual 
shapes  of  the  isopycnics  obtained  experimentally  appear  similar  to  those  pre- 
dicted by  Kutler  and  Shankar.  However,  the  results  from  Schneyer  do  not 
represent  the  physical  flow.  One  could  discard  the  isopycnics  shown  by  Schneyer 
using  the  following  argument . Schneyer 's  isopycnics  (even  if  their  lines  were 
taken  to  represent  the  shock  as  a result  of  artificial  viscosity  smearing) 
have  the  same  value  over  the  entire  length  of  the  reflected  shock  R.  This  means 
that  the  density  jimp  across  R is  constant  everywhere.  However,  since  the  angle 
of  incidence  between  the  flow  entering  R (in  a frame  of  reference  attached  to 
the  reflection  point  P)  decreases  as  R moves  a way  from  P,  the  strength  of  R 
should  increase  to  maintain  a constant  density  jump.  This  contradicts  both 
theory  and  experiment  where  the  shock-wave  strength  along  R decreases  as  the 
distance  from  the  reflection  point  increases,  as  shown  in  Fig.  68c.  Conse- 
quently, the  density  ratios  in  the  vicinity  of  3,  in  Fig.  68a  are  much  larger 
compared  with  the  actual  result  (Fig.  68c) . However,  just  the  range  of  density 
given  by  Schneyer  (Fig.  68a)  approximates  the  measured  range  (Fig.  68c) . Un- 
fortunately only  a qualitative  comparison  can  be  made  with  the  results  of  Kutler 
and  Shankar  since  their  values  corresponding  to  the  various  isopycnics  were  not 
given  on  the  figures  in  Ref.  63,  and  were  unavailable  from  private  communications. 

A comparison  between  the  actual  shock-wave  configuration  and  the  shapes 
predicted  by  Schneyer  and  Kutler  and  Shankar  is  shown  in  Fig.  68d.  Since 
several  isopycnics  of  Schneyer  represent  shock  waves  as  a result  of  artificial 
viscosity,  the  two  extreme  contours  were  reproduced  in  Fig.  63d.  The  distance 


between  the  incident  shock  wave  and  the  corner  is  normalised  for  all  shapes. 

It  can  be  clearly  seen  from  Fig.  68d  that  the  predicted  shapes  do  not  differ 
too  much  from  the  actual  shock-wave  configuration.  It  is  worth  mentioning 
that  Kutler  and  Shankar  predict  their  numerically-obtained  configuration  to 
be  slightly  larger  than  the  actual  one.  However,  their  explanation  that  this 
is  due  to  viscous  or  real-gas  effects  is  reasonable  for  the  latter.  It  can 
be  shown  analytically  that  vibrational  excitation  will  reduce  the  angle 
between  the  reflected  shock  wave  and  the  wedge  and  hence  will  result  in  a 
smaller  configuration.  The  angle  between  R and  the  wedge  surface  at  the 
reflection  point  is  16.34°  for  a perfect  gas  and  l6.04#  when  real-gas  effects 
are  considered.  Note  that  although  Schneyer  and  Kutler  and  Shankar  used  two 
different  computational  methods  which  disagree  in  the  prediction  of  the  entire 
density  field,  they  nevertheless  agree  in  the  shock  shapes  and  systems.  One 
can  only  conclude  that  the  isopycnics  are  much  more  sensitive  indicators  of 
the  accuracy  of  the  physical  flow  modelled  by  a specific  numerical  technique. 

Comparison  of  Case  2 

The  general  shapes  of  the  isopycnics  predicted  by  Kutler  and  Shankar 
(Fig.  69a)  do  not  agree  with  the  actual  ones  (Fig.  64a).  Their  predicted  shapes 
for  this  case  of  RR  are  almost  the  same  as  those  discussed  previously  (Fig.  68b) 
while  our  results  are  very  different  (Fig.  68c).  The  disagreement  between  the 
actual  and  numerical  isopycnics  may  arise  from  the  fact  that  Kutler  and  Shankar 
assume  a perfect  gas.  However,  for  Ms  = 4.71,  P0  = 15  torr,  T0  = 298  K,  real- 
gas  effects  cannot  be  neglected  as  the  vibrational  contribution  is  significant 
(e.g.,  pg/po  = 14.53  for  a perfect  gas  and  17-75  for  an  imperfect  gas) . Although 
this  will  change  the  absolute  numerical  values  of  the  isopycnics  it  might  also 
affect  their  shapes. 

Although  Kutler  and  Shankar  did  not  report  the  values  of  the  various 
isopycnics  on  their  reproduced  Fig.  69a,  we  deduced  representative  values  from 
their  density  distribution  along  the  wall  shown  in  Fig.  69b.  The  range  of  their 
densities  is  lower  than  the  present  owing  to  their  assumption  of  a perfect  gas. 
Note  that  the  strength  of  the  reflected  shock  wave  R decreases  from  the  reflec- 
tion point  to  the  shock-tube  wall  in  agreement  with  experiment  (Fig.  64a) . The 
numerically  calculated  density  distribution  along  the  wall  surface  predicted 
by  Kutler  and  Shankar  as  well  as  the  measured  values  close  to  the  wall  above 
the  thin  boundary  layer  are  shown  in  Figs.  69b  and  65a,  respectively.  The  fact 
that  the  actual  density  values  are  higher  than  those  predicted  by  Kutler  and 
Shankar  again  arises  from  their  perfect-gas  assumption.  For  their  case  of 
Ms  = 4.71  the  perfect-gas  value  p2/po  = l4-53j  while  for  the  present  case  the 
measured  value  of  pg/po  = l6.77>  and  it  lies  between  the  perfect  and  imperfect 
equilibrium  values  of  l4.48  and  17.65,  respectively. 

The  vortical  singularity  predicted  by  Kutler  and  Shankar  (V,  Fig.  69a) 
or  point  d (Fig.  69b)  cannot  be  seen  in  the  interferogram.  Therefore  the  curve 
from  b to  f (Fig.  65a),  unlike  the  curve  in  Fig.  69b,  has  no  discontinuity. 

If  one  existed  it  would  be  smeared  out  by  the  boundary  layer . Both  curves 
consist  of  a sharp  density  jump  P20  (f  - f)  at  the  reflection  point  followed 
by  a uniform-density  region  (which  is  longer  in  Fig.  69b)  terminated  by  an 
expansion.  At  the  detached  shock-wave  position  b another  sharp  rise  takes 
place  (the  rise  is  greater  in  Fig.  69b) . The  fringe  pattern  due  to  the 
shock-boundary-layer  interaction  at  the  corner  was  too  complex  to  analyse. 
Therefore  the  density  ratio  was  extrapolated  to  the  location  of  the  detached 
shock  wave  near  the  corner. 
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The  actual  shock  wave  configuration  and  the  one  predicted  by  Kutler 
and  Shankar  are  shown  in  Fig.  69c.  This  figure  agrees  with  their  statement 
that  "the  experimental  shock  location  would  fall  inside  the  numerical  solution". 

As  mentioned  earlier,  the  reason  for  this  lies  in  real-gas  effects  rather  than 
viscous  effects  suggested  by  Kutler  and  Shankar.  The  angle  between  the  reflected 
shock  wave  R and  the  wedge  surface  for  a perfect  gas  is  l4.80°  and  12.45°  for 
a real  gas.  Consequently  the  actual  shape  is  smaller  than  the  one  obtained 
from  a perfect-gas  model. 

Comparison  With  Case  3 

The  shapes  of  the  isopycnics  obtained  numerically  by  Schneyer  (Ref. 

6l)  using  the  Lagrangian  and  Eulerian  computer  codes  as  well  as  those  obtained 
experimentally  are  shown  in  Figs.  70a,  70b  and  70c,  respectively.  It  is  evident 
that  the  Lagrangian  and  Eulerian  results  differ  quite  considerably  from  each 
other  and  from  the  actual  isopycnics.  The  numerical  configurations  suffer  from 
the  same  artificially  viscous  spreading  of  the  incident  and  reflected  shock 
waves.  The  isopycnics  maintaining  the  same  values  along  R imply  a physically 
unrealistic  reflected  shock  wave  of  increasing  strength  as  it  moves  away  from 
the  triple  point  T.  Spurious  isopycnics  (n  = 4 and  5 , Fig.  70a)  appear  from 
the  Lagrangian  code  in  the  middle  of  the  reflected  wave  R.  There  are  no 
density  contours  generated  in  the  important  region  between  the  Mach  stem  M and 
the  slipstream  S.  The  Eulerian  contours  fail  to  predict  the  existence  of  a 
slipstream  altogether  (Fig.  70b).  If  a line  is  drawn  at  the  estimated  location 
of  the  slipstream,  it  intersects  the  isopycnics  implying  that  the  densities  on 
both  sides  of  the  slipstream  are  equal  in  violation  of  the  physical  condition 
that  the  slipstream  divides  two  thermodynamic  regions  of  different  densities, 
even  if  not  large.  Schneyer  (Ref.  6l)  attributes  the  appearance  of  the  spurious 
expansion  and  shock  waves  in  the  middle  ofR  in  the  Lagrangian  results  to  an 
inexact  choice  of  the  initial  velocity  profile.  He  explains  the  disappearance 
of  the  slipstream  (in  the  Eulerian  result)  as  being  "washed  out"  by  the  "effective 
[artificial]  viscosity". 

The  actual  shape  of  the  isopycnics  (Fig.  70c)  show  very  clearly  that 
the  densities  on  both  sides  of  the  slipstream  are  different.  The  densities 
behind  the  reflected  shock  wave  are  higher  than  those  behind  the  Mach  stem, 
as  expected  from  gasdynamic  sondierations . The  approximate  density  ratio  across 
the  slipstream  is  1.12  (3-67/3-27  near  the  triple  point)  shows  that  it  is  indeed 
a weak  discontinuity.  The  analytical  density  ratio  across  the  slipstream  in 
the  vicinity  of  the  triple  point  is  1.10.  A qualitative  comparison  between 
the  shapes  of  the  actual  isopycnics  and  those  of  Schneyer  indeed  show  poor 
agreement.  However,  as  in  the  previous  case  only  the  range  of  density  values 
approximate  those  obtained  experimentally. 

A comparison  between  the  normalized  predicted  shock  shapes  and  the 
actual  wave  configuration  is  shown  in  Fig.  70d.  The  Lagrangian  wave  system  is 
somewhat  larger  than  the  Eulerian.  The  actual  shock  wave  system  lies  close  to 
or  inside  the  Eulerian  shape  boundaries.  The  slipstream  predicted  by  the 
Lagrangian  code  agrees  reasonably  well  with  experiment.  It  is  worth  repeating 
that  also  in  this  case  the  Lagrangian  and  Eulerian  codes  predict  the  actual 
wave  system  quite  well.  However,  the  isopycnic  fields  and  the  varying  strength 
of  the  reflected  wave  are  poorly  represented. 

It  can  be  summarj  zed  that  the  present  experiments  show  that  all 
numerical  methods  provide  reasonable  predictions  of  the  wave  systems  and  their 
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shapes  for  the  two  analysed  cases  of  RR  and  SMR.  No  numerical  data  exist  for 
the  cases  of  CMR  and  DM*.  The  numerical  codes  predict  rather  poor  values  and 
locations  for  the  more  sensitive  indicators  of  the  flow  properties,  namely  the 
isopycnics.  Of  the  various  numerical  analyses  produced  so  far,  the  one  of 
Kutler  and  Shankar  (Ref.  63)  and  Shankar,  Kutler  and  Anderson  (Ref.  64)  are 
superior . Even  their  codes  require  a reassessment  and  perhaps  a new  approach 
in  the  light  of  the  disagreement  with  the  detailed  and  very  accurate  interfer- 
ometric data  presented  here . 


Undoubtedly,  numerical  codes  will  evolve  in  the  future  that  will 
reliably  predict  not  only  RR  and  SMR  but  also  CMR  and  DMR  in  real  gases.  The 
present  interferometric  data  of  all  these  cases  should  provide  a solid  base 
for  comparison . In  the  meantime,  those  laboratories  that  have  shock  tubes 
equipped  with  interferometers  will  benefit  from  experiments  in  nonstationary 
flows  in  order  to  check  their  numerical  analyses. 


9-  CONCLUSIONS 

The  criteria  for  the  formation  and  termination  of  RR,  SMR,  CMR  and 
DMR,  their  domains  and  transition  boundaries  in  both  nitrogen  and  argon,  in 
the  (Ms,  0w) -plane  for  nonstationary  flows  and  in  the  (Mo,  <fo) -plane  for 
steady  flows  have  been  established  analytically  and  verified  experimentally. 

The  equations  of  motion  for  RR  and  SMR  in  a perfect  and  an  inperfect 
gas  (nitrogen  in  dissociation  equilibrium  and  argon  in  ionization  equilibrium) 
have  been  solved,  and  the  significance  of  real-gas  effects  on  shifting  the 
boundary  lines  between  the  domains  of  different  reflections  was  shown  for  the 
first  time. 


Analytical  methods  for  predicting  the  values  of  x and  x'  were 
developed,  and  consequently  the  reflection  process  was  transformed  from  the 
(Mg,  ew) -plane  to  the  more  physical  (Ms,  ©„) -plane.  The  reflection  phenomenon 
in  the  (Ms,  0W) -plane  was  then  superimposed  with  the  flow  deflection  process 
in  the  (Ms,  0W) -plane  to  yield  the  overall  shock-wave-diffraction  phenomenon. 

It  is  shown  that  the  range  1 < Ms  < 10  seven  different  types  of  diffraction  are 
possible  in  nitrogen  and  six  in  argon  depending  on  the  incident  shock  wave 
Mach  number  Ms,  the  wedge  angle  q^,  and  the  initial  pressure  PQ  and  tempera- 
ture T0.  The  present  analyses  were  all  substantiated  by  over  120  interfer- 
ometric experiments  conducted  in  the  UTIAS  10  cm  x l8  cm  Hypervelocity  Shock 
Tube,  as  well  as  many  other  data  in  nitrogen,  oxygen  and  air  for  a diatomic 
gas  and  argon  and  helium  for  a monatomic  gas  obtained  from  other  sources.  This 
has  brought  new  order  and  understanding  of  the  various  results  from  different 
researchers . 


One  significant  problem  remains,  namely  a more  accurate  analytical 
method  for  calculating  the  triple-point  trajectory  angle  x than  the  one 
developed  earlier  by  Law  and  Glass  (Ref.  6l)  using  a graphical  technique  and 
presently  by  using  a more  accurate  analytical  method. 

The  density  fields  of  the  different  diffraction  processes  have  been 
deduced  from  the  corresponding  interferograms,  and  compared  with  available 
numerical  predictions . 


The  very  comprehensive  isopycnic  data  are  the  first  since  the  early 
pioneering  work  of  White  (Refs.  17  and  18)  who  first  noticed  the  four  types  of 
reflection.  The  results  provide  an  important  base  for  testing  available  and 
future  computational  codes  describing  such  complex  flows.  Although  the 
numerical  methods  can  satisfactorily  predict  the  gross  features  of  tne  wave 
system  and  shock  shapes  for  regular  and  single-Mach  reflections,  they  are  as 
yet  unsatisfactory  for  predicting  the  isopycnics  of  the  flow  (Ref.  73).  Ho 
computational  data  presently  exist  for  complex  and  double-Mach  reflections. 
Undoubtedly,  such  codes  will  evolve  in  the  near  future. 
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Table  1 


Diffraction  Regions  in  Nitrogen 
(Fig.  39) 


Table  2 

Diffraction  Regions  in  Argon 
(Fig.  40) 


Table  3 

Initial  Conditions  for  Obtainin, 


Various  Shock  Wave  Mach  Numbers  - Nitrogen 


M 

s 

P4l 

Driving  Gas 

2.0 

110 

co2 

3-7 

220 

He 

4.8 

690 

He 

6.2 

800 

He 

7.0 

1550 

H2 

7.8 

3320 

«2 

Table  4 

Initial  Conditions  for  Obtaini 


Various  Shock  Wave  Mach  Numbers  - Argon 


M 

s 

P4i 

Driving  Gas 

2.0 

60 

C°2 

3-0 

70 

He 

4.4 

350 

He 

5-2 

800 

He 

6.1 

520 

H2 

6.9 

1000 

H2 

7-9 

1900 

H2 

9w 

M 

s 

P 

0 

To 

X 

H 

Reflec- 

tion 

Exp. 

2 

1.95 

52.50 

297.2 

26.0 

SMR 

70 

2 

I.89 

53.50 

297.3 

26.5 

SMR 

71 

2 

1.85 

52.50 

297.4 

26.5 

SMR 

72 

2 

3.84 

15.19 

297.4 

23.5 

SMR 

73 

2 

4.15 

15.17 

297.3 

23.0 

SMR 

74 

5 

3-75 

15.25 

295.9 

20.5 

SMR 

37 

5 

4.71 

15.81 

295.2 

20.0 

SMR 

31 

5 

5.85 

15.18 

297.0 

19-5 

SMR 

33 

5 

6.01 

15 .19 

295.8 

18.5 

SMR 

34 

5 

6.86 

10.00 

295.6 

18.0 

SMR 

35 

5 

7.51 

5.17 

296.0 

17.5 

SMR 

36 

10 

2.01 

50.00 

295.8 

19.0 

SMR 

39 

10 

2.37 

35.44 

297.7 

18.5 

SMR 

88 

10 

2.61 

37.00 

297.8 

18.0 

SMR 

90 

10 

2.82 

30.34 

297.6 

18.0 

SMR 

89 

10 

3.62 

15.23 

295.4 

16.5 

SMR 

40 

10 

4.59 

15.16 

298.5 

16.2 

SMR 

5 

10 

4.72 

15.00 

295.0 

16.0 

SMR 

4l 

10 

5.92 

15.27 

295.0 

15.5 

SMR 

42 

10 

6.79 

10.21 

295.2 

15.0 

SMR 

43 

10 

7.58 

5.13 

294.8 

l4 .5 

SMR 

44 

20 

1.93 

51.00 

297.2 

12.5 

SMR 

50 

20 

3.74 

15.31 

297.4 

12.0 

CMR 

49 

20 

4.81 

15.29 

296.6 

11.5 

15-5 

CMR 

48 

20 

6.27 

15-33 

296.0 

11.2 

14.5 

CMR 

47 

20 

6.87 

10.12 

295.8 

11.0 

l4.o 

CMR 

46 

20 

7.71 

5.06 

296.O 

10.0 

11.5 

CMR 

45 

26.56 

2.01 

50.00 

296.6 

9.2 

SMR 

26 

26.56 

8.06 

5-10 

298.2 

9-0 

9-9 

DMR 

102 

30 

1-97 

51.00 

297.4 

8.5 

SMR 

51 

30 

3.68 

15.27 

297.3 

8.0 

10.0 

CMR 

52 

30 

4.68 

15.28 

297.4 

7.8 

9.5 

DMR 

53 

30 

5-93 

15.22 

297.4 

7.7 

10.0 

DMR 

54 

30 

6.96 

10.11 

297.4 

7.6 

9.8 

DMR 

55 

30 

7-97 

4.99 

297.4 

7.4 

9.0 

DMR 

56 

Continued 
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Initial  Conditions  for  the  Experiments  in  Nitrogen 


Table  6 


Initial  Conditions  for  the  Experiments  in  Argon 


e 

M 

P 

1 

W 

s 

o 

o 

2 

2.03 

50.00 

297.4 

28.5 

SMR 

84 

2 

3.02 

20.29 

297.8 

28.0 

SMR 

83 

2 

4.39 

15.00 

297.4 

27.5 

SMR 

75 

2 

5.19 

15.30 

297.2 

27.0 

SMR 

78 

2 

5-33 

5-04 

297.3 

76 

2 

5.42 

5.08 

297.3 

77 

2 

6.13 

15-33 

296.0 

SMR 

80 

2 

6.47 

15.32 

295.4 

26.0 

SMR 

81 

2 

7.77 

9.80 

297.6 

25.0 

SMR 

82 

10 

2.01 

50.00 

298.6 

21.0 

SMR 

85 

10 

2.96 

20.28 

299-0 

20.0 

SMR 

86 

10 

4.39 

15.32 

297.0 

19-5 

SMR 

87 

10 

5-22 

15.22 

298.4 

19.2 

SMR 

91 

10 

6 .06 

15.24 

299.O 

18.5 

SMR 

92 

10 

6.47 

15.27 

299.O 

18.5 

SMR 

93 

10 

7.88 

9.96 

298.6 

17.5 

SMR 

94 

20 

2.00 

50.00 

298.4 

15.0 

SMR 

101 

20 

2.82 

20.32 

299.O 

14 .5 

SMR 

100 

20 

4.4o 

15.26 

299-0 

l4 .0 

CMR 

99 

20 

5.20 

15.22 

299-0 

i4.o 

17.0 

CMR 

98 

20 

6.04 

15.27 

297.2 

13.7 

17.5 

CMR 

97 

20 

6.84 

15.22 

298.4 

l4 .0 

18.0 

CMR 

96 

20 

7.76 

9-84 

299.0 

i4.o 

17.5 

CMR 

95 

30 

2.03 

50.00 

299-6 

9-5 

SMR 

103 

30 

2.89 

20.24 

299.2 

9-5 

12.0 

CMR 

104 

30 

4.51 

15.25 

299-0 

10.0 

13.0 

CMR 

105 

30 

5.29 

15.21 

299-4 

10.0 

12.5 

CMR 

106 

30 

6.36 

15.27 

299-4 

10.2 

13.0 

CMR 

107 

30 

6.96 

15.00 

295.4 

10.0 

13.0 

CMR 

109 

30 

8.01 

9.80 

299-5 

10.0 

13.0 

CMR 

108 

4o 

2.05 

50.00 

297.8 

6.0 

11.0 

CMR 

116 

4o 

3-11 

20.34 

299.8 

5-5 

10.0 

CMR 

115 

4o 

4.44 

15.00 

299-1 

5-5 

10.0 

DMR 

114 

4o 

5.28 

15.29 

297.9 

5.5 

10.0 

DMR 

113 

4o 

6.12 

15.32 

297.6 

5-7 

10.0 

DMR 

112 

4o 

6.81 

15.23 

298.8 

5.7 

9.5 

DMR 

111 

4o 

7.53 

9-87 

297.0 

5.5 

9.5 

DMR 

110 

50 

2.04 

50.00 

298.2 

1.0 

CMR 

117 

50 

2.96 

20.50 

298.4 

1.5 

4.0 

DMR 

118 

50 

4.4o 

15.30 

299.2 

2.0 

4.0 

DMR 

119 

50 

5.27 

15.32 

298.2 

1.5 

4.0 

DMR 

120 

50 

6.27 

15.34 

299-4 

1.4 

4.0 

DMR 

121 

50 

7.03 

15.29 

299.4 

1.3 

3.5 

DMR 

122 

Continued. 
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Table  7 


Discovery  of  Four  Types  of  Oblique-Shock-Wave  Reflections 


Discovered  by 


Mach 
Ref.  1 
1878 


Mach 
Ref.  1 
1878 


Suggestion  of 

Criterion  for 

Formation 

Termination 

* 

Neumann 

Hornung  et  al 

Ref.  2 

Ref.  65 

19^3 

1976 

* 

Hornung  et  al 

White 

Ref.  65 

Ref.  17 

1976 

1951 

White 

Ben-Dor 

Ref.  17 

Present 

1951 

1978 

Ben-Dor 

Present 

*-** 

1978 

Smith 
Ref.  6 
19^5 


White 
Ref.  17 
1951 


* The  termination  criterion  of  RR  and  formation  criterion  of  SMR  for  non- 
stationary flows  were  established  by  Neumann  (1943)  and.  for  steady  flows 
by  Henderson  and  Lozzi  (1975)-  However,  Hornung  and  Kychakoff  established 
a more  general  criterion  which  holds  for  both  steady  and  nonstationary  flows 

**  Recall  that  CMR  and  DMR  do  not  exist  in  steady  flows. 

***  For  the  time  being  DMR  may  not  terminate.  However,  it  is  possible  that  for 
very  high  values  of  Mg,  the  flow  Mach  number  in  state  (5),  see  Fig.  19, 
will  become  so  high  that  a shock  wave  will  be  needed  to  prevent  it  from 
colliding  with  the  wall.  Under  these  circumstances  DMR  will  terminate, 
and  a triple-Mach  reflection  (TMR)  might  be  formed. 


Calculation  of  Boundary  Lines  Between  Reflection  Domains 


Monatomic  Gas  Diatomic  Gas 


Perfect 

Imperfect 

* 

/ SMR 

RR  / CMR 
/ DMR 

Law 

Ref.  48 
1970 

Ben-Dor 

Present 

1978 

Neumann 

Ref.  2 

19^3 

Ben-Dor 

Present 

1978 

DMR/CMR 

Ben-Dor 

Present 

1978 

Ben-Dor 

Present 

1978 

Ben-Dor 

Present 

1978 

Ben- Dor 
Present 

1978 

cmr/smr 

Ben-Dor 

Present 

1978 

Ben-Dor 

Present 

1978 

Bazhenova  et  al 
Ref.  55 

1976 

Ben-Dor 

Present 

1978 

* Recall  that  the  termination  of  RR  can  result  in  either  an  SMR,  CMR  or 
DMR  depending  on  the  value  of  Ms  (see  Figs.  39  and  40  for  details). 


Table  9 

The  Compared  Reflection  Cases  and  Their  Initial  Conditions 


Inperfect 

Gas 

Perfect  Gas 

Initial* 

Present 

Schneyer 

Kutler 

Type 

Conditions 

Experiment 

Eulerian 

Lagrangian 

& Shankar 

Gas 

N2 

7=1.40 

7=1.40 

0 0 
w 

63.43 

63.43 

63.41 

M 

2.01 

2.00 

2.00 

1 

s 

PQ  (torr) 

50.00 

760 

760 

s 

0 

T0  (K) 

298.6 

H 

PQ(g/cm  ) 

7. 57x10" 5 

10"3 

10"  3 

& 

« 

Gas 

N2 

°2 

1 

0 0 

60.00 

60 

$3 

Ph 

w 

M 

4.68 

4.71 

2 

s 

P ( torr) 

15-31 

15 

Tq  (K) 

298.1 

298.6 

Po(g/cm3) 

2. 31x10" 5 

2. 50x10" 5 

s 

Gas 

N2 

7=1.40 

7=1.40 

0 

H 

V 

26.56 

26.56 

26.56 

I 

M 

2.01 

2.00 

2.00 

gj 

3 

s 

w 

0 

PQ  (torr) 

50.00 

1 

a 

8 

H 

T0  (K) 

296.6 

PQ(g/cm3) 

7. 57x10" 5 

10"3 

10"  3 

U1 

* Initial  conditions  are  for  the  quiescent  gas  ahead  of  the  incident  shock  wave. 


NITROGEN 

(a)  REGULAR  REFLECTION  (RR) , 0w  = 

(b)  SINGLE-MACH  REFLECTION  (SMR) , 

(c)  COMPLEX-MACH  REFLECTION  (CMR) , 

(d)  DOUBLE -MACH  REFLECTION  (DMR) , 


Actual  Wedge  Angle  - 9\ 


Incident  Shock  Wave  Mach  Number  - Ms 

FIG.  3.  DOMAINS  AND  BOUNDARIES  OF  SHOCK  WAVE  REFLECTION  IN  (M  6 ) -PLANE 
(ENLARGE  REPRODUCTION  OF  FIG.  S,  REF.  62).  SOLID  LIN^S  W 
ARE  ANALYTICAL,  DASHED  LINES  ARE  EXPERIMENTAL. 


FIG.  15.  I AND  R-POLAR  COMBINATION  FOR  2.23  < M < 2.40,  y = 1.40 
M = 2.30,  8'  = 46.27°,  M * 1.59.  ° 


FIG.  17.  SCHEMATIC  DIAGRAM  OF  A COMPLEX-MACH  REFLECTION  (CMR) . 
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Trajectory  Angle 


FIG.  41.  EQUIVALENT  OBLIQUE  SHOCK-NAVE  REFLECTIONS  IN  STEADY  AND  NONSTATIONARY  FLOWS 

(a)  RR  IN  STEADY  FLOW 

(b)  RR  IN  NONSTATIONARY  FLOW 
(C)  SMR  IN  STEADY  FLOW 

(d)  SMR  IN  NONSTATIONARY  FLOW 


(b) 

FIG.  47.  THE  INTERRUPTION  OF  THE  FLOW  IN  STATE  (2)  BY  THE  GEOMETRICAL 
SET-UP.  (a)  Comer  expansion  wave,  (b)  Reflected  shock  wave. 


ATLANTIC  RESEARCH  U)-2 5 PRESSURE  TRANSDUCERS 
* KISTLLR  601-B  PRESSURE  TRANSDUCER 
P PREAMPLIFIER 
A AMPLIFIER 

CA  ICISTLER.  MODEL  504,  CHARGE  AMPLIFIER 


FIG.  51.  INSTRIM.NTATION  SET-UP  FOR  RECORDING  THE  SHOCK -WAVE  DIFFRACTION  PROCESS. 
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FIG.  52.  ERRORS  ASSOCIATED  WITH  THE  PRESENT  F.XPERIICNTS. 

(a)  Relative  error  in  shock-wave  Mach  niasber. 

(b)  Absolute  error  in  pressures. 
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26 
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Incident  Shock  Wave  Mach  Number  - Ms 

PIC.  S6.  VARIATION  OF  * WITH  M#  FOR  SMALL  \ AND  COMPARISON  WITH 
PRESENT  EXPERIMENTS . SOLID  LINES  - IMPERFECT  NITROGEN, 

Po  • 15  TORR,  T0  ■ 300  R.  DASHED  LINES  - PERFECT  DIATOMIC 
GAS  y • 1.40. 


FIG.  S7.  VARIATION  OF  x WITH  M,  FOR  A GIVEN  #„  AND  COMPARISON  WITH 
PRESENT  EXPERIMENTS.  SOLID  LINES  - IMPERFECT  ARGON, 

P0  • IS  TORR,  T0  - 300  K.  DASHED  LINES  - PERFECT  MONATOMIC 
GAS  y • S/3. 


Incident  Shock  Wave  Mach  Number  - Ms 

FIG.  58.  VARIATION  OF  x *ITH  M,  FOR  SMALL  #„  AND  COMPARISON  WITH  PRESENT 
EXPERIMENTS.  SOLID  LINES  - IMPERFECT  ARGON.  P0  - IS  TORR. 

T0  • 300  X.  DASHED  LINES  - PFRFECT  MONATOMIC  GAS  y • S/3. 


Actual  Wedge  Angle 


Incident  Shock  Wave  Mach  Number  - M3 

FIG.  S9.  EXPERIMENTAL  VERIFICATION  OF  ANALYSIS  OF  OBLIQUE  SHOCK-WAVE  RE- 
FLECTION IN  (M# ,6*) -PLANE  FOR  DIATOMIC  GAS.  LINES  ARL  FOP.  IMPER- 
FECT NITROGEN.  PQ  « 15  TORR,  T0  - 300  K.  EXPERIMENTAL  DATA:  AIR  - 
SMITH  A®  (REF.  6).  WHITE  *4®  (REF.  17),  KZHENOVA  El  ALBV 
(P.EF.  62)  ; OXYGEN  - LAW  AND  CLASS  CD  VAO( REF.  51)  ; NITROGEN  - 
BEZMENOVA  ET  AL  A (REF.  6.’)  AND  PRESENT  Ita# 


Incident  Shock  Wave  Mach  Number  - Ms 

FIG.  60.  EXPERIMENTAL  VERIFICATION  OF  ANALYSIS  OF  OBLIQUE  SHOCK-WAVE 

REFLECTION  IN  (Ms, ®w) -PLANE  FOR  MONATOMIC  GAS.  LINES  ARE  FOR 
IMERFECT  ARGON,  PQ  • 15  TORR,  Tp  - 300  K.  EXPERIMENTAL  DATA: 
HE LI l*  - LAW  ▼ (Ref.  48);  ARGON  - LAW  AND  GLASS  V (Ref.  51), 
BAZHENOVA  ET  AL  B (Ref.  62)  AND  PRESENT  |T4| 


FIG.  61.  EXPERIMENTAL  VERIFICATION  OF  ANALYSIS  OF  TEN  DOMAINS  AND  THFIP 
TRANSITION  BOUNDARIES  CF  XONSTATfONARY  OBtlQUF  SHOCK  M\T- 
Dtf  FPJVCTION.  E XPf.RI’IESTAI.  DATA:  AIR  - WHITE  0 44 
(Ref.  17);  NITROGIN  - PRESENT  IDVffAAliOl 
LINES  ARE  FOR  IMPERFECT  NlTPOGFN,  P0  « IS  TORR,  T0  • 300  K. 


Incident  Shock  Wave  Mach  Number  - Ms 

TIG.  62.  EXP!  R I Ml  NTAL  VERIFICATION  Of  MM  DOMAINS  AND  THI  1R  TRANSITION 
BOUNDARIES  OF  NONSTATIONARY  'BLIOUT  SHOCK-WAVI  DIFFRACTION. 
EXPERIMENTAL  DATA:  ARG0N-GV®A4  4 • 

PRLSEYT  RESUL1S.  UNIS  ARJ  ICR  IMITRIEH  ARGON  Pp  • IS  TORR 
Ta  • 300  K. 


298.2  k. 


(g)  M = 8.0b,  0 

2>  »N 


2b.5b°,  P =5.11  torr,  T = 
o o 


n p/Pg 


FIG.  64  ACTUAL  FLOW  ISOPYCNICS  CORRESPONDING  TO  THE  SEVEN  DIFFERENT  OBLIQUE 
SHOCK  WAVE  DIFFRACTIONS  SHOWN  IN  FIG. 63  (WITH  IDENTICAL  INITIAL  CON- 
DITIONS) . THE  IMPERFECT  AND  PERFECT  (IN  BRACKETS)  DENSITY  (p/PQ) 
VALUES  ARE  GIVEN  FOR  POINT  a,  b AND  c AT  THE  UCINITY  OF  THE  TRIPLE 


POINT  OR  THE  REFLECTION  POINT, 
(a)  a-S.32,  b-17.66  (4.88,  14.48) 


sity  ratio  distribution  along  the  wall-yO/yD( 


Distance  from  the  corner  along  the  wall  — r (mm) 

FIG.  65  (a)  INTERFEROMETRIC  DENSITY-RATIO  DISTRIBUTION  ALONG  THE  WALL  FOR 
DIFFERENT  OBLIQUE  SHOCK-WAVE  DIFFRACTIONS  SHOWN  IN  FIG.  63  (WITH 
IDENTICAL  INITIAL  CONDITIONS) 


Distance  from  the  corner  along  the  wall-r  (mm) 


FIG. 65  (c)  CONT. 
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Distance  from  the  corner  along  the  wall  -r  (mm) 

FIG. 65  (e)  CONT. 
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(c)  ACTUAL  FLOW  ISOPYCNICS,  NITROGEN.  M - 2.01.  0 • 2t.St V 


(d)  COWAR  ISON  OF  NORMALIZED  SHOCK-WAVE  SYSTEMS  OF  SCHNEYER  WITH  EXPERIMENT . 


APPENDIX  A 


(P,  9) -SHOCK  POLARS 


The  locus  of  all  possible  post-shock  states  that  result  from  a given 
flow  with  pressure  P,  temperature  T and  Mach  number  M,  when  it  passes  through 
any  oblique  shock  wave,  is  the  shock  polar-  (Fig.  A.l)  . 

Two  important  deflection  angles,  the  sonic  angle  es  and  the  angle 
of  maximum  deflection  (the  detachment  angle)  0m,  divide  the  shock  polar  (and 
the  deflection  process  as  well)  into  three  different  cases: 

(1)  0 < 0 < es;  the  deflection  is  through  a straight  and  attached  shock 
wave  (Fig.  5b)  behind  which  the  flow  is  supersonic. 

(2)  9S  < 0 < 9j£$  the  deflection  is  through  a curved  but  attached  shock 
wave  behind  which  the  flow  is  subsonic  (Fig.  5c). 

(3)  0 > 0m>  where  the  curved  shock  wave  is  detached  (Fig.  5d) . 

The  part  of  the  shock  polar  below  m is  called  the  weak- shock  solution  portion 
and  the  remainder  of  the  polar  is  for  strong  shocks.  Consequently,  for  every 
deflection  angle  0 the  pressure  ratio  is  double  valued.  However,  it  is  an 
experimental  fact  that  in  the  case  of  sharp  compressive  corners  the  weak  solu- 
tion occurs.  Note  that  in  the  presence  of  boundary  layer  or  blunt  bodies  the 
strong  solution  can  occur. 

Each  point  on  the  shock  polar  corresponds  to  a state  which  could 
have  been  obtained  by  passing  through  an  appropriate  oblique  shock  wave. 
Consequently,  there  are  two  other  points  of  interest  in  every  shock  polar, 
namely  the  limiting  cases  corresponding  to  a Mach  wave  (0o  = H = sin-1  l/u) 
and  a normal  shock  wave  (0O  = 90°) . The  former  represents  no  change  in  the 
flow  properties,  while  the  latter  signifies  the  greatest  changes  in  the  flow 
properties.  Note  that  the  flow  deflection  for  both  cases  is  aero  (0=0). 

It  is  well  known  (Ref.  27)  that  the  higher  the  flow  Mach  number 
Mo  the  greater  the  corresponding  maximum  deflection  angle  0m.  This  is  illus- 
trated in  Fig.  A. 2 where  the  deflection  angle  increases  from  43 .2°  to  49.0° 
as  Mo  increases  from  5 to  10.  Consequently,  a flow  with  Mo  = 10  would  negotiate 
a wedge  of  45®  through  an  attached  shock  wave,  while  the  deflection  of  a flow 
with  Mq  = 5 would  be  through  a detached  shock  wave.  A detailed  discussion  of 
the  influence  of  real-gas  effects  on  the  (P,  0)-polar  is  given  in  Refs.  48  and 
51.  A typical  illustration  is  shown  in  Fig.  A. 3.  It  is  seen  that  both  the 
pressure  ratio  P/PQ  across  the  oblique  shock  wave  and  the  angle  of  maximum 
deflection  0m  increase  (see  NSW  points  for  example)  as  a result  of  real-gas 
effects.  For  this  particular  illustration  of  Mq  = 10,  P0  = 15  torr  and 
To  = 300  K,  0m  = 42.7°  for  a perfect  diatomic  gas  and  0m  = 49.0®  for  nitrogen 
in  dissociation  equilibrium. 


APPENDIX  B 


ANALYTICAL  SOLUTION  OF  THE  EQUATIONS  OF  MOTION  FOR  RR 
(Eqs.  2.1  to  2.9) 


The  following  computer  program  follows  the  method  discussed  in 
Section  2. 2. 1.2  for  solving  Eqs.  2.1  to  2.9  with  and  without  real-gas  effects. 

Although  the  following  program  is  for  nitrogen  it  could  be  modified 
very  easily  for  any  other  gas  by  changing  the  following: 

(1)  SUBROUTINE  RGP  that  calculates  the  equilibrium  thermodynamic  properties 
of  the  gas  under  consideration  for  a given  pressure  and  temperature  - 

(2)  SUBROUTINE  GUESTP  that  supplies  appropriate  initial  guesses  of  tempera- 
ture and  pressure  ratios  across  a normal  shock  wave. 

Two  typical  outputs  follow  the  listing  of  the  program  for  the  two 
possible  modes  of  operation  ( see  comment  statements  in  the  following  computer 
program  for  details) . The  CPU  times  of  these  two  outputs  were  0.19  and  0.06 
sec,  respectively,  on  an  IBM  370  computer. 


B-l 


c— 


■T-t  FCLL--ING  PROGRAM  CALCuLATtS  STaUS  I A I ANw  121  Of  A Kt  UOL Ax 
REFllCUuN.  tvK  A SET  IMI1AC  CwUl'U^,  l.t., 

PO-PRtSSuRt  AHEAD  f Tut  INC  I Dt  N I SmoC*  »A»t.  iTGRKI 
TQ-TfcMPEkAluRE  A«tAU  of  In*  |*ClulM  S»,„CK  -A*t.  tM 
1 1 T A •wCOMMt  **  I UN  wtOGE  *‘iul  L • IUtu«Etl 

“S-lNClDENT  SmuC  a aAvt  MACm  NyMOtR. 

T-0  OPERATION  MODES  AHt  POSSIBLE. 

RHSuL-1.  ONLY  STAtt  til  li  CAlCuLAIEw  .ilAlt  Ul  1 i OlwtN 

DUMMY  .AkbUH 

RRSOL-2.  BOTH  STATES  111  AM,  i>  ARE  C AlCW. * T Ew  . 


CUMMCN/ABC/1  TYPE 
COMMGN/OE/RRSOLV 
REAL  MC.*1  .mj.ms 
INTtGER  RHSOl.RNSOLV 


•RRSOL-i  CALCULATE  STATE  111  UNLY, 

■R  R SOL  ■ 2 CALCULATE  SlATkS  111  AND  HI, 
RRSOL-2 


C— — CONVERTING  FACTOR*. 

C EROM  DEGREE  TO  RADIAN. 

O-O.0l7t.SJJ 

C— — FROM  TORR.  TO  OYN/Cm--* 
00-V8U.LAS-1 . 36 


C— -INITIAL  CONDITIONS  AH  t AS  F ULlUWS 

C po  |N  TORRS. 

P0-1S. Jl 
C— — -TO  IN  K. 

T0-J00. 

C TETA-D  IS  DEGREE. 

TETAwD-fcO. 

MS-*. 68 


WRITE!*. 1161 
116  FORMATIlMlI 

C 1 TYPE-1.  IDEAL  GAS  SOLUTION. 

C- -I  type-2,  real  GAS  solution. 

00  l ITYPE-1,2 
wRITEl*. 11SI 
US  FORMAT!  1*  .801  '•'  I I 


C- -CONVERT  uM ITS  TO  APPROPRIATE  WORKING  UNITS. 

TETAwR-TETAwO-O 
PO-PO-O Q 


C move  FRAME  OF  REFERENCE  TO  ThE  REFLECTION  POINT. 

MO-MS/COS ITETAwR I 
Pn 10-90. -TETAwO 


C— - -CALCULATE  STATE  I0>  FROM  TmE  GIVEN  PO  AND  TO 


CALL  RGPlPO.TO.ALO.ROO.HO.AOI 
U0-A0-M0 


C- — —CALCULATE  STATE  ID. 

CALL  OSwIPO.TO.ALO.ROO.HO.AO.MO.UO.PUT  l.ALl.ROL.Hl  .A.  .M».u*  .PHlU. 
10EL1DI 


C —GIVE  DUMMY  VALUES  TO  STATE  121. 

M2-0. 

P2-0. 

T 2-0. 

R02-0. 

PHllD-O. 

DEL2D-0. 

RSwAD-O. 


GO  T0I2. J I .RRSOL 

C CALCULATE  STATE  I 2 > 

C—  FIRST  GUESS  FOR  PmIID. 

J Pm | 10-60 • 

00  11  I-l.JO 

CALL  0SMIP1.T1  .Al1.R01.m1.A1,M1.u1.P2.T2.AL2.Ru2.h2.A2.M2.U2*Ph|1D 
1.0EL2DI 

F 1-DEL2D/ DEL  10-1 . 

IFIABSIF1 I .LE.l.E-OSIGO  TO  12 
PH|NO-l.01«PHl 10 

CALL  OSw  I PI  . T 1 .AL1  .ROl.Hl  .AuMi.ui  . P2  . T 2 . Al  2 . Ru.  , A.  ,mt  ,\jt  . Pm  I Nu 

1.0EL20I 

F 2-0EL2D/ DEL  10-1 . 

11  Pm| 10-Ph  J 10* II.-0.01-F1/IF2-F1 I » 

WRITE! 6. 999  I 

999  FORMATI 1*. '•••••  NO  CONVERGENCE  IN  MAIN  RPOGHAM  ••••••  I 

STOP 


C— CALCULATE  THE  ANGLE  BETWEEN  TmE  REFLECTED  SHUCK  AND  TmE  wEugE  SURFACE 
12  RSwAO-PhI 10-0EL1D 


C— — CHANGE  PRESSURES  BACK  TO  TORR. 
2 P O-PO/UU 
P1-P1/00 
P2-P2/00 


C -PRINT  RESULTS. 

GO  T0I6.S I .1  TYPE 

6 wR I TE l 6. 1 00  I 

100  F0RMATI2SK. 'PERFECT  GAS  SOLUTION1) 

WRITIIS.101I 

101  F ORMAT I 2Sx  *201 I » 

GO  TO  * 

5 WRITEIS.102I 

102  FORMAT (2SK. 'IMPERFECT  GAS  SOLUTION* » 

WRITEI6. 103) 

103  FORMAT  2SX.22I '-'II 

6 wRl TE I 6. 10S I 

105  FORMAT! 18K. 'REGION  O' • 7X . ' REG! ON  1 • . 7 X . ' REG  I UN  2'l 
wRl TE 16. 1 06  I MO .Ml .M2 

106  FORMATI 1*.' MACH  NUMBE R'.Tx.JiFS.V.lOXll 
WRITE!*. 10TIP0. Pi. P2 

107  FORMAT! IX. 'PRESSURE ' .iOX.JlF 7 . 2 . 8 X I .6 a . • TORR • I 
wRITEI6.10BiT0.Tl.T2 

108  FORMATilx.'TEMPERATuRE'.SX.J!F7.2,v*l • S X • ' K ' I 
wR| T E I 6 . 109IROO.RO1 .R02 

109  FORMATI lx. 'DENSITY* .10X.3IE10.6.SXI . 7X. 'G/CM--J* I 
WRJTEI6. 1101 

110  F Ohm a TI2JX.' I • . 2 jx . ' R ' I 
WRITE!*. 1 11 IPHlO.Pnt 10 

111  FORMATI IX .' INCIDENT  ANGLE • .6* • F 5 . 2 . I9x .F S . 2 • 20X • • DEGREE ' I 

WRITE!*. 112I0EL1D«0CL2D 

112  FORMAT! IX. 'DEFLECTION  ANGLE ' . 6X .F S . 2 . 1 9X .F S . 2 .20X . ' DEGRE E ' I 
WRITE!*.  UJlRSwAD 

111  FORMAT l 1* , 'ANGLE  BETWEEN  R AND  wEDGt  SURFACE  I S ' • IX »F S .2 » 1 X . • UEGRE 

IE  ' I 

1 CONTINUE 


wRITt  16,  U9I 
WRITE!*. 11*1 
STOP 
ENO 


B2 


SUbRUUT I Nt  OSW  <PA • T A .AlA • ROA  »HA » A A (MA ,uA .MB . Tb • Ac o • HUB • MO  * Ao • Mb  • UB 
1«PmIAO*OELBQ) 




C — — SUBRUUT  IN£  OS*  CALCULATES  T Mfc  CHANGE  S ACROSS  AN  ObLlUuE  ShOCa 
C wAVfc  ( OS m > • ThE  MARAVfc.Tfc.Kb  ARE » 

C PA-PRESSURE  AhEAD 

C TA-TtMP£WATu«fc.  AhEAD 

C ALA-DfcGRtt  UF  UlSSOClAT JuN  AntAU 

C ROA-DtNSJTY  AHtAO 

C HA-ENThAlPy  Ah  fc.AU 

C AA-SPfcED  OF  bOuNO  Anfc.AU 

C ma-v»ACh  NUMBER  AHEAD 

C UA— F LUw  VELOCITY  AhEAD 

C PB-PRESSORE  BEHIND 

C Tb- TEMPERATURE  bfc  h J ND 

C ALB-DEGRtE  OF  DlbbuCIATluN  btHjNO 

C ROB-DENSI TY  BEHIND 

C hB-ENThAlMY  bEnlND 

C AB-bPEED  OF  bOuND  BEHIND 

C M8-MACH  NUMBER  BEHIND 

C UB-FLOu  VELOCITY  BEHIND 

C Ph  J AD- I NC  I DENL  E ANGLE  BtTwEEN  SHwCA  AND  FuOw  AHtAD 

C DElBD-FlOw  DEFLECTION  ACROSS. 

c — — 

c Pi  IN  DYN/CM»*2 

REAL  F( 3 » »G< 31 *MA.Mb*MAN 



C— CONVERT  ANGLE  TO  WADjAN. 

(J«0.017<.603  KAD/DtO 

PHI AR-PM I AD*0 



C- -CALCULATE  NORMAL  COMPONENTb  OF  VELOCITY  AN u MACH  NUMBER* 

UAN«uA«b*MPHlAR  > 

MAN«uAN/ A A 

C — — — - — — — 

C- DETERMINE  CONSERVATION  EOuATIuN  CONbTANTb  F ROM  INITIAL  CONwITjUNb 

C l»ROA*uAN 
C2«PA»R0A«uAN*uAN 
C3»hAO.  b *uAN*uAN 



C -BEGIN  SdlvInG  THt  SYSTEM  ObINO  NtwTONS  MtTnuD 

C- — — — CHOObE  STARTING  VALUE  b F vH  T AND  P 
CALL  GUEbTPC MAN»PBA .TBA » 

ANEwP»P8A>PA 

ANEwT»TBA*TA 

c- 

C- — — HtGlN  ITERATION  PRuCtDuRt  » Yflx  Nu  » 2b 

DO  70  I T ER« 1 .2  b 

C- — — DtTERMJNt  NtCtbbARY  PARMA,.  DERIVATIVES  TO  SwlVE  wblNO  Ntw!u..S 

T-ANEwT 

C DEFINE  bTEP  bIZE  F sjR  P IN  DETERMINING  DtRlvATlvES 

STEP»n*10«ANEuP 
DO  60  Km  1 ,2 


C — MOwD  T CONSTANT.  FIND  F AND  G V UK  3 vALOtb  OF  P 

c — — a»2  hold  p Constant,  find  ► and  g fur  i values  of  t 

no  so  *3»n 

GO  TO  ( 32  *34  ) »«. 

32  P-ANEwP* t J-2 >»STEP 
GO  TO  30 

34  T .ANEwT*- « J-2  (•STEP 
36  PB«P 
TR«T 

CALL  RGO » PH. Tb *ALH»ROB»HH*AB ) 

c — — determine  particle  velocity  (shock,  co— oRos i from  Continuity  u 

U-Cl/RCB 

C- —MOMENTUM  ♦ CNIRGV  LUUAHONb  DtF  I Nt  FUNCTION.,  TO  CwtCA  VALUES 

F ( J J ■ ( pb ♦ C 1 •li-C 2 I /C2 
G ( J ) * (0. 5*U*U*Hb-C 3 » ✓ C 3 

C- — — — -IF  F AND  G MEET  ACCURACY  RtOU I KtMLN T S • t A 1 T *ITm  SOLUTION 
I F I AHS<  F (JII-l.E-5>46»48.50 
46  I F ( APS ( j i j l i -1 ,L-b i 6 J .HG.bu 
50  CONTINUE 

C CALCULATE  APPROPRIATF  OEKlVAt IvES 

DFDT-(F(  3>-F  (1)  » / ( bTtP*bTEP  > 

DG0T«(G<  3 l-G(  1 ) I /<  STtP^bTtP  « 

GO  TO  (54*A0I»«.  **  J* 

C ——FOR  i .DERlvAT  IvtS  ARt  *ITh  RtbPtCT  TO  P ▼/%  TA 

54  DFDP-DFDT  •Wjl  ^ 

DG0P-0G0T 

C OEFINE  STEP  SIZt  F JR  T IN  DETERMINING  utKlVATIVtb 

STEP*Q*Oi#ANfcwT 
PaANEwP 
60  CONTINUE 

C— — -DETERMINE  Nt»  VALUES  FOR  P AND  T USING  Nt  w T UNS  METHOD 

»RONb»DFOP*DGDT-DFLT  *DODP 
ANE wP«ANEwP-(F  i 2 »*DGOT-G(2 I«uFDT J/wKUNS 
ANE wT  »ANE  w T « I G ( 2 ) *DF  DP-F  (2  > *DGDP  » /»»RONb 
70  CONTINUE 

C GIVE  ERROR  Mt  Sb AGE  IF  NO  CONVERGENCE 

WRITE (6.2 J 

2 FORMAT ( * *N0  CONVERGENCE*') 

C — — 

C — — -CALCULATE  DYNAMIC  PROPERTIES. 

80  IF (PhJAD-90. I«2»81.81 

81  DELBR-0. 

GO  TO  83 

8 2 DElRR*Rh 1 AR-AT AN ( ROA/ROB*S I N ( PH  I AN  J /COS  < PM 1 AH ) » 

83  0B«C1/SI N (PHI AR-DELBR ) /ROB 
mr«ub/ab 

c- — 

C — CONVERT  ANGLES  TO  DEGREE. 

PMlAD-PMl AR/Q 

OElBD*0ElPR/U 

* — 

RETURN 
t NO 


SUBROUTINE  HOP  I P . I .ALPHA  iKU.H.A  i 
(OMMUN/ABC/ t T»Pt 
HEAL  HO.** 

DIMtNSJUN  t A I * I .EMI  4 1 «OA  I *<  .UMI4  I .BHi  hi  .It  TAM14 I ,04 

1 A>X  2 I 4.  I ,P»1I  } I 


C— — THE  OAi  CONSTANT*  ARt  AS  FuLlO*S 

c — — -k-unuersal  f.A$  const  am 
«»o.#n*.n  os 

C '--ATOMIC  ■ £ I GhT 

“•Ia.oco 

K.R/2./M 

RT*R*T 


OO  TOI  90  • 91  I . I T»PE 

C IDEAL  GAS  SECTION 

90  GAMA*i,4 
AlPhA.O. 

2*1. 

RO*P/RT 

E*1.4  IGA-.A-1  , I ART 
h*{ »RT 

AaSORT IOAmaaVT  l 

he  Torn 


C — *EAl  (.AS  SECTION 

91  CONST »;* J 11 .4 

C CONST*. '•  ( A**  IS  / i I*  J.i<./«»«2/NA)  ••!.* 

c-* — oo.m  sre-is  bolt/^an  constant 

c- M.o.SNiist -26  plan*  constant 

C — NA.J.bCV/SVt  .V.  A.J0A0MO  NUMBER 


c — —ATOMIC  'UNSTANTS 

C—  LA-NU^BER  JF  tNERO*  LEvILS  GIVEN  F UK  Tut  A T OM 

LA*> 

GAI 11*4. 

GA(2I*10. 

GAI 31*10. 

CAItl'ti 

GAISI.6. 

EAIH-0. 
tAI2l*276*6. 
t A ( J I *27856. 

EAI4I-4146S. 

EA(*l*4l4*S. 


C—  MOLECULAR  CONSTANT 

C— -LM-NuMBER  Of  ENtHOv  LEVELS  OIVEN  fOK  THt  MOLtCULE 
LM*2 
GMT  11*1. 

0MI2I-J. 

£*•(21*71*71.466 

EMIll-O. 

BMllJ-2. 661621 
BWT2 1*2.0619964 
OAMAMI 1 I *S. l6*S9*E-06 
0AMAMT21  *i.  74426U-0* 

TETAMili* jj*o.7U*6 
I E T AMI  2 1 *2060.4343 
TETAD-11 JJOO. 

TETAV-3340. 


CONST *C0NST* imi**1.S 
P*2T*0. 

00  1 1*1. LM 
IFIEMTH/T-1T0.T3.J.6 
J BHH*  1 /£  * P 1 1 **  I I I / I I 
OO  TO  > 

4 BHH* 1 /£  XR  < 1 70. i 

* 0T*2i  I 1*1  1 .OAMAMI  | | *T  I *hbB / I l-l/t APT TtTAMI  I I /T I | 

P*2 I I l*GMi | l *uT  « 2 I I i / H * i I i 

PA2T*PA2T*P«2> I I 

1 CONTINUE 
P«1T*0. 

00  2 1*1. LA 

If ((All l/T-170. 16.6.7 

6 CCC*m*PIEAl  1 l/T  I 
GO  TO  e 

7 CCC-l/t API  I TO.  I 
6 P*ll I l«SAI | l*ccc 

Pi  1 T*P« 1 T.P» 1 l 1 l 

2 CONTINUE 

If  I T-6J0.  19.9.10 
9 »0*0. 

GO  TO  U 

10  AAA*E*PI  TET AO/ I4*T I I 

*0*CONST*T ••!.*/ AAA /AAA*P*l ! «PX 1 T /P«2 T / AAA/ AAA 

11  AlPMA.SURT IKD/ If 0*4«P  I l 
2*1. »AL=HA 
R0«P/RT/2 
22-TETAv/T 

Ill’ll/ I EAPl/2 1-1.  i 

t *12* **0.S* ALPHA. ALPHA* Tt  TAO/T.I  1 .-ALPHA  I *222  l*RT 
H*( .P/RO 

CVIBH*E>P<22  I *222 **2 . 

C VR* 1 . *• I 1.. ALPHA |.(1. -ALPHA  I •( I..CV I BR I .ALPHA* ( 1 . -AlPhA I / ( 2 . -Al PH 
lAI*(TETAO/T*O.S-222 

CPR«2. s* I 1..ALPHA I . I 1 .-ALPHA  I • I 1 .*c V | HR  I .ALPHA* I 1.-AlPmA**2  . I /i  . « I 
1T£TA0/T*1.*-222I**2. 

OAMA-CPR/CVR 

GAMAE  «GAMA*2./(2 . .ALPHA— ALPHA*  *2 • I 
OAMAf • ( 7**3. • ALPHA  I / ( *.. AlPHA I 
A ■ SURT I 0 AMA E *2 *RT  I 

RETURN 

END 
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SUBROUTINE  GUEST P ( MAN , PH A . ! B A > 


C -SUBROUTINE  GUESTP  GuSSES  Thi  TEMpRATokE  ANU  PRESSURE  RATIO  ACROSS 

C A NORMAL  SHOCf  KAvE . 

C— — This  SUBROUTINE  IS  fOR  N | TROOt  N • 

C— — The  PARAMETERS  ARE. 

C MAN— f LO«  MACH  NUMBER 

C PBA-PRESSURE  RATIO 

C TB A- TEMPERATURE  RATIO. 


REAL  MAS 


COMMON/ABC/ I TVPt 

c— — first  guess  for  pressure  ratio  across  the  snuca 

01*1.4 

Bl-IGI-1, I/2./G1 
ALFA1*IG1*1.I/(G1  — l.l 
PH A* I MAN*  42 / B1  — 1 , I / AlF  A 1 
GO  TOI 90 S.ll.ITvPfc 


C— — -FIRST  GUESS  »l>«  TEMPERATURE  RATH)  ALRUSS  The  SHUCf 
1 If IPBA.L T.40C. IGC  TO  90C 
TBA*26. 

RETURN 

900  If  IPHA.lT.300. IGO  TO  902 
THA*23. 

RETURN  ■ 

902  If  IPRA.lT.200. IGO  TO  903 
THA-21. 
re  turn 

901  If (P0A.L T.100. IGO  TC  904 
TRA*  1 t • 

RETURN 

904  I f (PRA.l T .40, IGC  TO  90S 
TRA*  7. 

RETURN 

905  TRA.PRA*  I PRA.AlF  A|  I /!  1 . . l-HA  • Ai  » A I I 
RETURN 
I NO 
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1 


RERFECT  GAS  SOLUTION 


MACH  NUMBER 
PRESSURE 
TEMPERATURE 
OENSITY 


REGION  0 
9.36 
16.31 
300.00 
0.2293E-04 


REGION  1 
3.58 
388.95 
1560.00 
0.1120E-03 


REGION  2 
0.0 
0.0 
0.0 

0.0 

R 

0.0 
0.0 
DEGREE 


I 

INCIDENT  ANGLE  30.00 

DEFLECTION  ANGLE  23.26 

ANGLE  BETWEEN  R AND  WEDGE  SURFACE  IS  0.0 


IMPERFECT  GAS  SOLUTION 


MACH  NUMjER 
PRESSURE 
TEMPERATURE 
DENSITY 

INCIDENT  ANGLE 
DEFLECTION  ANGLE 
ANGLE  BETWEEN  R 


REGION  0 
9.36 
15.31 
300.00 
0.2293E-04 
I 

30.00 


REGION  1 
3.81 
396. 77 
1 460 .54 
0.1221E-03 


23.81 

AND  WEDGE  SURFACE  IS  0.0 


REGION  2 
0.0 
0.0 
0.0 

0.0 

R 

0.0 

0.0 

DEGREE 


TORR 

K 

G/CM**3 

DEGREE 

DEGREE 


TORR 

K. 

G/CM**3 

DEGREE 

DEGREE 


RERFECT 

GAS  SOLUTION 

REGION  0 

REGION  1 

REGION  2 

MACH  NUMBER 

9.36 

3.58 

2.13 

PRESSURE 

15.31 

388.95 

2147.41 

TORR 

TEMPERATURE 

300.00 

1560.00 

2907.71 

K 

DENSITY 

0.2293E-04 

0.1120E-03 

0.3319E-03 

G/CM**3 

INCIDENT  ANGLE 

30.00 

38.08 

DEGREE 

DEFLECTION  ANGLE 

23.26 

23.26 

DEGREE 

ANGLE  BETWEEN  r 

AND  WEDGE  SURFACE  IS  14.82 

DEGREE 

IMPERFECT  GAS  SOLUTION 

REGION  0 

REGION  1 

REGION  2 

MACH  NUMBER 

9.36 

3.81 

2.42 

PRESSURE 

15.31 

396.77 

2253.20 

TORR 

TEMPERATURE 

300.00 

1460.54 

2498.72 

K 

DENSITY 

0.2293E-04 

1 

30.00 

0.1221E-03 

0.4052E-03 

G/CM**3 

INCIDENT  ANGLE 

36.28 

DEGREE 

OEFLECT I ON  ANGLE 

23.81 

23.81 

DEGREE 

ANGLE  BETWEEN  R 

AND  WEDGE  SURFACE  IS  12.47 

DEGREE 
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APPENDIX  C 


ANALYTICAL  SOLUTION  OF  THE  EQUATIONS  OF  MOTION  FOR  SMR 


[Eqs . 2.17  to  2.30) 


The  following  computer  program  follows  the  method  discussed  in 
Section  2. 3.1.3  for  solving  Eqs.  2.17  to  2.30  with  and  without  real-gas  effects. 


Although  the  following  program  is  for  argon  it  could  be  modified  very 
easily  to  fit  any  gas  by  changing  the  following: 

(1)  SUBROUTINE  RGP  that  calculates  the  equilibrium  thermodynamic  properties 
of  the  gas  under  consideration  for  a given  pressure  and  temperature. 

(2)  SUBROUTINE  GUESTP  that  supplies  some  initial  guesses  of  temperature  and 
pressure  ratios  across  a normal  shock  wave. 

Note  that  unlike  SUBROUTINE  RGP  of  Appendix  B where  the  thermodynamic 
properties  are  calculated  from  the  partition  function  P(T)  itself,  here,  a fit 
of  the  form  P(T)  = f(T)  was  used  in  order  to  reduce  the  computing  time  and  hence 
operational  cost. 

A typical  output  follows  the  listing.  The  computing  time  for  this 
typical  output  was  3*6  seconds.  The  author  is  indebted  to  Mr.  B.  T.  Whitten, 
Research  Assistant  at  UTIAS,  for  svpplying  him  with  subroutines  RGP  (for  argon) 
and  MCNPS . 


C2 


SUBRUUT  1 Nt  OSw  (PA.TA.ALA.ROA.hA.AA.MA.UA  .Pd.  Ib.ALtt.WUb.Mb.AD.MD.UO 
1 .PHI AD.DELBD ) 


c— SUBROUTINE  OSW  CALCULATES  Th£  CHANGES  ACROSS  AN  OBLlUUE  SHOCK 

C WAVE  (OSw)*  TML  PARAMETERS  ARE. 

C PA-PRESSORL  AHEAD 

C TA-TEMPtRATURE  AHEAD 

C ALA-DEOREE  OF  DISSOCIATION  AHEAD 

C ROA-UENSITY  AHEAD 

C HA-ENTHALPY  AHEAD 

C AA-SPEED  OF  SOUND  AHEAD 

C MA-MACH  NUMBER  AHEAD 

C UA-FLOW  VELOCITY  AHEAD 

C PB-PWESSURt  BEHIND 

C TU-TEMPERATURE  BEHIND 

C ALB-DEGREE  OF  DISSOCIATION  BtHlND 

C RO0-DENSITY  BEHIND 

C HB-ENTHALPY  BEHIND 

C AB-SPEED  OF  SOUND  BEHIND 

C MB-MACH  NUMBER  BEHIND 

C UB-FLOW  VELOCITY  BEHIND 

c phi  ad-incidence  angle  between  shuck  and  flow  ahead 

C DELBD-FLOW  DEFLECTION  ACROSS* 


C pi  I*  DYN/CM*»2 

REAL  F( 3 ) *G( 3)  .MA.MB.MAN 

0*0*0174603  RAD/DEO 

PHI AR*PH I AD*0 

UAN*UA*SIN(PHIAR> 

MAN-UAN/AA 


C DETERMINE  CONSERVATION  EQUATION  CONSTANTS  FROM  INITIAL  CONDITIONS 

C 1 *RQA*U AN 
C2*PA4-R0A»UAN*UAN 
C3*HA*0.  5*UAN*UAN 

C- ——BEGIN  SOLVING  THE  SYSTEM  USING  NEWTONS  METHOD 

C— — — —CHOOSE  STARTING  VALUES  FOR  T AND  P 
CALL  GUESTP(MAN.PBA*TBA) 

ANEWP*P0 A*PA 
ANEWT*TB A*TA 

C — —BEGIN  ITERATION  PROCEDURE.  MAX  NU  ■ 26 

DO  70  ITER*1.25 

C ——DETERMINE  NECESSARY  PARTIAL  DERIVATIVES  TO  SOLVE  USING  NEWTONS 

T*ANEWT 

C — — -DEFINE  STEP  SIZE  FOR  P IN  DETERMINING  DERIVATIVES 

STEP*0*10«ANEWP 
DO  60  K* 1*2 

C- — — — K*1  HOLD  T CONSTANT.  FIND  F AND  G FOR  3 VALUES  OF  P 

C K*2  HOLD  P CONSTANT*  FIND  F AND  G FOR  2 VALUES  OF  T 


DO  SO  J*  1 » 3 » K 
GO  TO  ( 3 2 * 34  ) .K 
32  P*ANEWP+ ( J-2 >*STEP 
GO  TO  36 


34  T*ANEWT*< J-2)*STEP 
36  PB*P 
TB*T 

CALL  RGP( PB.TB.ALB.ROB.HB.AB) 

C— — — DETERMINE  PARTICLE  VELOCITY  (SHOCK  CO-ORDS ) FROM  CONTINUITY  E(J 
U*C1/RO0 

C— — ——MOMENTUM  ♦ ENERGY  EUUATIUNS  DfcFlNt  FUNCTIONS  TO  CHECK  VALUES 
F ( J)*(PB+C1*U— C2 ) /C2 
G ( J ) ■ ( 0. 5*U*U*HB— C3 ) /C3 

C— — — -IF  F AND  G MEET  ACCURACY  REQUIREMENTS.  EXIT  WITH  SOLUTION 
IF(ABS(F( J) )—l*E— 5)46.40.50 
48  I F ( ABS (G(J) )-l*E-5)80.80.50 
50  CONTINUE 

c— — CALCULATE  appropriate  derivatives 

DFDT* ( F ( 3 )-F ( 1 ) ) / < $TEP*STEP ) 

DGDT*(G( 3 ) — G ( 1 ) ) / ( STtP*$TEP ) 

GO  TO  ( 54  .60  ) . K 

C — FOR  K*l. DERIVATIVES  ARE  WITH  RESPECT  TO  P 

4 dfdp*dfdt 
dgdp*dgdt 

C- —DEFINE  STEP  SIZE  FOR  T IN  DETERMINING  DERIVATIVES 

STEP*0*01*ANEWT 
P*ANEwP 
60  CONTINUE 

C— — — DETERMINE  NEW  VALUES  FOR  P AND  T USING  NEWTONS  METHOD 
WRONS*DF  DP*DGDT-DFDT*OGDP 
ANEwP*ANEWP-( F ( 2 ) *DGDT -G ( 2 ) *DFDT ) /WRONS 
ANEwT*ANEWT-(G ( 2 ) *DFDP-F ( 2 ) *DGDP ) /WRUNS 
70  CONTINUE 

C GIVE  ERROR  MESSAGE  IF  NO  CONVERGENCE 


WRI TE ( 6. 2 ) 

2 FORMAT ( ' *N0  CONVERGENCE*') 

80  I F ( PHI AD-90 • >62.81 .81 

81  DELBR-O. 

GO  TO  83 

82  DELBR*PH I AR-AT  AN ( ROA/ROB*S I N ( PH  I AR ) /COS ( PH  I AR ) ) 

83  UB-C1/SI N (PHIAR-DELBR ) /ROB 
MB*UB/AB 

PHI  AD* PH  I AR /O 


DELBD*DELBR/Q 

PBA-PB/PA 

tba*tb/ta 

RETURN 

END 
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SUBROUTINE  RGPIP.T .ALPHA ,R0 »H » A ) 
HEAL  CA(  lwi  ,d  ( 14)  .MA.NA.Nfc 
COMMON/ABC/ I TYPE 


1 


— ARGON  GAS  CONSTANTS 

—DETERMINE  EQUILIBRIUM  PROPERTIES  of  THE  MQNATOMC  GAS.  given  t a no  p 
—ENTER  FITTING  COEFFICIENTS  OF  PARTITION  FuNCTS  VS  TEMP  FOR  ATOMS  AND  IONS 


DATA  CA  /ll..  0.101*73*  01.  0.2*W31E-ul  • V.  Ik /31OE-01 . AH 

1 0.S122*.bE-02»  0.9S0966E-0*.  -0 . 6 i 3 330E-U  3 . -0 . 2009  3 7E-0  3 . AH 

2 0.2B7S63E-03.  0 . 32  76SOE-U3  • 0. 1 B360SE-V3 . 0. 662  *B2E-0*  • Ah 

3 0.106J82E-03*  -0.1319UE  01/  AH 

DATA  Cl  /11m  0.SS7626E  01.  0.3*7bU0E  Oi.  -0.1333SSE  00.  AH* 

1 0.S66B9SE-01  . -0  . i J6u*«t  - J 1 . 0.fl67t><.3t-02  * -0.  32  1 b 1 *E-02  • AH* 

2 0.122572E-02.  -U.*67ulbE-u3 . 0«  161S1  bE-L3 . -U.b 773b2E-0* • AH* 

3 0.106382E-03.  -0.1319UE  01/  AR* 


VA-6.633AVE-23 
W »2 .OH  12  9t  6 
t i-ia/av j. 

K T *H*T 

GO  TUI  1,2  I .1 TYPt 


— IDEAL  GAS  CONSIDEhAT  luN 
1 GAVA.S./3. 

ALPMA»0. 

2-1. 

RO-P/RT 

E»1./(GAMA-1.)*KT 

h*E*RT 

ASSORT (GAMA»RT I 
RETURN 


— REAL  GAS  consideration 

2 ZA-l.O 
DZA*0.C 

IF  IT.GT. 1UOOO. )CALL  MCNP&t 2 A »D2A ,1 .LA) 

CALL  MCNPSI2 I »D2 I ,T ,C 1 ) 

D»0.0 

IF l T.OT.2S00. IU*0.66to7l<.b*Z  l*T**/.b  / UA*P*E*PI  T 1/ T) I 
ALPhA-SQRT (D/I l.*DI » 


R0»P/(R*T«( 1,*ALPmA) ) G/CMj 

NE*R0*ALPhA/MA  /CM3 

NA*RO* ( 1 .-ALPHA ) /MA  /LM3 

M»R*  ( 2 . b*  I 1 .0*  ALPHA  ) «T*T  J*ALPmA*T*T*I  ( I.u-AlPhA  ) *DZA/2  A* ALPHA *uZ  I / 

121))  ERGS/G 

DALOT»0. b*( 1 .-ALPHA* ALPHA )« ALPHA* l2.b*Tl/I*T*(UZl /ZI-DZA/ZA  * » / T 

NEGLECT  EXCITED  STATE  CONTRIBUTION 

CPR-R* I 2 . b«T*T I ) *DALOT  EWGS/G-A 

CP«2.b*R» I 1 .0* ALPHA ) *CPR  EHGS/u-a 

A "SORT ( ( b .*AlPhA* ( 1 .-ALPHA )*«2.b*TI/T)**2l*R*T/( 3 • -1 • b« ALPHA* 

1 ALPHA •( 1 .-ALPHA ) / ( 1 .*ALPHA ) • ( 1 . b* T I / T ) • •* ) I CM/ SEC 

RETURN 


E'.L' 


SUBROUTINE  VCNPSIV.OY.X.C  I 
DIMENSION  C( I A ) 

V * 0 * 0 
DY*0 .0 
N«C  I 1 I 

I F ( N— 1 1 J , 2 . 1 

1 A«G*C  r»*2  ««K*C  < N*  3 1 

1 F I AbS IAHGI— 2.000)6.6.3 

2 v *C I 2 ) 

3 RETURN 

6 ARG-ARG*  ARC. 

H 1 *0 , 0 
M0*0.0 

OMl *0.0 
UHO*0.0 
DO  7 1*1. N 
A *N* 1 — I 
"?*Hl 
H 1 »M0 
0H2-0H1 

DmI-umO 

0Hr*ARG*D"l-DH2*2»Hl 

7 M0*ARG«Hl-H2*C IA*1 ) 
Y*0.b*lC l 2 I — h2*hO ) 

DY*j.b»C I N *2 )• IUHU-0H2  > 
RETURN 

END 


SUBROUTINE  GUESTP(PAN.PPA.TBA) 


C— —SUBROUT  1 NE  GuESTP  GuSSES  The  TEMPRATuRE  ANU  PRESSURt  RATIO  ACROSS 

c A NORMAL  SHOCK  -ave. 

C- -THIS  SUBROUTINE  IS  FOR  ARGON. 

C— — The  PARAMETERS  are. 

C man-FlOw  MACH  NUMBER 

C P9A-PRESSURE  RATIO 

C TBA-TEMPERATURE  RATIO. 


REAL  man 

COMMON/A0C/ I TYPE 

C FIRST  GUESS  FOR  PRESSURE  RATIO  ACROSS  THE  SmOCa 

G1-S./3. 

B l* (Gl-1 . 1/2. /G1 
ALFA1*IG1*1.  I/IG1-1. ) 

PBA* (MAN* *2 /Bl— 1 . ) /ALFA  1 


GO  TOI 800 .1 ) .ITYPE 

C- FIRST  GUESS  FOR  TEMPERATURE  RATIO  ACROSS  THE  SHOCK 

1 lF(PBA-i.ASieOQ»6QO*b01 

800  TBA-PBA*  t PBA*ALFA1 ) /( l.*PBA*ALFAl I 
RETURN 

801  TBA*PBA*( PBA*ALFAi) /( l.*PRA*ALFAl !•( 3./(MAN*A,| ) 
RETURN 

END 
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APPENDIX  D 


BffT.AIIVE  MPT  ION  OF  THE  SECOND  TRIPLE  POINT  (OR  KINK) 
WRT  THE  FIRST  TRIPLE  POINT 


The  self- similarity  (pseudo- steady)  of  the  reflection  phenomenon, 
i.e.,  the  fact  that  any  point  on  the  wave  configuration  having  a radius  vector 
f with  the  corner  as  origin  is  transformed  to  a new  point  C?  where  C is  a 
scalar  constant  implies  that  the  angles  between  the  shock  waves  do  not  change. 
Consequently,  the  second  triple  point  Tq  (or  kink  K)  moves  along  a line  formed 
by  the  straight  portion  of  the  reflected  shock  wave  (i.e.,  the  portion  between 
T and  Ti  or  K)  wrt  the  first  triple  point  T. 

Consider  Fig.  D.la  in  which  the  triple  point  T moves  with  the  velocity 
VTL  along  a line  inclined  at  an  angle  x to  the  wedge  (wrt  the  laboratory  frame 
of  reference) . 

The  second  triple  point  Ti  moves  along  the  line  defined_by  the 
straight  portion  of  the  reflected  shock  wave  R with  the  velocity  Vit  (i  stands 
for  Tq  or  K, respectively)  wrt  the  triple  point  T.  In  the  laboratory  frame  of 
reference  Tq  moves  along  a line  inclined  at  an  angle  x'  to  the  wedge  surface 
with  a velocity  VjLL- 

Single  vector  analysis  yields: 

|VTL  l = UsCOsec<t>0  (D-1) 


(Note  |VTL|  s UQ  by  definition.) 


If  the  distance  travelled  by  the  incident  shock  wave  I from  the 
moment  it  collided  with  the  wedge  corner,  t = 0 (the  moment  when  the  reflec- 
tion started)  is  Ls,  then  the  time  passed  is: 


At 


(D.2) 


and  the  velocity  of  T]_  wrt  T is: 


(D-3) 


where  Li  is  the  length  of  the  straight  portion  of  the  reflected  shock  wave  R, 
and  L is  defined  as  Li/Ls. 

The  velocity  of  Ti  in  the  laboratory  frame  of  reference  is  then: 


(D.4) 


D-l 


The  x and  y canqponents  of  Vjx  can  be  expressed  as  follows: 


V^L  - |VTL  |sin0o  - |»  jsi.nu 


(D-5) 


= I^TL  lCOS0o  + t^iT  lCOSW 


(D.6) 


where  u is  the  angle  between  the  incident  and  reflected  shock  waves  (I  and  R). 

The  velocity  of  T^  in  the  laboratory  frame  of  reference  now  becomes: 

P1Ll  - tO'J/  + (^l)2i1/2 

With  the  help  of  Eqs.  D-3  and  D.4  the  last  equation  can  be  rewritten 
as: 

|Vii  | = Us[cosec20o  + 2Lcosec0ocos(0o  + u)  + L2]1^2  (D.7) 

If  the  frame  of  reference  is  attached  to  Ti  the  angle  of  incidence  between 
the  flow  in  state  (0)  and  the  incident  shock  wave  -0^  is: 

V* 

0'  = tan-1  — 

° v* 

iL 

and  again  by  inserting  Eqs.  D.3  and  D.4  one  obtains: 


0'  = tan 
o 


-1  T 1 - Li 
cot 0O  + 


- Lsinu 
i + Lcosu 


(D.8) 


The  angle  u can  be  expressed  using  Fig.  D.lb  as  follows: 


w = l80°  - (01  + 0Q  - e1) 


(D.9) 


Thus  Eqs.  D.7  and  D.8  can  be  rewritten  in  the  form: 


|Vii  | = Ug[cosec20o  + 2Lcosec0^cos(0^  - 01)  + L2]1//2  (D.10) 


■if  1 " Lsin^l  + ' 91^ 


0o  tan  cot0Q  - Lcos(0^ 


*o  - V • 

+ *o  “ el)  . 


( D.ll) 


D-2 


It  should  be  mentioned  here  that  states  (l) , (2)  and  (3)  are  found 
by  solving  Eqs.  2.17  to  2.30  which  holds  here  too.  Consequently,  all  the 
parameters  (except  L which  should  be  measured  from  an  actual  interferogram) 
are  known.  Note  that  one  could  measure  (from  an  actual  interferogram)  the  value 
of  x'  [the  second  triple  point  (or  kink)  trajectory  angle]  rather  than  L (=Li/Ls, 
Fig.  D.la)  and  perhaps  simplify  the  above  equations.  However,  the  accuracy  in 
measuring  angles  is  in  general  much  smaller  than  that  of  lengths,  and  hence  the 
present  method  ensures  a greater  accuracy. 

It  is  worthwhile  mentioning  that  Law  and  Glass  (see  Ref.  51  for  details) 
developed  an  analytical  method  by  which  the  location  of  Ti  can  be  predicted  in 
fairly  good  agreement  with  experiments  in  the  range  Qy,  < 40°(Ref . 62) . Applying 
their  method  will  result  in: 


L = -2  cosec(01  + 0Q  - 9^ 


(D.12) 


Inserting  Eq.  D.12  into  Eqs.  D.10  and  D.ll  yields: 


IVnl  = us  cosec 20q  - 2 cosec^  + 0Q  - e1)cosec0ocos(01 

/ p0  \2  2 1 1/2 
+ — J cosec  (0X  + 0Q  - 01) 


(D.13) 


0'  = tan 
*o 


1 ~ Pp/Pl 

po 

CO t0o  - — COt(01  + 0O  - 0^ 


(d.i4) 


Equation  D.13  can  be  rewritten  in  a slightly  different  form: 


M'  = M 
o o 


^ Po  sin0ocos(01  - 0^ 

1 " 2 sin(01  + 0Q  - 0^  + 


po  \2 


si n20 


1 / sin  (0L  + 0Q  - 01) 


(D-15) 


where  Mq  is  the  Mach  nunfcer  of  the  flow  in  state  (0)  wrt  the  triple  point  T^. 

The  Mach  number  of  this  flow  wrt  the  incident  shock  wave  in  a frame  of  reference 
attached  to  is: 


M 

m+  = s 
o sin0^ 


Using  the  relation  Ms  = Mosin0o  finally  results  in: 


+ sin^o 

Mo  = Mo 


(D.16) 


i 


ik 


It  is  worth  mentioning  that  all  the  variables  involved  in  Eqs.  D.13  to  D.l6 
consist  either  of  initial  conditions  Us,  Mo,  0o  811(1  Po  or  variables  that 
result  from  the  solution  of  Eqs.  2.17  to  2.30,  i.e.,  0i,  0i,  pq.  Consequently 
|ViL | > 0 o>  Mi  and  M+  can  all  be  calculated  once  the  first  triple  point  is  solved. 
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APPENDIX  E 


METHOD  OF  REDUCING  ISOFYCNICS  FROM  I NT  ERF  PROGRAMS 


The  present  scheme  for  analyzing  interferograms  was  developed  by 
Whitten  and  reported  in  detail  in  Refs.  71  and  75-  It  requires  the  transforma- 
tion of  the  fringe's  location  into  digital  form  wrt  a chosen  reference  point. 

The  digitized  results  of  two  interferograms,  flow  and  no- flow,  are  shown  in 
Figs.  E.la  and  E.lb,  respectively.  The  same  physical  point  is  chosen  for  all 
the  interferograms  (of  cne  experiment)  as  the  reference  point  and  it  is  assigned 
to  be  the  origin  (0.0);  see  Figs.  E.la  and  E.lb.  Each  fringe  can  then  be  thought 
as  a line  [or  locus  of  (x,y) -point]  having  the  value  of  fringe  number  (inter- 
ference order).  With  all  four  interferograms  (1  flow  and  1 no- flow  for  each 
wavelength)  of  one  experiment  digitized  in  this  fashion  wrt  the  same  reference 
point  (origin),  the  actual  fringe  shift,  i.e.,  change  in  the  value  of  the  fringe, 
at  any  (x,y) -coordinate,  can  be  easily  determined. 

To  transform  the  fringes  into  digital  form,  a digitizer  was  used. 

First,  it  involved  picking  up  a common  reference  point  in  all  four  interferograms 
(usually  a cross  hair  that  was  put  in  the  field  of  view  prior  to  the  experiment) 
and  then  a hand-held  cursor  was  used  to  trace  the  fringes,  and  their  location 
wrt  the  reference  point  were  recorded  on  a magnetic  tape.  The  recorded  data 
were  later  transferred  onto  computer  cards.  An  IBM  1800  plotter  was  then  used 
to  reproduce  the  digitized  version  of  the  interferograms  (Figs.  E.la  and  E.lb). 

By  comparing  the  digitized  and  original  interferograms,  all  the  points  misaligned 
due  to  errors  in  the  digitizing  process  were  corrected  and  rechecked  until  a 
satisfactory  reproduction  was  obtained.  The  original  interferograms  were 
enlarged  to  about  four  times  their  actual  size  in  order  to  make  the  whole 
digitizing  correction  procedures  much  easier  and  more  accurate. 

Once  the  reproduction  was  satisfactory,  all  the  points  were  linearly 
interpolated  into  uniform  steps.  Fringe  numbers  were  also  interpolated  and 
assigned  to  each  (x,y) -coordinate.  With  the  flow  and  no- flow  pictures  of  the 
two  wavelengths  digitized  and  fringe  numbers  assigned,  fringe  shift  at  any 
point  could  be  calculated  yielding  in  turn  the  changes  in  flow  properties 
(density  and  degree  of  dissociation)  at  that  point.  Consequently,  the  entire 
density  field  was  then  known. 

Once  the  above  procedure  was  complete,  the  density  distribution  along 
lines  of  constant  x or  y was  found.  Typical  density,  distributions  (i.e.,  p/pn 
vs  y where  pn  is  a normalizing  density)  along  two  different  constant-x  lines, 
x = 0.00  and  -5-00  (on  Figs.  E.la  and  E.lb)  are  shown  in  Figs.  E.2a  and  E.2b 
respectively.  By  changing  x in  small  steps  (0.5  - 2 mm  usually)  a set  of 
these  lines  of  density  distribution  was  obtained.  From  these  density-distribu- 
tion figures,  (x,y) -coordinates  that  correspond  to  the  same  p/pjj  were  deduced. 

For  example,  Figs.  E.2a  and  E.2b  suggest  that  p/pn  = 0.9  at  x = 0.00,  y=-  -32.42 
and  x = -5.00,  y = -31.12.  Sets  of1 (x,y) -coordinates  for  different  values  of 
p/pn  were  obtained  from  the  set  of  density- distribution  figures  (Figs.  E.2a  and 
E.2b,  for  example)  in  a similar  way  and  finally  the  isopycnics  were  drawn.  Note 
that  pu  ^ pQ  and  hence  p/pN  was  multiplied  by  pj i/po  in  order  to  get  the  isopyc- 
nics in  the  form  p/p0- 


E-l 


In  selecting  the  values  of  p/pQ  (i.e.,  p/pn)  for  which  sets  of 
(x,y) -coordinates  were  obtained,  care  was  taken  that  the  step  Ap/po  will  be 
greater  than  the  maximum  possible  error  (Eqs.  8.1  and  8.2)  associated  with 
the  reduction  of  density  from  interf erograms . 


Density  Distribution  Along  Constant  X - P/ PH  Density  Distribution  Along  Constant  X - P/ PH 


Location  of  the  Wedge  Surface 


X * 0.00  mm 


^ Jump  Across  the 
\J  Reflected  Shock  Wave 


PN  - §* 


I Uniform  Region  Behind 
! Incident  Shock  Wave 


065 1 

-50 


Distance  Along  Y (mm) 


FIG.  B-2(a) . DENSITY  DISTRIBUTION  ALONG  X * 0.00  OF  FIG.  E-l. 
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FIG.  E-2(b) . DENSITY  DISTRIBUTION  ALONG  X - - 5.00  OF  FIG.  E-l. 
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